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IV.  SHALLOW  SHELL  ANALYSIS  OF  GRAVITY  LOAD 


A.  SHALLOW  SHELL  SOLUTION 

A  paraboloidal  shell  of  revolution  is  considered  shallow  if 

rdy3\2  2  , 

7  «  1 


so  that 


1  +  v 


1 


(IV-1) 


( IV  —  2) 


If  terms  involving  transverse  shears  are  neglected  in  the  tangential  equilibrium  equations  and 
terms  involving  tangential  displacements  are  neglected  in  the  expression  for  the  rotations,  a 
consistent  application  of  (IV-2)  leads  to  the  following  two  simultaneous  linear  partial  differential 
equations 


AV^V2$  =  —  v^w  +  A(1  —  v)  v^P 


dv2v2w  =  Yi  v2*  -  j-  p  +  p^ 


where 


92(  )  1  9(  )  1  92(  ) 

V  ()  ~  — ?  +7  -9F  +  “2  — 2 

3r  r  3© 


and  P  is  a  load  potential  such  that 
3P 

3r  ^r 

1  3 P  = 
r  90  PG 

F^or  gravity  load 


(IV -3) 
(IV-4) 


(IV-5) 
( IV  —  6) 


P  =  2fph^y  sin©  sin^  —  Xg-  cos^j 


<i>  is  a  stress  function  from  which  the  membrane  stress  resultants  are  derivable 

1  34>1 


N  - 
r 


Ne  = 


i  p 

r2  90 2  r  9r 


2 

a  * 


3r 


N~  -  I  *  9* 


r© 


1  924> 


2  90  r  9r90 


( IV  —  7 ) 


(IV -8) 


(IV -9) 


(IV-10) 


The  transverse  shears  and  moment  resultants  are  related  to  the  transverse  deflection  w  by 

.2  / ,  „  J 


M  =  -  D 
r 


3  w  .  / 1  3w  ,  1  3  w\l 

1?  +  *V  a?  +  7  ^)1 


(IV-ll) 
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Fig.  13.  Master  flow  chart. 
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(IV- 12) 


a2w  ,  . 

fi 

3w  .  1 

92w\1 

v  T  1 
9r 

u 

SF  +  “2 

r 

ae2'l 

Mre--D(l-,)  £  (ig) 

Qr—U  g|(V2w) 

%=-d7^(v2w)  • 

To  be  consistent  with  the  shallow  shell  approximation, 
follows  [see  Eqs.  (11-28)  and  (11-29)]: 


( IV  — 13) 
( IV- 14) 
(IV-15) 

the  rotations  must  be  approximated  as 


Ow 

wr  Or 


(IV- 1  6) 


1  Ow 
r  oe 


(IV-17) 


Having  4>  and  w,  the  two  tangential  middle-surfaee-displaeement  components  u  and  v  ean 
be  obtained  by  integrating  the  strain-displaeement  relations  for  a  shallow  shell  corresponding 
to  (11-25),  (11-26),  and  (11-27). 

The  solution  to  (IV-3)  and  (IY-4)  involves  Thomson  functions  of  order  zero  and  their  deriv¬ 
atives,  W  hile  these  functions  have  been  tabulated  in  various  sources  (most  extensively  in  Ref.  6), 
none  of  these  tabulations  is  complete  for  our  purpose.  For  any  particular  shell,  we  may  need 
values  of  these  functions  not  tabulated  in  the  given  reference.  Subroutines  have  been  written  in 
double  precision  arithmetic  to  generate  the  values  of  these  functions  for  any  argument  (within 
the  limitations  of  the  IBM  7094).  Because  of  the  double  precision  arithmetic  and  the  nature  of 
the  computations  involved,  these  subroutines  are  extremely  time  consuming.  On  the  other  hand, 
the  nature  of  the  Thomson  functions  is  such  that  they  contribute  significantly  to  the  complete 
solution  only  in  the  neighborhood  of  the  edges  of  the  shell.  In  these  regions,  the  magnitude  of 
these  functions  varies  rapidly  over  a  short  interval  of  the  independent  variable.  To  get  sufficient 
information  about  the  shell  behavior,  it  is  necessary  to  have  smaller  increments  between  the 
y-tabular  points  in  these  regions.  In  contrast,  the  solution  in  the  interior  varies  slowly  over 
the  span  of  the  shell,  so  that  we  do  not  need  the  same  density  of  y-tabular  points.  A  waste  of 
computer  time  can  and  will  be  avoided  by  changing  the  increment  of  the  y-tabular  points,  once 
w°  move  away  from  the  edge. 


B,  DIGITAL  COMPUTER  PROGRAM 


The  general  scheme  of  the  program  for  a  shallow  shell  analysis  is  outlined  in  the 
flow  chart  (Fig.  13).  When  it  fails  to  locate  additional  input  at  step  (1)  or  step  (8),  the 
exits  automatically.  The  program  listing  is  presented  in  Table  V. 


master 

program 
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TABLE  V 

PROGRAM  LISTING 


*  SHEA  ,  ELA I  NE  SHALLOW  SHELL  WITH  GRAVITY  LOAD 

*  XEQ 

C  SHEA  *  FLA  I NE  SHALLOW  SHELL  WITH  GRAVITY  LOAD 

C 

c 

c 

DIMENSION  R( 50)  .  Y( 50 ) .BPSRI 50) .bBEI C  50) , AAKER ( 50  I  *  AAK E I C  50  I .XI 2  I 
1 ,  AAKERPI  50)  .AAK.EIPI  50) ,BBRP( 50) ,Bt>IP( 50) ,SBER( 50) .S3EI ( 50) ,SKER( 50 

2  )  ,SkEI  (  50  )  *  TBER  (  50  )  .  TBEI  (  50)  .  T'.ERI  50  )  .TICE  I  (  50  )  ,WW>  50  I  ,W3(  50)  . 

3  W ( 50, 10 ) ,WP(  50*  10  )*  VI  50*10)tll(50*10)  *CMR(  60*10)  *CMT  (  50.10  I  . 

4CMRT (50.10) . CNR ( 50,10’ *CNT( 50* 10) *CNRT ( 50 . 1 0 ) »QR t  50 » 10 t ,QT I  50 . 10 ) . 
5VV(  60) »UU  ( 50) » Z  3 ( 8 ) *  DMR ( 50) »DMT ( 50  » , DMRT ( 50 ) *DNR ( 50 ) *DNT I  50 1 .DNRT I 
650) ,OQR (50) » QQT ( 50 ) *P ( 50 ) . W1 ( 50 ) . U ( 50 ) *EMR( 50) *EMT (501 ,  ENR  t  501 .ENT 
7(5C),Q1(50),A(8,8)»S(8),AA(4,4),SS(4)»22(4I,B(6,6|,BB(6),B9(6),D(3 
8 .3 ) ,DD( 3 ) *D3( 3 ) ,RG( 50 ' 

DIMENSION  P2( 50) *W2( 50 ) *FMR( 50) *FNR( 50) .02(50) *FMT (50) *FNT ( 50) .03 
1(50) »W4(50)»W5(50) » GMR (50) .GMT (50) ,GMRT ( 50 ) *0RR ( 50) .OTTI50) » GNR ( 50 
2 ) ,GNT ( 50) ,GNRT( 50 ) ,U2 ( 50 ) . V2 ( 50 ) .WT ( 50 ) .WFLEXT (50) ,WPT ( 50 ) *UT t 50  I 
DIMENSION  THETA! 20) .THETRI20) .FMT1 ( 12 ) ,FMT2( 12 ) 

DIMENSION  PWB( 50 ) ,WA5 ( 50) .WFLEX ( 50, 10 ) ,PWB2( 50) ,WA54 ( 50) 

COMMON  W.WP.V.U.CMR.CMT ,CMRT .CNR.CNRT .CNT 

EQUIVALENCE  (  BBER.P2  )  ,  (  BBt  I  ,  W2  )  ,  (  AAKER.FMR  )  ,  (  AAK.E  I  ,  FNR  )  ,  (  A  AK.ERP  • 
102) » (AAKEIP.FMT ) , ( BBRP.FNT ) , ( BB I P ,03 ) , ( SBER , W4 ) , ( SBE I . W5 I . ( SKER . 
2GMR ) , ( SKEI , GMT ) » ( TBER.GMRT ) , ( T BE  I . ORR I , ( TK ER » OT T ) . ( TKE I , GNR ) 

C 

C  -  -  (1)  INPUT  I  (CONTROL. GEOMETRY, MATERIAL. EDGE  CONSTRAINT) 

C 

98  READ  INPUT  TAPE  2 , 101 . M , MU , ML , MC 

READ  INPUT  TAPE  2 » 101 . N , NORM , K1 , K2 , K 3 
READ  INPUT  TAPE  2.102.RHO.E.PRT 
READ  INPUT  TAPE  2 , 103 , R 1  .R2  .R3 , R4 
READ  INPUT  TAPE  2 , 103 . THE T A3 , THE T A4  ,  F  ,H 
IFIN1182.182.181 

131  READ  INPUT  TAPE  2 . 205 . FMT 1 . FMT 2 
182  WRITE  OUTPUT  TAPE  3,10 
C 

C  -  -  (2)  GENERATE  THE  INDEPENDENT  VARIABLES 
C 

F2*  2 • *F 
F1-1./F2 

C»SORTF( 12. *( 1 • -PRT**2  )  ) 

C-SORTF(C) 

C*C/SQRTF( F2*H) 

PI-3.1415926535 

RAD-57.2957795 

R5-0. 

36*0  • 

IF!MC-1)76.77,77 

77  R5-C2*H 
R5»2.*SQRTF( R5 ) 

EM-MU 

DRq-R5/EM 

RIM) -R4 
DO  1  1*1. MU 

J*k'-  ’ 

1  R( J)-R( J+l )-DRB 
I F ( ML-1 ) 76. 76,78 

78  EM-ML 
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TABLE  V  (Continued) 


DRB-R5/EM 
R (  1) *R3 
DO  931  I  *  1 » ML 

931  R ( I +1 ) *R ( I ) +DRB 
R6  =  R5 

76  MM=M-ML-MU 
M 1 *MM~ 1 
EMl*Ml 

DRB-R4-R3-R5-R6 
DRB-DRB/EM1 
MM1-ML+1 
RIMM1 )*R3+R6 
J*MM  1 

DO  932  1*1. Ml 
J  =  J  +  1 

932  R<  J)*R(  J-D+DRB 
IF(NORM)680. 680.682 

680  DO  681  I  *  1 »M 
RGU  )*R<  I  )/F2 

681  Y (  I  )  *R (  I ) /F2 
GO  TO  684 

682  DO  683  1  *1. M 
RG ( I  )  *R ( I  >/F2 

683  Y( I )*R< I ) 

684  G*E*H**3/ ( 12.*( l.-PRT**2  )  ) 
DRB-DRB/F2 

PM*PRT-1. 

PP-PRT+1 . 

P 1  *  1 . -PRT 

T*H/SQRTF ( 12. *( l.-PRT**2 ) ) 

I F ( N ) 84 . 84 . 46 

46  THETA3*THETA3*PI 
THETA4*THETA4*PI 
Nl-N-1 

ANl'Nl 

THETA5*THETA4-THETA3 
DELTT-THETA5/AN1 
THETA! 1 ) *THE  T  A3 
DO  47  J« 1 , N 1 
J1*J+1 

47  THETAIJ1 )«THETA( JI+DELTT 

C 

C  -  -  (3)  GENERATE  THE  THOMSON  FUNCTIONS 

C 

84  DO  5  1*1. M 
X=C*R ( I ) 

Xl*l./X 

I F ( X-6 • 0  )  2 . 2 . 3 

2  LONG* 1 0 
GO  TO  4 

3  LONG* 16 

4  BBER (  I  ) *BER ( X • LONG ) 

BBEI ( I ) *  B  E I ( X .LONG ) 

BBRP ( I ) *BRP ( X . LONG ) 

BBIP! I ) *B I P ( X • LONG ) 

AAK.ER  (  I  )  *  AKER  (  X  .LONG  ) 

AAKEI ( I ) * AKE I (X.LONG) 

A AKERP ( I >=AKERP(X.LONG) 
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TABLE  V  (Continued) 


AAKEIPI I )  ”AKE I P I  X  *  LONG ) 

SBER ( I ) »-BBE I ( I >-Xl*BBRP< I ) 

SBE I ( I ) “BBER ( I ) -Xl *BB I P (  I  ) 

SICER  (  I  >»-AAlCEI  I  I  ) -X1*AAXERP  (  I  ) 

SICE  I  (I  )»AAICER<  I  > -X  1*AAICE  I P  I  I) 

TBER ( I ) »BBE I ( I )/X-BBlP( I >+2.*BBRP(  I  I  /X**2 
TBE I ( I ) “BBRP I I ) -BBER ( I )/X+2.*BBIP( I ) /X**2 
TXER  (  I  )”AAICE  I  (  I  l/X-AAICEIP  I I ) +2.*AAXERP ( I ) /X**2 
5  TKEI  (  II-AAXERP  (  I  ) -AAXER  (  I  ) /X  +  2.  *AAICE  I P  (  I  I /X**2 
C 

C - (4)  CALCULATE  THE  CONSTANTS  OF  INTEGRATION  FOR  THE  AX  I -STMHETR I C 

C  DEFORMATIONS 

C 

I F I  NORM ) 32  *  32 .31 

31  CNO-1. 

CNO 1 » 1  • 

CN02-1. 

CN03-1. 

CN08»1. 

GO  TO  202 

32  CNO»F2*RHO*H 
CN01”E/(F2**2*RH0) 

CN02”F2*RH0*H**2 

CN03”RH0*H*SQRTF(F2*H> 

CN08”RH0*F2**2/(E*SQRTF(F2*H) ) 

202  Z”R3 
X”C*Z 
Xl-l./X 
X2”C/Z 
L-l 
MM”  1 

IF(ICl-2)29»28»85 
85  I  F ( *1-3)28.28.34 
26  Z-R4 
X”C*Z 
Xl”l./X 
X2-C/Z 
MM”M 

IFIX2-2129.28.28 
C  W(R)”0. 

28  BIL. 1 )”BBER (MM) 

B ( L  *  2 ) ”BBE I (MM) 

BIL»3)”AAICER(MM) 

B  I  L  <  4  ) ”AAXE I ( MM ) 

B ( L»  5 ) ”1  • 

B ( L  *  6 ) *0 • 

BB(L)-RH0*P1*Z**2*.5/E 
L-L  +  l 

C  U(R)”0. 

BIL. 1 )”X1*PP*BBIP(MM) 

B(L.2)”-X1*PP*BBRP(MM> 

BIL.  3  )»X1*PP*AAXEIP(MM) 

B(L.4)”-X1*PP*AAXERPIMM) 

B I L . 5  )  ”1  • 

BIL. 6)—  PP*X1**2 
BB(L)-2.*F**2*RHO*PM/E 
L”L  +  1 

IFIL-3I203.203.56 
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TABLE  V  (Continued) 


203  IF(ICl-3)30.29,34 
56  I F ( K2-3 130,29,34 
C  WP ( R ) «0. 

30  B(L,1 )«C*BBRP(MM> 

B(L,2)«C*BBIP(MM> 

B(L,3)»C*AAICERP(MM) 

B(L.4)»C*AAKEIP(MM) 

B ( L , 5 ) *0 • 

B<L,6)«0« 

BB(L)«RH0*P1*Z/E 
L*L  +  1 

IFIL-4 126,26, 36 
C  MR(R)»0. 

29  B(L,1)»-(C**2> *BBE I ( MM ) +PM*X2*BBRP ( MM ) 
B(L.2)=C**2*BBER(MM>+PM*X2*BBIP(MM) 
3(L.3)=-(C**2)*AAKEI  ( MM ) +PM*X2* AAK ERP ( MM > 
B(L.4)«C**2*AAKER(MM)+PM*X2*AAKEIP(MM) 

B ( L. 5  I «0. 

B(L.6)«0. 

BB ( L ) »RH0*H**3/ ( 1 2  .*G  ) 

L«=L  +  1 

IF(L-4)35,26.73 
73  IF(L-5)35,35.36 
C  NR(R)-»0. 

35  B(L,1 ) »-X2 *BB I P ( MM ) 

B(L.2)*X2*BBRP(MM) 

B(L,3)=-X2*AAKEIP(MM) 

B(L,4)*X2*AAXERP(MM) 

B  (  L  ,  5  )  =0« 

B ( L 1 6 ) *  1 • /Z**  2 
BB ( L ) «-F*RH0/ ( E*T  ) 

L  =L  +  1 

C  OR ( R ) =0 • 

B(L.  1  >»-C**3*BBIP(MM> 
b  (  L  *  2  )  »C**3*BBRP(MM> 
B(l.3)»-C**3*AAKEIP(MM) 
B(L.4)=C**3*AAKERP(MM) 

B  (  L , 5 ) =0« 

B(L.6)=0. 

BBIL )-0. 

L  =L  +  1 

IF(L-4>26.26.36 

36  DET»1. 

MATR  =  XSIME0F (6*6* 1*B*BB*DET  *B9) 

GO  TO  (6.7,8) .MATR 
6  B1»B< 1.1) 

B2*B ( 2 , 1 ) 

B3*B( 3 .1 ) 

B4»B ( 4,1 ) 

B5=B( 5.1 ) 

B6«B (6*1 ) 

GO  TO  42 
34  X2«1./(C*R2) 

X3*C/R2 

C  W ( R  2 ) *0  « 

0 ( 1.1 ) =BBER ( M  ) 

0(1.2) »BBE I(M) 

D  (  1 »  3  )  =  1  • 
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TABLE  V  (Continued) 


DD( 1 )-RHO*Pl*R2**2*.5/E 
C  U(R2)«0. 

D(2,1)*X2*PP*BB1P(M) 

D(2,2)«-X2*PP*BBRP(M) 

D(  2  *  3 ) *  1  • 

DD(2  )«2.*F**2*RH0*PM/E 
1 F( K2-2 ) 37. 37  *38 

37  0(3,1 )«C*BBRP(M) 

0(3,2) *C*BB 1 P ( M ) 

D ( 3 »  3 ) *0 • 

DD(3)»RH0*P1*R2/E 
GO  TO  43 

C  MR ( R2 ) *0 • 

38  D( 3. 1 )»-C**2*BBEl ( M ) +PM*X3*BBRP ( M ) 
D(3.2)«C**2*BBER(M)+PM*X3*BB1P(M> 

0(3*  3)*0» 

00(3) «RH0*H**3/ ( 12.*G) 

43  DET* 1 • 

MATR*XS1MEQF( 3,3.1 ,D. DO. DET. D3) 

GO  TO  (44,7.8)  .MATR 

44  B 1 *0 ( 1 » 1 ) 

B2*D (2.1) 

B  3*0  • 

B4*0» 

65*0(3.1 ) 

B6*0  • 

C 

C  -  -  (5)  GENERATE  THE  DEPENDENT  VARIABLES  FOR  THE  AX  I -S YMMET R I C  OEFORMA 
C  T IONS 

C 

42  DO  9  1*1. M 
X«C*R ( 1  ) 

Xl*l./X 

P(  I  ) *B1*BBER(  I  ) +B2*BBE I ( I )+B3*AAKER( I )  +  B4*AAICEI  ( I )+B5-RH0*Pl 
1  *R  (  I  )**2*.5/E 
PWB ( I ) * ( P ( I ) -B5 ) *CN01 
P ( I ) *P ( I )*CN01 

Will )«B1*C*BBRP( I ) +B2*C*BB I P ( I ) +B 3*C*AAKERP ( I ) +B4*C* A AXE  I P  < I l-RHO* 
IP 1*R ( I ) / E 
W1  (  I  ) *W1 ( I  ) /CN08 

CMR1«B1*(-C**2*BBEI ( I )+PM*C/R( I )*BBRP( I ) ) +B2* ( C**2*BBER ( I )+PM*C/R( 
1  I ) *BB I P (  I  )  ) 

EMR ( 1 )«-G*(CMRl+B3*(-C**2*AAKEI ( I )+PM*C/R< I )*AAXERP( I ) )+B4*(C**2*A 
1 AKER ( 1 ) +PM*C/R ( I )*AAKE I P ( I ) ) ) +RH0*H**3/ 1 2 . 

EMR ( 1 ) *EMR ( I I/CN02 

ENR( I )«E*T*H*(B2*C/R( I )*BBRP( I ) -B 1 *C/R ( 1 )*BBIP< I ) -B3*C/R ( I ) *AAKE I P 
1 ( 1 )+B4*C/R ( I ) *AAKERP ( I )+B6/R( I ) **2 ) +F*RHO*H 
E  NR ( I ) -ENR ( I l/CNO 

OKI )*-G*(-Bl*C**3*BBIP( I ) +B2*C** 3*BBRP ( I ) -B3*C**3*AAXE I P ( I )+B4*C* 
1 *3* AAXERP ( 1 ) ) 

01(1 >>01(1 I/CN03 

CMT1«B1*(-PRT  *C**2*BBEI ( I )+Pl*C/R ( I )*BBRP( I ) )+B2* (PRT  *C**2*B 
1BER ( 1 )+Pl*C/R( I ) *BB 1 P ( 1 ) ) 

EMT ( I )«-G*(CMTl+B3*(-PRT  *C**2*AAKEI (I)+P1*C/R(I )* AAKERP ( I ) )+B4* 
1 ( PRT  *C**2*AAXER( 1 )+Pl*C/R( I )*AAKE1P( I ) ) )+RHO*H**3/12. 

EMT ( I ) *EMT ( 1 )/CN02 

ENT1-B1* (-BBERI I ) +X1*BB I P ( I ) ) -B2* ( BBE I < I )+Xl*BBRP( I ) ) +B3* ( -A AKER ( I 
1 ) +X 1 *AAKE I P (  1  )  ) 
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TABLE  V  (Continued) 


ENT (  I  )=E*H*F1*(ENT1-B4*(AAKEI ( I )+Xl*AAKERP< I ) ) -B6/X**2 ) +F*RHO*H 
ENT ( I  )*ENT( I >/CNO 

U 1 »8  1  *X1 *8B I P ( I  >-B2*Xl*B8RP( I ) +B3*X 1 *AAKE I P ( I ) -84*X 1 *AAKERP (  I )+85/ 
1PP-86*X1**2 

0(  I  )«R< I )*F1*PP*U1+F*RH0*P1*R( I )/E 
0 ( I ) *0 ( I ) *CN0 1 
9  CONTINUE 

IF(K3)402.402.45 

C 

C  -  -  (6)  CALCULATE  THE  CONSTANTS  OF  INTEGRATION  FOR  THE  ANT  I -SYMMETR I C 

C  deformations 

C 

45  Z*R3 
X*C*Z 
Xl*l./X 
L=  1 
MM  *  1 

I F ( <1-2 ) 54, 53.87 
48  Z=R4 
X«C*Z 
XlM./X 
MM  =  M 

I F  t K2-2  >54,53.53 
87  IF(Kl-3)53,53,60 
C  W  t  R 1 ) *0 • 

53  A  i  L  »  1  )  *BBRP ( MM ) 

A(L.2)*BBIP(MM) 

A(L,3)=AAKERP(MM) 

AIL,4)=AAKEIP(MM) 

A  (  L  ,  5  >  *  X  1 
A  (  L  ,  6  )  «  X 
A  (  L  ,  7  )  =  0  • 

A (L, 8 1=0. 

S ( L ) =0 • 

L  =  L  + 1 

C  U ( R ) =0 • 

AIL, 1 ) =PP*S8E I (MM ) 

A  I  L , 2  )  =-PP*SBER I  MM ) 

A  ( L , 3  >  =PP*  SKEI(MM) 

A  (  L , 4  )  =-PP*  SKER(MM) 

A(L,5)«2.*LOGF(X) 

A ( L , 6 ) =X**2*  •  5 
\ I L , 7 ) =PP*X 1**2 
ML, 81*1, 

S (L  >  =  .2  5*X*Z*F*RHO*( l.+5.*PRT  > /E 
L  *L  + 1 

C  V(R)«0. 

A  (L, 1  )  «PP*X1*BBIP(MM) 

A  I  L , 2 >  =-PP*Xl*BBRP (MM) 

A (L, 3) =PP*X1*AAKEIP(MM ) 

A ( L , 4  )  =-PP*Xl*AAKERP(MM) 

A( L, 5  )  *2 • *LOGF ( X ) 

A ( L , 6 ) =-X**2* • 5 
A(L,7)=-PP*X1**2 
A  (  L ,  8  )  *  1  ♦ 

S ( L  )  «.2  5*X*Z*F*RHO* (ll.  +  7.*PRT)/E 
L*L  +  1 

IF(L-4)57,57,58 
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TABLE  V  (Continued) 


57  I  F I K 1 -3 ) 55 1  54 , 60 

58  IF( <2-3)55.54,60 
C  WPIRI-0. 

55  AIL.1 )«SBER ( MM ) 

A  I  L  » 2 ) *SBE I  I  MM ) 

A(L,3»«SKER(MM) 

A  (L.4)*S<EI I MM) 

A ( L . 5 ) *-X 1**2 

AIL. 61*1. 

A  I L . 7 ) *0 • 

A  I L »  8 ) *0. 

S I L ) *0 • 

L*L+1 

I F I L-5 ) 48  .48  »  6 1 
C  MRIRI-O. 

54  A (L, 1 )«TBER (MM >+PRT  *X 1  *  I SBER i MM ) -x 1 *BBRP ( MM ) ) 

A (L. 2 )«TBEI I  MM) +PRT  *X 1* ( SBE I ( MM ) -X 1 *BB I P ( MM )  ) 
A(L«3)*TKER(MM)+PRT  *X1* I SKER ( MM ) -X1*AAKERP ( MM ) ) 

A(L.4)»TKEI (MM>+PRT  *X 1* I SKE I  I  MM ) -X 1 *A AK E I P I  MM )  ) 

A  I L. 5 ) *P 1*X 1**3 
A  I L . 6 ) *0  » 

A  I L , 7 ) *0. 

A ( L , 8 ) =0  » 

S I L ) *0. 

L*L  +  1 

I F I L-5 ) 62 . 48 . 63 
63  1 F ( L-6  >62.62.61 
C  NR(R)»0. 

62  A (L. 1 ) -BBER I  MM ) -2  .*X 1 *BB I P ( MM ) 

A(L.2)»BBEI (MM)+2.*X1*BBRP(MM) 
A(L.3)«AAKER(MM)-2.*X1*AAKEIP(MM) 

A  I L . 4 ) «AA<E 1 (MM)+2.*X1*AAKERP(MM) 

A(L.5)*-1./PP 
A (L.6  )*0. 

A  I L »  7 ) «2  «  *X 1**2 
A ( L » 8 ) *0« 

SIL) «-X*Z*RHO*F*.5/E 
L*L  +  1 

C  QR(R)+1./ALPHA2*DMRT/DTHETA=0. 

AIL. 1 ) «G*C*SBE 1 (MM)-(G*P1*X1/Z)*(BBEI I  MM ) +2 . *X 1 *B0RP I M )  ) 
A(L.2)«-G*C*SBER(MM)  +  (G*P1*X1/Z)*(BBER(MM)-2.*X1*BBIP(M)  ) 

A (L.3 ) «G*C*SKEI I  MM ) - I G*P 1 *X 1 /Z ) *  I AAKE I  I  MM ) +2 . *X 1*AAKERP I  MM  )  ) 
A!L,4)«-G*C*SKEI(MM)  +  (G*P1*X1/Z)*(AAKER(MM)-2«*X1*AAKEIP(MM)  ) 
A(L.5)«-PM*G*C**2*F1*X1/(S0RTF( 1.+<Z*F1 )**2) ) **3*2 . *X 1 **2 
A  I L . 6 ) *0 » 

A  I L.7 )’0. 

A  I  Ltd ) =0* 

S(L)*0. 

L  *L  +  1 

C  NRT I R ) *0» 

AIL.1 )»BBERIMM)-2.*X1*BBIP(MM) 

A(L.2)*BBE1 (MMI+2  .*X1*BBRP I  MM ) 
A(L.3)*AAKER(MM)-2.*X1*AAKEIP(MM) 

A  I L . 4 ) *  AAKE 1  I MM)+2.*X1*AAKERP(MM) 

A(L.5)«1./PP 
A  I L . 6 ) *0 • 

A(L.7)«2.*X1**2 
A  I L, 8 ) *0. 
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TABLE  V  (Continued) 


S(L)«1.5*RH0*F*X*Z/E 
L  SL  + 1 

IF(L-5)48, 48.61 
61  DET* 1 • 

MATR*XS I MEQF ( 8*8*1 .A.S.DET.Z3) 

GO  TO  ( 206.7*8 ) .MATR 
206  A  1  =  A ( 1,1) 

A2* A ( 2 • 1 ) 

A3-A(3.1 ) 

A  4* A ( 4.1  ) 

A5  *  A  ( 5*1) 

A6* A ( 6 • 1 ) 

A  7  *  A  (  7  *1  ) 

A  8  *  A  (  8 . 1  ) 

GO  TO  68 

60  I F ( K2-2 )64. 64*65 
C  WP ( R2 ) * 0  • 

64  AAt  1 .1  ) =  SBER ( M ) 

A  A  t  1,2)»SBEI(M) 

AAt  1*3)«1, 

AAt  1  *  4 ) e  0 • 

5St 1 ) *0 • 

GO  TO  66 
C  MR ( R2 ) *  0 • 

6  5  A  A  t  1  .1  )=TBER( M)+PRT  /  t  C*R2)*(SBER(M)-BBRP(M)  /  (  C*R2  )  ) 
A A  t 1  .2 ) «TBEI ( M ) +PRT  / ( C*R2 )*(SBEl(M)-BBIP(M)/( C*R  2 )  ) 
AA ( 1 . 3) *0. 

A A ( 1*4) *0 • 

sst 1 ) *0 • 

C  W(R2)«0. 

66  AAt  2*1  )  - BBR  P ( M ) 

AA ( 2 , 2 ) s BB I P ( M ) 

A  A ( 2 . 3 ) 3C  *R  2 

AA ( 2 .4  )  =0. 

SS ( 2 ) -0 • 

:  V ( R2  )  =0 • 

AAt  3.1  )«PP*BBIP(M) / t  C*  R2 ) 

A A ( 3 .2 ) *-PP*BBRP( M) / ( C*R2 ) 

AA ( 3 .3 ) =-.5*( C*R2 )**2 
A  A  (  3 , 4  )  *  1  • 

SS ( 3 ) *.2  5*C*R2**2*F*RHO* ( 11 .+7.*PRT  )  /  E 
C  U(R2)»0. 

A  A ( 4  *  1 ) =PP  *  SBE I  (M) 

AA (4.2 ) =-PP#SBER(M) 

A A ( 4 • 3 ) =«5#(C#R2 )##2 
AA ( 4 .4  )  «1  . 

SS (4 ) ».2  5*C*R2**2*F*RHO*( 1 •  +  5.*PRT  )  /E 
OET  *  1  • 

MATR  =  XS I MEQF ( 4 . 4  *  1  * AA . SS * DE T » ZZ ) 

GO  TO  (67.7*8)  *  MAT  R 

67  A 1 * A A ( 1.1) 

A2  * AA (2.1) 

A3-0. 

A4*0  • 

A5  *0  • 

A6*AA (3.1 ) 

A 7  *0  • 

A8*AA (4,1 ) 
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GO  TO  68 

7  WRITE  OUTPUT  TAPE  3.120 
GO  TO  99 

8  WRITE  OUTPUT  TAPE  3.121 
GO  TO  99 

C 

C  -  -  (7)  GENERATE  THE  DEPENDENT  VARIABLES  FOR  THE  ANT  I -SYMMETR  I  C  DEFORM 
C  ATIONS 

C 

68  DO  71  1*1 .M 
X«R< I )*C 
X1M./X 

WW<  I  )  »A1*BBRP < I  )  + A2*BB I P  < I  )  +A3*AAKERP ( I  )+A4*AAKEIP< I  )+A5*Xl+A6*X 
W A5 ( I  )  * ( WW ( I  ) - A6*X  ) *CN01 
WWI I ) *WW ( I ) *CN01 

W3(I)*C*(A1*SBER(I ) +A2*SBE I ( I )+A3*SKER< I  )+A4*SKEI <  I ) -A5*Xl **2+A6 ) 
W3 ( I  )  *W3( I > /CN08 

DMR1*A1* ( TBER ( I l+PRT  *X1*(SBER< I  )-Xl*BBRP< I  111 +A2* < T BE  I < I l+PRT 
1  *Xl*(SBEI ( I >-Xl*BBIP< I ) ) )+A3*(TKER( I l+PRT  *X 1* ( SKER ( I ) -X 1*AAK 
2ERPI I ) ) ) 

DMR (  I  )=-G*C**2*<DMRl+A4*<TKEI (  I  l+PRT  *X1*( SKEI ( I  )-Xl*AAKElP<  I )  ) 

1  1+A5*  <Pl*Xl**3 1 ) 

DMR ( I ) *DMR ( I ) /CN02 

DMT  1*A 1* ( PR  T  *TBER( I  )+Xl*(SBER< I )-Xl*BBRP(  I  1  1  )+A2*(PRT  *TBE I 

1 ( I  1  +  X 1  * ( SBE I (  I  )  -  XI *BB I P (  I  ) )  1  +  A3* ( PRT  *TKER < I  1 +X 1* < SKER < I  1  -XI* 

2AAKERP ( I ) ) ) 

DMT  (  I  )  «-G*C**2*(DMTl  +  A4*  (PRT  *TKE  I  (  I  )+Xl*(SKEI  (  I  ) -X  1  *AAICE  I  P  <  I  )  )  ) 

1-A5* ( P 1*X 1 **3 ) ) 

DMT ( I ) *DMT ( I ) /CN02 

DMRT1*-A1*(BBEI ( I ) +2. *X 1 *BBRP ( I ) )+A2*(BBER( I ) -2 . *X 1 *BB  I  P (  I  ) >-A3*( 
1AAXEI (  I  )+2.*Xl*AAXERP(  I  )  ) 

DMRT ( I  )*PM*G*C**2*X1*(DMRT1+A4*(AAkER( I  )-2.*Xl*AAKElP( I  )  ) -A5* ( 2 . * 
IX 1 **2  )  ) 

DMRT ( I ) *DMR  T ( I 1/CN02 

OOR(  I ) =-G*C**3*(-Al*SBEI  ( I )+A2*SBER( I )-A3*SKE I ( I )+A4*SKER ( I  )  ) 

OUR ( I ) =QQR ( I ) /CN03 

OOT  (  I  )=-G*C**3*Xl*(-Al*BbIP( I )+A2*BBRP( I  )-A3*AAKElP( I )+A4*AAKERP 
1(1)) 

OOT ( I ) “GOT ( I ) /CN03 

DNR1*A1*(BBER(  I  ) -2  .*X 1 *BB I P ( I))+A2*(BBEI(I)+2.*Xl*BBRP(I))+A3*(AA 
1KER ( I )-2.*Xl*AAKEIP( I ) )+A4*(AAKEI ( I ) +2 . *X 1 *AAkERP ( I ) ) 

DNR ( I )=-E*H*Fl*Xl*(DNRl-A5/PP+A7*( 2.*X1**2 ) >-RHO*H*R (I ) / 4 • 

DNR ( I >=DNR< I l/CNO 

DNT  1*-A1*TBE I ( I  ) +A2*T BER ( I ) - A3*TXE I ( I  )  +  A4*TKER< I ) -A5 /PP*X 1  +  A7* ( 2. 
1*X 1**3 ) 

DNT ( I )«E*H*Fl*DNTl+5 • *RHO*H*R ( I )/4. 

DNT ( I ) =DNT ( I )/CNO 

DNRT ( I  ) *E*H*F 1*X 1  * ( DNR 1 +A5/PP  +  A7* ( 2 . *X 1 **2 ) )-3.*RHO*H*R ( I )/4. 

DNRT ( I  )  =DNR T (  I  )/CNO 

U1=-PP*(-A1*SBEI ( I )+A2*SBFR( I )-A3*SKEI ( I )+A4*SKER< I ) ) +A 5* ( 2 . *LOGF 
1 ( X ) ) +A6*X**2*.5+A7*(PP*Xl**2)+A8 
UU(  I  ) -Fl/C*Ul-RHO*( l.  +  5.*PRT I *R ( I ) **2/ ( 8 . *E ) 

UU< I  )*UU( I  ) *CN01 

V1*-PP*X1* ( -A  1 *BB I P ( I ) +A2*BBRP ( I )-A3*AAKEIP( I )+A4*AAKERP( I  ) ) +A5*< 
12. *LOGF ( X ) ) -A6* ( . 5 *X**2 ) -A7* ( PP*X 1**2 ) +A8 
VV(I )«F1/C*V1-RH0*(11.+7.*PRTI*R( I )**2/(8.*E) 

71  VV (  I  )  =W<  I  ) *CN01 
C 
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C  -  - 
C 

402 

C 

C  -  - 

c 


183 

184 

185 

13 

14 

15 

17 

22 

16 

18 

19 

21 

33 


3000 


3001 


3002 


(8)  INPUT  II  ( CONTROL  »  LOAD ) 

READ  INPUT  TAPE  2  ,404 , PS  I ♦MOOT  0 


(9)  OUTPUT  I  (C ON TROLtGEOMETRY, MATERIAL, EDGE  CONSTRAINTS) 


K  3  *  K  3- 1 

I F ( NORM ) 184,184.183 
WRITE  OUTPUT  TAPE  3,118 


GO  TO  185 
WRITE  OUTPUT  TAPE 
IF( Kl  — 2  )  13. 14.15 
WRITE  OUTPUT  TAPE 
GO  TO  16 

WRITE  OUTPUT  TAPE 
GO  TO  16 
I F ( K 1-4 ) 17.22,22 
WRITE  OUTPUT  TAPE 
GO  TO  16 

WRITE  OUTPUT  TAPE 
IF( K2-2 ) 18.19,21 
WRITE  OUTPUT  TAPE 
GO  TO  33 

WRITE  OUTPUT  TAPE 
GO  TO  33 

WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
TPS  I =2 • *PS I  - 1 • 
PSIDEG=PSI#RAD*PI 
PSI3PSI*PI 
T3DEG»THETA3*RAD 
T4DEG*THETA4*RAD 


3,119 

3*129 

3*130 


3.131 
3,136 
3*134 

3.132 

3.133 
3,128 

3 , 1 22 , RHO , E . H , F  ,PRT 
3,11 


AR  3  3 • / ( SQR  T  F ( 1 ,+(R4/F2 ) **2 ) ** 3-SQRT F ( 1 . +  ( R 3 /F 2 ) **2 )**3 ) 


FM1-0. 
FM2  3  0  * 
FM43  0 • 


RM 1  *  0 ♦ 

RM2  =  0 • 

RM43  0 * 

DO  3000  1*1, M 

RM 3 ( P ( I  )*COSF(PSI  )  )**2+(WW(  I )*SINF(P$I )  )**2/2. 
RMl=RMl+RM*4.*RG( I ) /SQR TF ( 1 • +RG ( I ) **2 ) 

FM * ( P WB (  I ) *COSF ( PS  I  )  ) **2+ ( WA5 (  I  )  *S I NF ( PS  I  )  )**2/2. 
FM13FM1+FM*4.*RG( I ) /SORTF ( 1 .+RG( I )**2 ) 

DO  3001  1=2, Ml 

RM= ( P (  I  )*COSF(PSI  )  )**2+( WW(  I  )*SINF (PSI)  ) **2 /2 • 
RM2=RM2+RM*4.*RG( I ) /SQR TF ( 1 . +RG ( I )**2) 

FM3(PWB(  I  )*COSF(PSI )  ) **2+ ( WA5 (  I  ) *S I NF ( PS  I  )  )  **2/ 2 . 
FM23FM2+FM#4*#RG( I ) /SQRTF ( 1 • +RG( I )**2) 

DO  3002  1  =2  , Ml, 2 

RM  =  (P(  I  ) *C  OSF ( PS  I  )  )**2+(WW( I  )#SINF (PSI )  )**2/2, 
rM4=RM4+RM*4.*RG( I )/SQRTF( l.+RG( I )**2) 

FM= ( PWB ( I ) *COSF (PSI ) )##2+(WA5( I )#SINF(PSI ) )*#2/2« 
FM4  =  FM4  +  FM*4.*RG(  I ) /SQRTF  C 1 • +RG C  I  )  **2  ) 

FMl 3AR*DRB/3«#( FM1+FM2+2«*FM4 ) 

RMS=AR*DRB/ 3 • * ( RM1+RM2+2«*RM4 ) 
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487 

187 


186 

188 


189 


190 

C 

C - 

C 


3005 

451 

452 


255 


450 

23 

403 


FMI*SQRTF (FMI ) 

RMS*SQRTF ( RMS ) 

WRITE  OUTPUT  TAPE  3.127 

WRITE  OUTPUT  TAPE  3 . 1 2 3 . R 1 . R2 . R3 . R4 . PS  I  DEG 

WRITE  OUTPUT  TAPE  3.11 

WRITE  OUTPUT  TAPE  3.3003 

WRITE  OUTPUT  TAPE  3 . 3004  .  RMS . FM I 

WRITE  OUTPUT  TAPE  3.11 

I F ( MOOT  0-2 1186.487,48  7 

I F ( MOOT  0-3 ) 187,187 ,186 

WRITE  OUTPUT  TAPE  3,111 

WRITE  OUTPUT  TAPE  3 . 1 1 2 , T3DEG . T40EG 

WRITE  OUTPUT  TAPE  3.11 

IFITPSI 1188.189.188 

WRITE  OUTPUT  TAPE  3.173 


WRITE  OUTPUT  TAPE  3.174.B1 .B2.B3.B4.B5.B6 
WRITE  OUTPUT  TAPE  3.11 
IFIPSI 1190.190.189 
WRITE  OUTPUT  TAPE  3.175 

WRITE  OUTPUT  TAPE  3 . 1 76 , A 1 , A2 , A3 . A4 , A5 . A6 , A7 , A8 
WRITE  OUTPUT  TAPE  3,11 
WRITE  OUTPUT  TAPE  3,10 


(10)  OUTPUT  II  (DEPENDENT  VARIABLES  AT  THE  TABULAR  POINTS) 

IF( MGOTO-3)450. 451 .3005 
I F ( K3 ) 98 , 98 , 402 
IFITPSI 1450.452.450 
WRITE  OUTPUT  TAPE  3.140 
DO  255  J=1.N 
DO  255  1*1. M 

W(  I. J)*WW( I  )*S INF (THETAI J) ) -P ( I  ) 

WFLEX ( I , J ) *WA5 ( I >*SINF( THETA ( J) )-PWB(I) 

WP( I  ,J  )*W3 ( I )*SINF ( THETA (J) ) -W 1(1) 

U(  I  , J  )  *UU ( I )*SINF (THETA( J) )-Q(  I  ) 

V  (  I  ,  J  )  *VV ( I ) *COSF ( THE  T  A ( J ) ) 

SKIP=1 

GO  TO  94 

IFIPSI 196,96.23 

IFITPSI 1403.95.403 

DO  405  1*1. M 

P2( I  )*P( I ) *COSF ( PS  I ) 

PWB2 (  I  ) *PWB ( I  ) *COSF ( PS  I ) 

W 2  (  I  )  *Wl ( I ) *COSF ( PS  I ) 

FMR ( I ) *EMR ( I ) *COSF ( PS  I ) 

FNR (  I  )  *ENR ( I ) *COSF (PS  I ) 

02  (  I  )  *01 ( I ) *COSF ( PS  I ) 

FMT ( I )*EMT ( I )*COSF (PS  I  ) 

FNT ( I ) *ENT ( I ) *COSF ( PS  I  ) 

Q3( I )*Q( I ) *COSF ( PS  I  ) 

W4( I ) *WW ( I ) *S I NF ( PS  I ) 

WA54 ( I ) *WA5 ( I ) *S I NF ( PS  I  ) 

W 5 ( I ) *W3 ( I ) *S I NF ( PS  I ) 

GMR(  I ) *DMR ( I ) *S I NF ( PS  I  ) 

GMT ( I ) -DMT ( I ) *S I NF ( PS  I  ) 

GMRT(I)*DMRT(I ) *S I NF ( PS  I  ) 

QRR ( I ) *QQR ( I ) *S I NF ( PS  I  ) 

QT  T (  I  )  *QQT ( I ) *S I NF ( PS  I  ) 
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GNR ( I ) =DNR ( I ) *S I NF ( PS  I ) 

GNT (  I  ) =  DNT ( I ) *S 1 NF ( PS  I ) 

GNR  T (  I  ) *DNR  T ( I  )  *S I NF ( PS  I ) 

U2 I  I ) -UU  < I ) * S I NF ( PS  I > 

WT  (  I  )  *W4 ( I ) +P2 ( I ) -P ( 1 ) 

WFLEX  T ( I ) =WA54 ( I ) +PWB2 ( 1 ) -PwB ( I ) 

WPT (  I  )  =  W5 ( I ) +  W2 ( I ) -W 1 (  I ) 

UT<  I  )  =  U2 ( I ) +03 ( I ) -0 ( I ) 

405  V21 1 ) «VV(1)*SINFIPS1 ) 

1 F ( MGOTO- 1 ) 96  *  796  *69 

69  DO  70  J=1.N 
DO  70  I  *  1 *M 

W(I.J)=W4(I)*SINF(THETA(J)  ) +P2 (  I  > 

WFLEX (  I  *  J ) =W A54 ( I  )*SINF(THETA(J)  ) +PWB2(  I  ) 

WP (  I .J ) =W5( I ) *SINF ( THETA ( J  )  ) +WZ  I 
CMR(  I >J > =GMR (  I )*SINF(THETA(J) ) +FMR ( I  ) 

CMT  (  I  ,J)=GMT<  I  )  *S  I  NF  (  THETA  (  J  )  )  +FMT  (  I  )' 

CMRT (  I  * J  )  *GMRT ( I )*COSF ( THETA! J)  ) 

OR! I . J) =QRR ( I )*SINF(THETA(J) ) +02 ( I ) 

OT ( 1  » J  ) =QTT (  I  )*COSF( THETA! J)  ) 

CNR!  I  * J  )  =GNR ( I ) *S INF! THETA! J) >  +FNR ( I ) 

CNT (  I  »  J ) =GNT (  I  )*SINF(THETA(J)  )+FNT ( I  ) 

CNRT ( I . J ) =  GNRT ( I ) *COSF ( THETA! J) ) 

J!  1 »  J  )  =U2 (  I  )*S INF (THETAI J)  )+Q3( I  ) 

70  V! I,J)=V2< I )*COSF(THETA(J) ) 

SKIP=0 

IF (MGOTO-2 ) 9  4 *94*215 

215  DO  217  J* 1 »N 
DO  217  1=1.  M 
U(  I . J ) =U ( I .J)-0( I  ) 

W( I *J ) «W(  I • J  > -P ( I ) 

WFLEX! 1 »J)=wFLEX( 1 ,J)-PWB( I ) 

217  WP ( I  .J ) *WP ( I . J)-Wl (  I  ) 

WRITE  OUTPUT  TAPE  3.140 
SK I P  *  1 
GO  TO  94 

96  IF! PSI >230.230.231 

231  WRITE  OUTPUT  TAPE  3.138 

230  IF(NORM>240. 240.241 

240  WRITE  OUTPUT  TAPE  3.324 
GO  TO  242 

241  WRITE  OUTPUT  TAPE  3.124 

242  I F I  PS  I >520.520.521 

520  WRITE  OUTPUT  TAPE  3 . 1 2 5  .  ( Y ( I > . P (  I >  .  W 1 (  I  )  .ENR (  I  >  ,ENT ( I > . I • 1  .M  I 
GO  TO  522 

521  WRITE  OUTPUT  TAPE  3 . 12 5  .  ( Y ( 1 > . P2 ( I >  .  W2 ( I )  .FNR ( I ) . FNT ( I >  .  I • 1 .M > 

522  WRITE  OUTPUT  TAPE  3.10 
I F I  PS  I  >232.232.233 

233  WRITE  OUTPUT  TAPE  3.138 

232  1  F ( NORM  >243.243.244 

243  WRITE  OUTPUT  TAPE  3.326 
GO  TO  245 

244  WRITE  OUTPUT  TAPE  3.126 

245  IFIPSI >523.523.524 

523  WRITE  OUTPUT  TAPE  3 . 1 25  ♦  ( Y ( I > . EMR ( I  >  .  EMT ( I ) .0 1 ( I > .0 ( I  > . I • 1 .M ) 
GO  TO  525 

524  WRITE  OUTPUT  TAPE  3  ♦  1 25  ♦  ( Y (  I  > . FMR ( I )  .FMT ( I > .02 ( I > .03 ( I > ♦ I » 1 . M ) 

525  WRITE  OUTPUT  TAPE  3.10 
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IF( PSI) 900*900*901 

901  WRITE  OUTPUT  TAPE  3*138 

900  IF(N0RM)902*902*903 

902  WRITE  OUTPUT  TAPE  3.380 
GO  TO  904 

903  WRITE  OUTPUT  TAPE  3.480 

904  IFtPSI >905,905,906 

905  WRITE  OUTPUT  TAPE  3 .9 25 . ( Y ( I ) . PWB ( I ) • I » 1 *M ) 

GO  TO  907 

906  WRITE  OUTPUT  TAPE  3 • 925 . ( Y II) . PW82 ( I ) , I - 1 ,M ) 

907  WRITE  OUTPUT  TAPE  3*10 
IFtPSI 199.99.95 

95  IFtJPSI 1234,235*235 

234  WRITE  OUTPUT  TAPE  3,138 

235  IF(N0RM)246,246,247 

246  WRITE  OUTPUT  TAPE  3,313 
GO  TO  248 

247  WRITE  OUTPUT  TAPE  3,113 

248  IF( TPSI  1527,526,526 

526  WRITE  OUTPUT  TAPE  3 , 1 2 5 , < Y ( I  1 , WW ( I  1 , W3 ( I  1 , VV ( I  1 ,UU ( I  1 , 1 - 1 , M ) 

GO  TO  528 

527  WRITE  OUTPUT  TAPE  3 , 1 2 5  * ( Y ( I ) , W4 ( I ) , W5 ( I ) , V2 ( I ) ,U2 (I) , I - 1 ,M ) 

528  WRITE  OUTPUT  TAPE  3*10 
IF<  TPSI  1236,237,237 

236  WRITE  OUTPUT  TAPE  3*138 

237  IF(NORM)249,249.250 

249  WRITE  OUTPUT  TAPE  3*314 
GO  TO  251 

250  WRITE  OUTPUT  TAPE  3,114 

251  IFtTPSI 1530*529*529 

529  WRITE  OUTPUT  TAPE  3  *  1 2 5  * ( Y ( I ) * DMR t I ) ,DMT ( I ) * DMR T ( I ) *QQR ( I ) *  I » 1  * M ) 
GO  TO  531 

530  WRITE  OUTPUT  TAPE  3 . 1 2 5 , ( Y t I  I  * GMR t I ) ,GMT ( I ) * GMR T (I) *QRR t I) *  I -1 *M ) 

531  WRITE  OUTPUT  TAPE  3,10 
IFtTPSI 1238.239.239 

238  WRITE  OUTPUT  TAPE  3,138 

239  I F ( NORM  1252*252*253 

252  WRITE  OUTPUT  TAPE  3,315 
GO  TO  254 

253  WRITE  OUTPUT  TAPE  3.115 

254  IFt TPSI 1533*532*532 

532  WRITE  OUTPUT  TAPE  3  *  1 2 5 . t Y t I ) . DNR t I ) , DNT t I ) *DNR T t I ) ,QQT t I ) . I »1 . M ) 
GO  TO  534 

533  WRITE  OUTPUT  TAPE  3 *  1 25 . t Y t I ) . GNR t I ) ,GNT t I ) , GNR T t I ) ,QTT t I ) *  I - 1 *M ) 

534  WRITE  OUTPUT  TAPE  3,10 
IFtTPSI 1908,909,909 

908  WRITE  OUTPUT  TAPE  3,138 

909  IF(NORM)910*910*911 

910  WRITE  OUTPUT  TAPE  3,380 
GO  TO  912 

911  WRITE  OUTPUT  TAPE  3,480 

912  IFt TPSI1914.913.913 

913  WRITE  OUTPUT  TAPE  3 . 92 5 , t Y t I  1 . WA5 t  I  1 , 1  - 1  ,M  1 
GO  TO  915 

914  WRITE  OUTPUT  TAPE  3 . 925 . t Y t I  1 . WA 54 ( I ) , I - 1 ,M 1 

915  WRITE  OUTPUT  TAPE  3,10 
IFt  K3) 98 *98, 402 

796  WRITE  OUTPUT  TAPE  3,140 
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IF(NORM)797*797»798 

797 

WRITE  OUTPUT  TAPE 
GO  TO  799 

3.327 

798 

WRITE  OUTPUT  TAPE 

3.328 

799 

WRITE  OUTPUT  TAPE 

3.329. 

(Y( I  I .WT( II . WFLEXT ( I ) . 

WPT  < I ) .UT ( I ) » V2 ( I ) . 

11=1. M) 

WRITE  OUTPUT  TAPE 
GO  TO  99 

3.10 

94 

DO  160  J=1.N 

160 

THETR(J)*THETA(J)*RAD 

I  F ( NORM  1601 .601 .604 

601 

IF  (  SK.IPJ602  *602*603 

602 

WRITE  OUTPUT  TAPE 
GO  TO  607 

3.361 

603 

WRITE  OUTPUT  TAPE 
GO  TO  607 

3.461 

604 

IF(SK.IP)605»605.606 

605 

WRITE  OUTPUT  TAPE 
GO  TO  607 

3.161 

606 

WRITE  OUTPUT  TAPE 

3.261 

607 

WRITE  OUTPUT  TAPE 

3.FMT1 

.(THETR(J) *  J*  1 . N ) 

WRITE  OUTPUT  TAPE 

3.11 

WRITE  OUTPUT  TAPE 

3.FMT2 

»(Y(I).(W(I»J).J«1*N), 

1=1. M) 

WRITE  OUTPUT  TAPE 

3.10 

IF (NORM)608. 608.611 

608 

IF(SK.IP)609,609.610 

609 

WRITE  OUTPUT  TAPE 
GO  TO  6 1 A 

3.362 

610 

WRITE  OUTPUT  TAPE 
GO  TO  61 A 

3.462 

611 

I F ( SKIP  1612*612 .613 

612 

WRITE  OUTPUT  TAPE 
GO  TO  6 1 A 

3.162 

613 

WRITE  OUTPUT  TAPE 

3.262 

614 

WRITE  OUTPUT  TAPE 

3.FMT1 

. ( THETR ( J ) .J«l .N) 

WRITE  OUTPUT  TAPE 

3.11 

WRITE  OUTPUT  TAPE 

3.FMT2 

.  ( Y (  I )  .  (WFLEXI I . J ) . J* 1 

*N) . I =1 *M) 

WRITE  OUTPUT  TAPE 

3.10 

I F ( NORMJ615.615.618 

615 

I  F  < SKIPl616.616.617 

616 

WRITE  OUTPUT  TAPE 
GO  TO  621 

3.363 

617 

WRITE  OUTPUT  TAPE 
GO  TO  621 

3.463 

618 

IFISKIP1619.619.620 

619 

WRITE  OUTPUT  TAPE 
GO  TO  621 

3.163 

620 

WRITE  OUTPUT  TAPE 

3.263 

621 

WRITE  OUTPUT  TAPE 

3.FMT1 

.(THETR(J).J-l.N) 

WRITE  OUTPUT  TAPE 

3.11 

WRITE  OUTPUT  TAPE 

3.FMT2 

. ( Y ( I ) . ( WP ( I . J ) . J«1 .N 1 

. 1=1. MI 

WRITE  OUTPUT  TAPE 

3.10 

IFIN0RM1622.622.625 

622 

IFISK.IP1623.623.62A 

623 

WRITE  OUTPUT  TAPE 
GO  TO  628 

3.364 

624 

WRITE  OUTPUT  TAPE 
GO  TO  628 

3.464 
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625 

IFISK.IPI626.626.627 

626 

WRITE  OUTPUT  TAPE 
GO  TO  628 

3  •  1  64 

627 

WRITE  OUTPUT  TAPE 

3.264 

628 

WRITE  OUTPUT  TAPE 

3.FMT1 

.  (THETR(J) »J«1 .N) 

WRITE  OUTPUT  TAPE 

3.11 

WRITE  OUTPUT  TAPE 

3 .FMT  2 

.  m  i  i ,  <uu  .  j i . j-i  ,n i .  i  • 

1  *  M ) 

WRITE  OUTPUT  TAPE 

3.10 

IF (NORM)629.629.632 

629 

IF( SKIPI630.630.631 

630 

WRITE  OUTPUT  TAPE 
GO  TO  635 

3.365 

631 

WRITE  OUTPUT  TAPE 
GO  TO  635 

3.465 

632 

IFISKIPI633.633.634 

633 

WRITE  OUTPUT  TAPE 
GO  TO  635 

3.165 

6  34 

WRITE  OUTPUT  TAPE 

3.265 

635 

WRITE  OUTPUT  TAPE 

3  .FMT 1 

. ( THETRl J  ) » J« 1 • N ) 

WRITE  OUTPUT  TAPE 

3.11 

WRITE  OUTPUT  TAPE 

3.FMT2 

.(Y(I).(V(I.J) • J« 1 , N ) . I« 

l.M) 

WRITE  OUTPUT  TAPE 

3.10 

I F ( SKI P) 599.599.99 

599 

IF(NORM)267t267.268 

267 

WRITE  OUTPUT  TAPE 
GO  TO  269 

3.382 

268 

WRITE  OUTPUT  TAPE 

3.582 

269 

WRITE  OUTPUT  TAPE 

3 . FMT 1 

. ( THETRl J ) • J» 1 »N) 

WRITE  OUTPUT  TAPE 

3.11 

WRITE  OUTPUT  TAPE 

3 .  FMT  2 

. I Y< I  ) . (CMR ( I ,J)  ,J«1.N)  , 

I-l.M) 

WRITE  OUTPUT  TAPE 

3.10 

IFINORMI270.270.271 

270 

WRITE  OUTPUT  TAPE 
GO  TO  272 

3.366 

271 

WRITE  OUTPUT  TAPE 

3.166 

272 

WRITE  OUTPUT  TAPE 

3 .  FMT 1 

. I THETRl J ) >  J* 1 . N ) 

WRITE  OUTPUT  TAPE 

3.11 

WRITE  OUTPUT  TAPE 

3 .  FMT 2 

. I Y I  I ) . I CMT I  I . J ) . J* 1 . N ) . 

1*1 t M) 

WRITE  OUTPUT  TAPE 

3.10 

IF( NORMJ273.273.274 

273 

WRITE  OUTPUT  TAPE 
GO  TO  275 

3.367 

274 

WRITE  OUTPUT  TAPE 

3.167 

275 

WRITE  OUTPUT  TAPE 

3 . FMT  1 

. I THETRl J ) . J*1 .N ) 

WRITE  OUTPUT  TAPE 

3.11 

WRITE  OUTPUT  TAPE 

3 .  FMT  2 

. I Y I  I ) . (CMRTI I .J) . J* 1 . N ) 

.I-l.M) 

WRITE  OUTPUT  TAPE 

3.10 

IF (NORM)276.276.277 

276 

WRITE  OUTPUT  TAPE 
GO  TO  278 

3.368 

277 

WRITE  OUTPUT  TAPE 

3.168 

278 

WRITE  OUTPUT  TAPE 

3.FMT1 

. (THETRl J ) » J«1 »N) 

WRITE  OUTPUT  TAPE 

3.11 

WRITE  OUTPUT  TAPE 

3.FMT2 

. I Y 1 1) . (OR  I  I . J) » J»1 *N)  ,  I 

-l.M) 

WRITE  OUTPUT  TAPE 

3.10 

IF(NORM>279.279.280 

279 

WRITE  OUTPUT  TAPE 
GO  TO  281 

3.369 
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280 

WRITE  OUTPUT 

TAPE 

3.169 

281 

WRITE  OUTPUT 

TAPE 

3.FMT1 

WRITE  OUTPUT 

TAPE 

3.11 

WRITE  OUTPUT 

TAPE 

3.FMT2 

WRITE  OUTPUT 

TAPE 

3.10 

I F I  NORM ) 282  » 

282.283 

282 

WRITE  OUTPUT 

TAPE 

3.370 

GO  TO  284 

283 

WRITE  OUTPUT 

TAPE 

3.170 

284 

WRITE  OUTPUT 

TAPE 

3.FMT1 

WRITE  OUTPUT 

TAPE 

3.11 

WRITE  OUTPUT 

TAPE 

3.FMT2 

WRITE  OUTPUT 

TAPE 

3.10 

I F I  NORM  1285* 

285.286 

285 

WRITE  OUTPUT 

TAPE 

3.371 

GO  TO  287 

286 

WRITE  OUTPUT 

TAPE 

3.171 

287 

WRITE  OUTPUT 

TAPE 

3.FMT1 

WRITE  OUTPUT 

TAPE 

3.11 

WRITE  OUTPUT 

TAPE 

3.FMT2 

WRITE  OUTPUT 

TAPE 

3.10 

I F I  NORM ) 288  t 

268.2 

89 

288 

WRITE  OUTPUT 

TAPE 

3.372 

GO  TO  290 

289 

WRITE  OUTPUT 

TAPE 

3.172 

290 

WRITE  OUTPUT 

TAPE 

3.FMT1 

WRITE  OUTPUT 

TAPE 

3.11 

WRITE  OUTPUT 

TAPE 

3.FMT2 

WRITE  OUTPUT 

TAPE 

3.10 

500 

I F ( MGOTO-3 ) 99  *  99  * 

215 

99 

IF  ( K3 ) 98  »  98  * 

402 

101 

FORMAT (715) 

102 

FORMAT ( 3E2  0* 

8  ) 

103 

FORMAT (4F15. 

9  ) 

10 

FORMAT ( 1 H 1 ) 

11 

FORMAT(////) 

1 1  1 

FORMAT ( 51H 

112 

FORMAT ( 2F24. 

2  I 

113 

FORMAT ( 116H 

1  OMEGAR 

114 

FORMAT  t 12 iH 

1  MTHET A  t IN 

.-lb. 

/IN.  I 

2////  ) 

115 

FORMAT  1 121H 

1  NTHETAtIN./ 

LB.) 

>1  .N) 


♦  J ) »  J* 1 » N ) *  I  * 1 » M ) 


♦  J  *  1  ♦  N ) 


t (CNR (  I  *  J  ) * J*1 .N ) *  I  *  1 *M ) 


■1  •  N  ) 

!  * J)  * J* 1  * N ) .  I-l.M) 


'1  »N ) 


. ( CNR T ( I  ♦  J) .J*l .N ) .1-1 .M) 


THE  T  A3 ( DEG •  ) 

R(  IN.  ) 

V( IN.  ) 

R (  IN.) 

MRTHETA ( IN. -LB. /IN. ) 
R (  IN.) 

NRTHETAt IN. /LB. ) 


THETA4 ( DEG.  )////  ) 

W ( IN.) 

U( IN. )////) 
MR ( IN. -LB. /IN. ) 

QR  ( IN. /LB. ) 

NR ( IN./LB. ) 

QTHE  TA (  IN. /LB. ) 


if  n  n 

118  FORMAT ( 99H  A 
1LIZED  RESULTS 

119  FORMAT ( 99H  A 
1  ED  RESULTS  BY 

120  FORMAT(30H 
FORMAT ( 30H 


121 

122 

123 

124 


PARABOLOIDAL  SHELL  SUBJECTED 
BY  A  SHALLOW  SHELL  ANALYSIS 
PARABOLOIDAL  SHELL  SUBJECTED 
A  SHALLOW  SHELL  ANALYSIS 
UNDERFLOW  OR  OVERFLOW 
THE  MATRIX  IS  SINGULAR 


TO  GRAVITY-  -  -UN-NORMA 
////) 

TO  GRAVITY - NORMAL  I  2 

////) 

////) 

/// /) 


FORMAT ( 5F24. 5 ) 

FORMAT (4F24.5.F18.2 ) 
FORMAT ( 1 2 1H 
1  OMEGAR 

1////  ) 

125  FORMAT (F24.4.4F24.8) 


R(  IN.  ) 

NR( LB./IN. ) 


W ( IN.) 

NTHETA( LB./IN.  ) 


67 


TABLE  V  (Continued) 

126  FORMAT ( 12 1H 

R  < IN.) 

MR! IN.-LB./IN. ) 

1  MTHETA! IN.-LB./IN. 

1//// ) 

)  ORILB./IN.) 

U ( IN.) 

127  FORMAT ( 1 1 6H 

RUIN.) 

R2 ( IN. ) 

1  R  3  C IN.) 

R4( IN. ) 

PSI (DEG. )////) 

128  FORMAT ( 120H  WEIGHT 

DENS  I TY ( LB./ I N. 3  ) 

YOUNGS 

MODULUS« LB./IN. 

12)  THICKNESS! IN. ) 

FOCAL  LENGTH! IN. ) 

POISSONS  RATIO/ 

2///) 

129  FORMAT ( 51H  THE  SHELL 

IS  FREE  AT  R 1 

////  ) 

130  FORMAT ( 51H  THE  SHELL 

IS  CLAMPED  AT  R1 

////  ) 

131  FORMAT (  51H  THE  SHELL 

IS  SIMPLY  SUPPORTED 

AT  R 1 

////  ) 

132  FORMAT { 5 1H  THE  SHELL 

IS  CLAMPED  AT  R2 

////  ) 

133  FORMAT ( 5 1 H  THE  SHELL 

IS  SIMPLY  SUPPORTED 

AT  R2 

////  ) 

134  FORMAT {  5 1 H  THE  SHELL 

IS  FREE  AT  R2 

////  ) 

136  FORMA T { 5 1 H  THE  SHELL 

IS  CLOSED  AT  THE  APEX 

////  ) 

138  FORMAT ( 51H  THE  PORTI 

ON  OF  THE  RESULTS  INDEPENDENT  OF  THETA////) 

140  FORMAT ( 95H  THE  DISTORTION  OF  THE  SHELL  GIVEN  BELOW 

IS  MEASURED  REL 

1  AT  I VE  TO  THAT  OF  THE  FACE-UP  POSITION  ////) 

173  FORMAT (7X.2HB1.13X.2HB2 

174  FORMAT ( 6E 1 5 . 6 ) 

.13X.2HB3.13X.2HB4. 

13X.2HB5, 

1 3X  *  2HB6 • / / / / ) 

175  FORMAT  ( 7X . 2HA 1 . 1 3X .2 HA2 , 1 3X . 2HA3 . 1 3X . 2H A4 

.13X.2HA5 

*13Xi2HA6*13X» 

12HA7.13Xf 2HA8*//// ) 

176  FORMAT ( 8E15. 6 ) 

205  FORMAT ( 12A6 ) 

313  FORMAT ( 1 16H 

gamma 

W* 

1  OMEGAR* 

V* 

U*  ////) 

314  FORMAT ( 1 16H 

GAMMA 

MR* 

1  MTHETA* 

mrtheta* 

QR*  ////) 

315  FORMAT ( 1 16H 

gamma 

NR* 

1  NTHETA* 

nrtheta* 

QTHFTA*//// ) 

324  FORMAT ( 1 1 6H 

GAMMA 

W* 

1  OMEGAR* 

NR* 

NTHETA*//// ) 

326  FORMAT ( 1 1 6H 

GAMMA 

MR* 

1  MTHETA* 

QR* 

U*  ////) 

327  FORMAT ( 12 1H 

GAMMA 

W  TILDE* 

WFLEX 

1TILDE*  OMEGAR  TILDE*  U  TILDE* 

V  TILDE* 

2/// /  ) 

328  FORMAT ( 1 2 1H 

R(  IN.  )  W 

TILDE! IN. 

)  WFLEX  T I 

1LDEUN.)  OMEGAR  TILDEIIN*)  U  TILDEUN.) 

V  T I LDE ( IN. ) 

2//// ) 

329  FORMAT (F20.4.5F20.8) 

380  FORMAT ( 51H 

404  FORMAT ( F15 .9  *  15) 

GAMMA 

WFLEX*  ////) 

480  FORMAT ( 51H 

RUN.  ) 

WFLEX (IN.)  ////) 

925  FORMAT (F24.4.F24.8) 

3003  FORMAT ( 51H 

EPI SILON 

EPISILONFLEX  ////) 

3004  FORMAT ( 2F24. 8  ) 

161  FORMAT (50H 

- displacement 

W!  IN. ) - 

/////) 

162  FORMAT (50H 

- DISPLACEMENT 

WFLEX! IN 

•)  ////) 

163  FORMAT (50H 

- displacement 

OMEGAR - 

/////) 

164  FORMAT (50H 

- DISPLACEMENT 

U(  IN.  ) - 

/////) 

165  FORMAT ( 50H 

- displacement 

V!  IN.  ) - 

/////) 

582  FORMAT (60H 

- BENDING  MOMENT  MR! IN.- 

LB./ IN. ) - 

1  /////) 

166  FORMAT ( 60H 

- BENDING  MOMENT  MTHE TA ( I N . -LB . / I N . ) - 

1-  /////) 

167  FORMAT ( 60H 

- BENDING  MOMENT  MR THET A ( I N. -LB . / I N. ) — 

1 —  /////) 
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168  FORMAT (50H 

- TRANSVERSE  SHEAR  QRILB./IN.) - /////) 

169  FORMAT ( 60H 

- TRANSVERSE  SHEAR  QTHET A ( LB . / IN . ) - 

1  /////) 

170  FORMAT ( 60H 

- STRESS  RESULTANT  NR ( LB. / I N • > - 

i  nun 

171  FORMAT ( 60H 

- STRESS  RESULTANT  NTHET A ( LB . / I N . ) - 

1  /////) 

172  FORMAT ( 60H 

- STRESS  RESULTANT  NRTHET A ( LB . / I N .  ) - 

i  nun 

261  FORMAT! 50H 

- DISPLACEMENT  W  TILDE ( IN. ) -  /////) 

262  FORMAT! 50H 

- DISPLACEMENT  WFLEX  TILDEIIN.) -  /////) 

263  FORMAT ( 5  OH 

- DISPLACEMENT  OMEGAR  TILDE -  /////) 

26 A  FORMAT (50H 

- DISPLACEMENT  U  TILDEIIN.) -  /////) 

265  FORMAT! 50H 

- DISPLACEMENT  V  TILDEIIN.) -  /////) 

361  FORMAT (75H 

- normalized  displacement  W*(DIMENSIONL 

1ESS) - 

///// ) 

362  FORMAT ( 75H 

- NORMALIZED  DISPLACEMENT  WFLEX*  (DIMEN 

1 S I ONL  ESS ) - 

// / / ) 

363  F  ORM AT i 7  5H 

- NORMALIZED  DISPLACEMENT  OMEGAR* ( D I ME N 

1 S I ONLESS ) 

/////) 

36 A  FORMAT (75H 

- NORMALIZED  DISPLACEMENT  U* ( D I  MENS  I ONL 

1  ESS ) - 

//  /  /  /  ) 

365  FORMAT ( 75H 

- NORMALIZED  DISPLACEMENT  V*(DIMENSI0NL 

1  ESS ) - 

/////) 

382  FORMAT (75H 

- normalized  bending  MOMENT  MR*(DIMENSI 

10NLESS ) - 

/////) 

366  FORMAT ( 7 5H 

- normalized  bending  MOMENT  MTHETA*(DIM 

1ENSI ONLESS ) - 

/////) 

367  FORMAT ( 75H 

- NORMALIZED  BENDING  MOMENT  MRTHETA* ( DI 

1MENS I ONLESS ) - 

/////) 

368  FORMAT (75H 

- NORMALIZED  TRANSVERSE  SHEAR  OR* ( D I MEN 

IS  IONLESS ) - 

/////) 

369  FORMAT ( 75H 

- NORMALIZED  TRANSVERSE  SHEAR  QTHE TA* ( D 

1 I  MENS  I ONLESS ) - 

/////) 

370  FORMAT ( 75H 

- NORMALIZED  STRESS  RESULTANT  NR* ( D I  MEN 

1 S I ONLESS ) - 

/////) 

371  FORMAT ( 75H 

- NORMALIZED  STRESS  RESULTANT  NTHET  A* ( D 

1 I  MENS  I ONLESS ) - 

/////) 

372  FORMAT ( 7  5H 

- NORMALIZED  STRESS  RESULTANT  NRTHETA* ( 

1DIMENSIONLESS) - 

/////) 

A61  FORMAT ( 75H 

- NORMALIZED  DISPLACEMENT  W  T I LDE* ( D I  MENS  I 0 

INLESS ) - 

/////) 

A62  F ORMAT ( 7 5H 

- NORMALIZED  displacement  WFLEX  TILDE*  (DIM 

1 ENS I ONLESS ) - 

/////) 

A63  FORMAT ( 75H 

- NORMALIZED  DISPLACEMENT  OMEGAR  T I LDE* (DIM 

1  ENS I ONLESS ) - 

/////) 

A6 A  FORMAT ( 75H 

- NORMALIZED  DISPLACEMENT  U  T I LDE* ( D I  MENS  I 0 

INLESS) - 

/////) 

A65  FORMAT ( 75H 

- NORMALIZED  DISPLACEMENT  V  T I LDE* ( D I  MENS  I 0 

INLESS) - 

/////) 

END 

*  LABEL 

CFCTRL 

FUNCTION  FCTRL( INTGER) 
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D  FCTRL*1.0 

IF  (INTGER-1)  2  *  2  » 3 
3  DO  4  N«2» INTGER 
D  AN*N 

D  4  FCTRL»FCTRL*AN 
2  RETURN 
END 


*  LABEL 

CBER 

FUNCTION  BER  (X.LONG) 

D  PI»3. 141592653589793 

IF(X-17.0)5*5»8 
D  5  BER*1.0 

DO  16  K»l.LONG 
D  AKXK 

D  TOPBER-  (X/2.0)**(2.0*AK> 

D  TERBER»( < TOPBER/FCTRL < 2*K > )**2)*( <-1.0)**K) 

D  16  BER»BER+TERBER 

GO  TO  40 

D  8  COEFl»SQRTF(PI/(2.*X) ) *EXPF ( - ( X/SORTF ( 2 . ) ) ) 

D  COEF2*EXPF (X/SORTF ( 2  • > ) /SQRTF ( 2. *P I *X ) 

D  ALPHA*X/SORTF(2. J-PI/8. 

D  BETA«X/SQRTF(2. )+PI/8. 

D  FIRST*-SINF(BETA) 

D  SECOND=COSF( ALPHA) 

D  PROD*l. 

D  BER*COEF2*SECOND 

DO  53  1  *  LONG 

D  AK«K 

D  CCPSS*COSF( PI »AK/4. )*COSF( ALPHA ) +S iNF ( P I *AK/4 . )  »S I NF ( ALPHA ) 

D  PROD=PROD*(2.*AK-l.)»*2 

D  SCFM«PROD/FCTRLU> 

D  SCFMxSCFM/(8.*X)**(AK/2. > 

D  SCFM»SCFM/(8.*X)**(AK/2. ) 

D  SUM2«(SCFM*CCPSS)*COEF2 

IF(X-27.0)60*60»53 
D  60  SCFM*SCFM* (-1. )**K 

D  SCMCS— S INF (PI *AK/4. )*COSF( BETA) -COSF < P I *AK/4 . )*SINF( BETA) 

D  BER*BER-C0EF1*FIRST  /PI 

D  SUM  1 «  (SCFM*SCMCS )*C0EF1/PI 

D  53  BER*BER+SUM2-SUM1 
40  RETURN 
END 


*  LABEL 

CBEI 

FUNCTION  BE  I  (X.LONG) 
D  PI»3. 141592653589793 

IFIX-17.015.5.8 
D  5  BE  I *0  »  0 

DO  20  K*1 .LONG 
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TABLE  V  (Continued) 


D  AK  *K 

D  TPABEI » (X/2.0 >**< 2.0*AK ) 

D  TPBBEI»(X/2.0)**(2.0*AK-2.0> 

D  TER BEI= ( TPABEI /FCTRLI 2 *K-1 ) )  *  I TPBBE I /FCTRL I 2*K-1 ) ) *< < -1 .0 ) **  < K-l ) ) 

D  20  BEI-BEI+TERBEI 
GO  TO  40 

D  8  COEFl-SQRTFl PI / (2 .*X ) ) *EXPF<-( X/SQRTFI 2. >  )  > 

D  COEF2-EXPFIX/SQRTFI2. >  ) /SORT F ( 2 . *P I *X ) 

D  ALPHA-X/SQRTFI2. (-PI/8. 

0  BETA=X/SORTF(2.)+PI/8. 

D  FIRST*COSF (BETA) 

0  SECOND-SINFIALPHA) 

D  PROD= 1 • 

D  BEI*COEF2#SECOND 

DO  51  K  *  1  *  LONG 
D  AK*K 

D  CSMSC*COSF I  PI *AK/4. ) *S I NF ( ALPHA ) -S I NF ( P I *AK /4 . > *COSF ( ALPHA > 

D  PROD “PROD* (2«#AK-1.)**2 

D  SCFM-PROD/FCTRLIK) 

D  SCFM-SCFM/ <8.*X  >**< AK/2  .  ) 

D  SCFM-SCFM/ <8.*X)**<AK/2.  ) 

D  SUM2» ( SCFM*CSMSC > *COEF2 

IFIX-27. 0)60*60. 51 
D  60  SCFM=SCFM*(-1. )**K 

D  CCMSS*C0SF(PI*AK/4. >*COSF( BETA) -S INF (PI *AK/4. ) *S I NF ( BETA) 

D  SUM1  =  <SCFM*CCMSS)*CCEF1/PI 

D  BEI=BEI+C0EF1*FIRST/PI 

D  51  BEI=BEI+SUM1+SUM2 
40  RETURN 
END 


*  LABEL 

CAKER 

FUNCTION  AKER  (X.LONG) 

D  PI=3. 141592653589793 

D  GL2MIG=. 1159315156584124 

IFIX -8 .89)5*5*8 
D  5  SAME*GL2MIG-LOGF ( X ) 

D  PI4TH-  PI/4.0 

D  AKER*SAME 

DO  25  K«l.  LONG 
D  AK«K 

D  TOPBER=  (X/2.0)**(2.0*AK ) 

D  T  ERBER= ( ( TOPBER/FCTRL ( 2*K ) )**2>*( (-1.0) **K ) 

D  TPABEI«(X/2.0)**(2.0*AK> 

D  TPBBEI«(X/2.0)**(2.0*AK-2.0) 

D  TERBEI= I TPABEI /FCTRLI 2*K-1 ) ) *( TPBBEI /FCTRL ( 2«K-1 ) ) *  I I -1 .0 ) ** I K-l ) ) 

D  TERM1«SAME*TERBER 

D  TERM2»PI4TH*TERBEI 

I  «K-1 

IF  I  I )  7.6.7 
D  6  SUM1 = 1.5 
GO  TO  11 
D  7  PART *0.0 
57  DO  10  L=I »K 


71 


TABLE  V  (Continued) 


D  AL-L 

D  10  PART-PART+I  1.0/IAX+AL)  ) 

D  110  SUMl-SUMl+PART 
D  11  T0P3-I (X/2.0>**< 2.0*AX ) ) 

D  2 A  TERM3-I ( ( T0P3/FCTRL ( 2*X ) )**2) *SUM1*( ( -1 ,0)**X ) ) 

D  25  AXER-AXER+TERM1+TERM2+TERM3 
GO  TO  40 

0  8  C0EF1»SQRTF(PI/(2.*X) ) *EXPF ( -( X/SORTF ( 2 . ) ) ) 

0  ALPHA-X/SQRTFI2. ) -P I / 8 . 

D  BETA-X/SQRTFI 2. )+PI/8. 

D  FIRST  -COSF ( BETA ) 

D  PROD* 1 • 

D  AXER»COEFl*FIRST 

DO  50  X* 1  *  LONG 
D  AX-X 

D  PROD»PROD*( (2.*AX-1.)**2) 

D  CCMSS-COSFl PI*AX/4.)*C0SF(BETA)-SINF(PI*AX/4. )*S INF (BETA) 

D  SCFM*PROD/FCTRL (X ) 

D  SCFM-SCFM/ ( 8.*X )**( AX/2. ) 

D  SCFM«SCFM/(8.*X)**(AX/2.)*( (-l.)**X) 

D  SUM  1  *  (SCFM*CCMSS)*C0EF1 

D  50  AXER-AXER+SUM1 
40  RETURN 
END 


*  LABEL 

CAXEI 

FUNCTION  AXE  I  (X.LONG) 

D  GL2MIG-. 1159315156584124 

D  PI-3.141592653589793 

IFIX-8. 89) 5,5.8 
D  5  SAME-GL2MIG-L0GFIX ) 

D  PI4TH-PI/4.0 

D  AXEI »-P I 4TH 

DO  25  X* 1 .LONG 
p  AK-X 

D  TOPBER-  (X/2.0)**(2.0*AX) 

D  TERBER-I ( TOPBER/FCTRL ( 2*X )  )**2I* ( (-1.0)**X) 

D  TPABEI-(X/2.0)**(2.0*AX) 

D  TPBBEI-(X/2.0)**(2.0*AX-2.0) 

D  TERBEI-I TPABEI /FCTRL(2*X-1 ) )*( TPBBE I /FCT RL ( 2*X-1 ) )*( (-1.0)**  <  X  — 1 ) ) 

D  TERM1»SAME*TERBEI 

D  TERM2-PI4TH*TERBER 

I  *  X-l 

IF  ( I )  86.86,88 
D  86  SUM2-1.0 
GO  TO  80 
D  88  PART-0.0 

DO  94  L-I .X 
D  AL-L 

D  94  PART -PART+ ( 1 . 0/ ( AX+AL- 1 . 0 ) ) 

D  SUM2-SUM2+PART 

D  80  TERM3-SUM2* ( (-1.0)**(X-1) )*( ( (X/2.0)** ( 2.0* AX)  l/FCTRL ( 2*X-1  ))*((( X 
D  2/2.0 )**( 2.0* AX-2.0) ) /FCTRL ( 2*X-1 ) ) 

D  25  AXEI -AXEI+TERM1-TERM2+TERM3 
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TABLE  V  (Continued) 


GO  TO  40 

D  8  C0EF1*SQRTF ( PI / ( 2 .*X ) ) *EXPF ( -X/SQRTF ( 2. ) ) 

0  ALPHA  =  X/SQRTF( 2. ) —  P I / 8 • 

D  BETA=X/SORTF( 2. )+PI /8. 

D  FIRST* -SINF( BETA) 

D  PROD  *  1 • 

D  akei =coefi*first 

DO  52  K«1 .LONG 
D  AK  =  K 

D  SCMCS*-SlNF(PI*AK/4.)*C0SF(BETA)-C0SF(PI*AK/4.  )  *S  INF  (BETA) 

D  PROD  =  PROD* (  ( 2 • *AK - 1 • )**2) 

D  SCFM*PROD/FCTRL( K ) 

D  SCFM=SCFM/(0**X)**(AK/2.) 

D  SCFM=SCFM/(8**X )  ** ( AK / 2 • )  *  ( (  - 1  • )  **IC ) 

D  SUM  1  *  ( SCFM*SCMCS )*C0EF1 

D  52  AXEI =AKEI+SUM1 
40  RETURN 
END 


*  LABEL 

CBRP 

FUNCTION  BRP  (X.LONG) 

D  PI=3. 141592653589793 

IF (X-17.CU5.5.8 
D  5  BRP*  0  *0 

DO  65  K * 1 .LONG 
D  A  K  =  K 

D  P  A  R  T  A  *  (  (X/2*0)**(2*0*AK)  )/FCTRL(2*K-l) 

D  PARTB*((X/2*0)**(2.0*AK-1„0)) /FCTKL ( 2*K  ) 

D  TERBRP*PARTA*PARTB*( (-1.0 )**K ) 

D  65  BRP=BRP+TERBRP 
GO  TO  40 

D  8  C0EF1=-(SQRTF(PI/(2.*X>  )  *{fXPF  (  -X/SQRTF  (  2  .  )  )  ) 

D  COEF2*EXPF ( X/SQRTF ( 2.  )  ) /SORT F ( 2 • *P I *X ) 

D  ALPHA* X/SQRTF ( 2 • ) -PI /8 • 

D  BETA*  X/SQRTF(2. J+PI/8. 

D  FIRST*-SINF( ALPHA) 

D  SECOND*COSF( BETA) 

D  BRP=COEF2*SECOND 

D  PROD  =  1*0 

DO  55  K  =  1  *  LONG 
D  AK*K 

D  SSMCC*-S INF ( PI*AK/4. )*SINF( BET  A ) -COSF ( P I  * AK /4 • )#COSF (BETA) 

D  J=K-1 

D  A  J  =  J 

D  I F  C JJ70.70.71 

D  70  PROD  *  1*0 
GO  TO  72 

D  71  PROD«PROD*(2.*AJ-l.)**2 
D  7  2  SCFM  =PROD* ( 4. *AK**2-1 • ) /FCTRL  UC  ) 

D  SCFM*SCFM/( 8.*X )**( AK/2. ) 

D  SCFM*SCFM/ (8.*X>** ( AK/2. ) 

D  SUM2*SCFM*SSMCC*COEF2 

IF (X“27«0)60*60*55 
D  60  SCFM=SCFM* ( -1 . ) **K 
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T  A  HLK  V  (('nntinut.'ri) 


D  SCPCS=SINF (Pl*AK/4. ) *COSF ( ALPHA ) +COSF ( P I *AK /4 . )*S INF (ALPHA ) 

D  SUM1=SCFM*SCPCS*F IRST/P I 

D  BRP=BRP-C0EF1*F IRST/P I 

D  55  BRP=BRP-SUM1+SUM2 
40  RETURN 
END 


*  LABEL 

CBIP 

FUNCTION  BIP  ( X  t  LONG ) 

D  P  I  =  3  •  141592653589753 

IF (X-17.0)5,5f8 
D  5  B I P  =  0 • 0 

DO  68  K  =  1  *  LONG 
D  AK*K 

IF  ( K - 1 )  66  *  66 ♦ 67 

D  66  PARTC- ( 1 .0/ ( X/2 .0 )  ) /FCTRL ( 2*K-1  ) 

D  PARTD=(  (X/2#0)**(2«0#AK)  ) 

GO  TO  74 

D  67  PAR  TC  * ( (X/2«0)**(2»0*AK-3»0) )/FCTRLt2*K“l) 

D  73  PARTD=(  (X/2.0)**(2.0*AK)  )/FCTRL(2*K“2) 

D  74  TERBIP=PARTC*PARTD*( (-1 .0)**( K-l ) ) 

D  68  BIP=BIP+TERBIP 
GO  TO  40 

D  8  COEFls~(SQRTF(PI/(2**X)  )*EXPF(-X/SQRTF(2«)  )  ) 

D  COEF  2  *EXPF ( X/SQRTF ( 2. )  ) / SORT  F ( 2.*PI *X ) 

D  ALPHA-X/SQRTF(2. )-PI/8. 

D  BET  A  =  X/ $OR  T  F ( 2 • )+PI/8. 

D  FIRST=COSF( ALPHA) 

D  SECOND=S I NF ( BETA ) 

D  BIP=COEF2*SECOND 

D  PROD* 1 • 0 

DO  57  K=1 » LONG 
D  AK*K 

D  SCMC5  =  S INF ( PI *AK/4. ) *CO$F (BETA) -COSF ( P I *AK /4  • ) *SI NF ( BET  A ) 

D  J*K-1 

D  A  J  =  J 

D  IF (J)70t70*71 

D  70  PROD  =  1 • 0 
GO  TO  72 

D  71  PROD=PROD*(2.*AJ-l.)**2 
D  72  SCFM«PROD*(4«*AK**2-l. )/FCTRL(K) 

D  SCFM=SCFM/(8.*X)**(AK/2*) 

D  SCFM*SCFM  /  (8«#X)**(AK/2«  ) 

D  SUM2*SCFM*SCMCS*COEF2 

IF(X-27*0)60»60f57 
D  60  SCFM=SCFM*(-1. )**K 

D  SSMCC  *S  I  NF ( PI*AK/4. ) *S I NF ( A L PH A ) -COSF ( P I  * AK /4 .  ) *COSF (ALPHA ) 

D  SUM1*5CFM*S5MCC*C0EF1/PI 

D  B I P -B I P+F IRST*C0EF1/PI 

D  57  BIP=BIP+SUM 1 +  5UM2 
40  RETURN 
END 
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TABLE  V  (Continued) 


*  LABEL 

cakerp 

FUNCTION  AKERP(X.LONG) 

D  GL2MIG*.l 159315156584124 

D  PI=3. 141592653589793 

IFIX-8.89 J5.5.8 
D  5  SAME*GL2MIG-LOGF(  X ) 

D  PI 4TH»P 1/4.0 

D  CONST  =  1  • 0/X 

D  AKERP “-CONST 

DO  25  X  *  1  *  LONG 
D  Ak  =  K 

D  TOPBER=  (X/2.0)**(2.0*AX) 

D  T  ERBE  R= ( ( TOPBER / FCTRL < 2*K ) )**2)*( (-1.0)  **K  ) 

D  PART  A  * ( (X/2.0) **(  2.0*AK ) ) /FCTRL ( 2*K_1  ) 

D  PARTB»< (X/2.0 >**(2.0*AK-1.0) ) /FCTRL <2*K ) 

D  TERBRP=PARTA*PARTB*( (-1 .0) **M 

IF  ( K - 1 )  66.66.67 

D  66  PARTC=( 1 .0/ ( X/2.0  )) /FCTRL (2*X-1  ) 

D  PARTD= ( (X/2.0 ) **< 2.0*AX  )  ) 

GO  TO  74 

D  67  PARTC= ( (X/2.0 )**( 2. 0*AK-3.0) ) /FCTRL (2*X-1 ) 

D  73  P  AR  T  D“  ( (X/2.0) **( 2 • 0*AK  )  )/FCTRL(2*K“2) 

D  74  TERBIP  =  PARTC*PARTD*(  (- 1 *0 ) *» I K-l  )  ) 

D  TERM1*SAME*TERBRP 

D  TERM2 “CONST* TERBER 

D  TERM3  =  PI4TH*TERB  IP 

I  “K  — 1 

IF  ( I )  7.6.7 
D  6  SUM  1*1.5 

GO  TO  11 
D  7  PART  =  0.0 

57  DO  10  L“  I  *  K 

D  AL  =  L 

D  10  PART=PART  +  ( 1 .0/ ( AX  +  AL )  ) 

D  110  SUM1=SUM1+PART 

D  1 1  TERM4»SUM1*( 2.0*AK )* ( ((X/2.0)**(2.0*AX))/FCTRL(2*K))*(((X/2.0)**(2 
D  Z  .0*AK-1 .0 )  )  /  FCTRL ( 2*K)  )*  (  (-1 .0) **K ) 

D  25  AKERP-AKERP+TERM1-TERM2+TERM3+TERM4 
GO  TO  40 

D  8  COEF1  —  (SORTFIPI/ (2.*X)  )  * EXPF < -X/ SQR TF < 2 . )  )) 

D  ALPHA«X/SQRTF(2. l-PI/8. 

D  BETA=X/SQRTF(2. l+PI/8. 

D  FIRST-COSF(ALPHA) 

D  AKERP«C0EF1*FIRST 

D  PROD“ 1 . 0 

DO  54  K«1 .LONG 
D  AK“K 

D  SSMCC-SINFI PI*AK/4. ) *S I NF ( ALPHA ) -COSF ( P I *AK /4 . ) *COSF ( ALPHA ) 

D  J  =  K-1 

D  A  J  =  J 

I  F  t JI70.70.71 
D  70  PROD= 1 . 0 
GO  TO  72 

D  71  PROD“PROD*( 2.*AJ-1 . ) **2 
D  72  SCFM“PROD*(4.*AK**2-l . ) /FCTRL ( K ) 
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TABLE  V  (Continued) 


D  SCFM-SCFM/(8.*X)**(AX/2.) 

D  SCFM-SCFM/(8.*X)**(AX/2.  )*(-l. )**X 

D  SUM1»SCFM*SSMCC*C0EF1 

D  54  AXERP-AXERP+SUM1 
40  RETURN 
END 


*  LABEL 

CAxEIP 

FUNCTION  AXE  I P ( X  » LONG ) 

D  GL2MIG-. 1159315156584124 

D  PI-3.141592653589793 

I F ( X-15.5 ) 5 » 5  » 8 
D  5  SAME-GL2MIG-LOGF (X ) 

D  CONST-1. 0/X 

D  PI4TH-PI/4.0 

D  AXEIP-0.0 

DO  25  X - 1 .LONG 
D  AX-X 

D  TPABEI«(X/2.0)**(2.0*AX) 

D  TPBBEI«(X/2.0)**(2.0*AX-2.0) 

D  TERBE1- (TPABEI/FCTRL(2*X-1 ) )  *  I TPBBE I /FCTRL I 2*X- 1 ) ) *  I  I -1 .0 ) ** I X-l )  ) 

D  PART  A- ( (X/2.0)**(2.0*AX ) ) /FCTRL ( 2*X-1 ) 

D  PARTB- ( (X/2.0 )**( 2.0*AX-1 .0) ) /FCTRL ( 2*X  ) 

D  TERBRP-PART A*PARTB* ( l-1.0)**X) 

IF  (X-l)  66 *66.67 

D  66  PARTC-I 1 .0/1 X/2.0) ) /FC TRL ( 2*X- 1 ) 

D  PARTD-I (X/2.0)**(2.0*AX)  ) 

GO  TO  74 

D  67  PARTC-I (X/2.0)**(2.0*AX-3.0) ) /FCTRL < 2*X- 1 ) 

D  73  PARTD-I (X/2.0)**(2.0*AX) ) /FCTRL ( 2*X-2 ) 

D  74  TERBIP  =  PARTC*PARTD*I I - 1 . 0 ) ** I X-l  )  ) 

D  TERM1-SAME*TERBIP 

D  TERM2»CONST*TERBE I 

D  TERM3-P I 4TH*TERBRP 

I  -X-l 

IF  (I)  86.86.88 
D  86  SUM2-1.0 
GO  TO  80 
D  88  PART-0.0 

DO  94  L- I . X 
D  AL-L 

D  94  PART-PART+I 1.0/1 AX+AL-1.0) ) 

D  SUM2-SUM2+PART 

D  8  0  TERM4«SUM2*I2.0*AX-1.0)*I  I  - 1 . 0  )  **  I X- 1 )  )  *  I  I  I  X /2 .0  )  **  I  2 . 0*AX  )  ),' FCTRL 
D  2  I 2*X- 1 ) >*( I IX/2.0)**(2.0*AX-3.0) ) /FCTRL I 2*X-1  )  ) 

D  25  AXEIP«AXEIP+TtRMl-TERM2-TERM3+TERM4 
GO  TO  40 

D  8  COEF1—  (SORTF(PI/(2.*X  )  )  *EXPF  I -X/SQR  TF  I  2  .  ))) 

D  ALPHA-X/SQRTF 12. I-PI/8. 

D  BETA-X/SQRTFI2. I+PI/8. 

D  FIRST--SINF(ALPHA) 

D  PROD- 1 . 0 

D  AXEIP»C0EF1*FI RST 

DO  56  X-l. LONG 
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TABLE  V  (Continued) 


0  A  K  *K 

0  SCPCS=SINF ( PI*AK/4. ) *COSF ( ALPHA > +COSF ( P I *AK /4 . }*SINF (ALPHA } 

0  A  J  -  J 

D  J*K-1 

IF(J)70 *70*71 
D  70  PROD* 1*0 
GO  TO  72 

0  71  PROD*PROD* ( 2.*A J-l • ) **2 

D  72  SCFM*PROD*(4.*AK**2~l. )/FCTRL(K) 

0  SCFM*SCFM/ (8.*X)**(AK/2.) 

0  SCFM=SCFM/ (8.*X)**(AK/2.)*(-l.)**K 

D  SUM1*SCFM*SCPCS*C0EF1 

D  56  AKEIP*AKEIP+SUM1 
40  RETURN 
END 
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1.  Input 

The  input  to  our  program  as  indicated  by  step  (1)  and  step  (8)  in  the  master  flow  chart  con¬ 
sists  of  a  number  of  records  prestored  on  machine  tape  A2. 

Record  1  Control  Parameters  I  (415) 

This  card  contains  four  non-negative  fixed-point  variables. 

I  M  I  MU  I  ML  I  MC  I 

These  parameters  control  the  number  and  spacing  of  the  y-tabular  points  between  R3  and  R4  (to 
be  defined  later  under  Record  4). 

Total  number  of  y-tabular  points  between  R3  and  R4. 

Number  of  (evenly  spaced)  intervals  in  the  upper  edge  zone 
(a  region  2(2fh )^/2  from  R4). 

Number  of  (evenly  spaced)  intervals  in  the  lower  edge  zone 
(a  region  2(2fh)1/^  from  R3). 

A  control  parameter  to  be  set  equal  to  0,  1,  or  2. 

MC  =  0  y-tabular  points  are  evenly  spaced. 

MC  =  1  Tabular  points  are  spaced  differently  (usually 
denser)  only  in  the  upper  edge  zone.  In  par¬ 
ticular,  there  will  be  MU  +  1  tabular  points 
in  this  edge  zone.  This  option  is  used  prima¬ 
rily  when  the  shell  is  closed  at  the  apex  so  that 
no  edge  effect  appears  in  the  lower  edge  zone. 

MC  =  2  Tabular  points  are  spaced  differently  (usually 
denser)  in  both  edge  zones.  There  are  MU  +  1 
and  ML  +  1  tabular  points  in  the  upper  and  lower 
edge  zones  respectively. 

Record  2  Control  Parameters  II  (515) 

This  card  contains  five  non-negative  fixed-point  variables 

I  N  I  NORM  I  K1  I  K2  I  K3  I 

N  (<:  10)  Number  of  equally  spaced  0-tabular  points  between  THETA3  and 
THETA 4  (to  be  defined  later  under  Record  5)  along  any  fixed  lati¬ 
tude  of  the  shell.  If  only  the  superscripted  quantities  in  (11-85) 
and  (11-86)  are  desired,  N  must  be  set  equal  to  zero  and  Records 
6  and  7  omitted. 

NORM  Control  parameter  specifying  whether  normalized  results  are 

desired. 

NORM  =  0  Normalized  results. 

NORM  >  0  Un-normalized  results. 


M  «50) 
MU 

ML 

MC 
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Control  parameter  specifying  the  edge  constraints  at  R1  (to  be 
defined  later  under  Record  4). 

K1  =  1  Shell  is  free  at  Rl. 

K1  2  Shell  is  clamped  at  Rl. 

K1  -  3  Shell  is  simply  supported  at  Rl. 

K1  4  Shell  is  closed  at  the  apex. 

Control  parameter  specifying  the  edge  constraints  at  R2  (to  be 
defined  later  under  Record  4). 

K2  -  1  Shell  is  free  at  R2. 

K2  -  2  Shell  is  clamped  at  R2. 

K2  =  3  Shell  is  simply  supported  at  R2. 

Different  combinations  of  PS1  and  MGOTO  (see  Record  8)  may  be 
used  for  each  structure.  If  more  than  one  such  combination  is 
to  be  run  for  the  same  structure,  only  Record  8  needs  to  be 
changed.  K3  then  specifies  the  total  number  of  such  combinations 
for  the  same  structure.  Set  K3  =  0,  if  the  only  run  wanted  is  for 
$  =  0. 

Record  3  Material  Parameters  (3E20.8) 

This  card  contains  three  E-type  floating-point  variables. 

I  RHO  I  E  I  PRT  I 

RHO  Volume  weight  density  of  the  shell  (lb/in.  ). 

E  Young’s  modulus  (lb/in.  ). 

PRT  Poisson’s  ratio. 

Record  4  Geometrical  Parameters  I  (4F15.9). 

This  card  contains  four  F-type  floating-point  variables. 

I  Rl  I  R2  I  R3  I  R4  I 

Rl  (in.)  Value  of  r  at  the  lower  edge  of  the  shell.  If  the  shell  is  closed 
at  the  apex,  then  Rl  =  0. 

R2  (in.)  Value  of  r  at  the  upper  edge  of  the  shell. 

R3  (in.)  Smallest  value  of  r  at  which  stresses  and  displacements  are  to 

be  calculated  (smallest  y-tabular  point).  Generally,  this  is  the 
same  as  Rl.  However,  if  the  shell  is  closed  at  the  apex,  we 
set  R3  >  0  to  avoid  possible  complications  arising  from  calculat¬ 
ing  the  desired  quantities  at  the  apex. 

R4  (in.)  Largest  value  of  r  at  which  the  output  is  to  be  given  (i.e,,  the 
largest  y-tabular  point).  Generally,  this  is  the  same  as  R2. 


K 1 


K2 


K3 
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Record  5  Geometrical  Parameters  II  (4F15.9) 

This  card  contains  four  F-type  floating-point  variables. 

I  THETA3  I  THETA4  I  F  I  H  I 


THETA3  Smallest  value  of  0  in  fractions  of  7r  at  which  output  is  desired 
(the  smallest  0-tabular  point). 

THETA4  Largest  ©-tabular  point  in  fractions  of  7r  at  which  output  is 
desired. 


F  Focal  length  of  the  paraboloidal  surface  (in.). 

II  Shell  thickness  (in.). 


Record  6  Variable  Format  Statement  I(12A6) 

This  card  provides  a  format  statement  for  the  set  of  G-tabular  points,  for  example, 

I  19H  R(IN.)/THETA(DEG.),  NF11.2  I 

Record  7  Variable  Format  Statement  II  (12A6) 

This  card  provides  a  format  statement  for  the  dependent  variables  in  the  output,  for 
example, 

I  F19.6,  NF11.6  I 


Record  8  Control  Parameters  III  (F15.9, 15) 

This  card  contains  one  F-type  floating-point  variable  and  one  fixed-point  variable. 

I  PSI  I  MGOTO  I 


PSI 

MGOTO 


The  pointing  angle  ip  in  fractions  of  it. 


Control  parameter  specifying  the  type  of  output  desired. 

MGOTO  =  0  Stresses  and  displacements  for  symmetric 
or  antisymmetric  deformations  or  the 
portion  of  the  results  independent  of  0 
are  given  in  a  concise  form. 

Only  the  zero-superscripted  tilde  quantities 
(and  the  corresponding  w^ex  displacement) 
are  given. 

Usual  set  of  structural  quantities  are  given. 
Only  the  tilde  displacement  quantities  (and 
the  corresponding  displacement)  are 

given. 

Only  the  root-mean-square  of  the  phase 
error*  is  given. 


MGOTO  -  1 


MGOTO  =  2 
MGOTO  =  3 


MGOTO  -  4 


*To  be  discussed  in  a  subsequent  chapter. 
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2.  Output 

The  output  to  the  shallow  shell  program  is  arranged  in  a  manner  very  similar  to  that  of  the 

membrane  analysis.  The  first  part  of  the  output  reproduces  the  input  to  the  program.  It  states 

the  problem,  the  edge  constraints,  the  various  geometrical  and  material  parameters  and  whether 

the  numerical  results  have  been  normalized.  The  second  part  of  the  output  gives  the  value  of  the 

physical  quantities  for  the  net  of  tabular  points.  If  4  =  0,  the  behavior  of  the  shell  is  axisymmet- 

rio.  For  this  case,  we  have  v  =  N  „  =  M  ~  =  0,  while  u,  w,  oj  ,  N  ,  N^,  M  ,  M~,  and 

r0  0  r0  r  r  0  r  0 

CJ  are  all  functions  of  y  only.  The  output  in  this  instance  is  given  in  a  concise  form  as  in  the 
membrane  program.  For  0  <  4  <  0.5,  we  have  several  options  of  output  which  are  very  similar 
to  those  of  the  membrane  program: 

(1)  List  only  the  zero-superscripted  quantities  for  the  set  of  y-tabular 
points  with  the  displacement  quantities  being  those  relative  to  the 
undeformed  shell  (set  MGOTO  =  0). 

(2)  List  only  the  zero-superscripted  displacement  quantities  u,  v,  and 
w  relative  to  those  of  the  face-up  position  (set  MGOTO  =  1). 

(3)  List  all  the  actual  stresses  and  displacements  for  the  entire  net  of 
tabular  points  with  the  displacement  quantities  being  those  relative 
to  the  undeformed  shell  (set  MGOTO  -  2). 

(4)  List  only  the  actual  displacements  for  the  net  of  tabular  points  relative 
to  those  of  the  face-up  position  (set  MGOTO  =  3). 

(5)  List  only  the  root-mcan-square  of  the  phase  error*  (set  MGOTO  -  4). 

3.  Operational  Information 

The  program  is  coded  in  FORTRAN  language  for  a  32K  IBM  7090  (and  7094)  and  is  compiled 
by  a  FORTRAN  II  compiler.  It  requires  9  FORTRAN  functions,  which  calculate  the  values  of 
the  Thomson  functions  and  their  first  derivatives  using  double-precision  arithmetic  as  described 
in  IBM  Bulletin  J28-6114  "32K  709/7090  FORTRAN:  Double-Precision  and  Complex  Arithmetic.” 
All  other  routines  used  appear  as  library  routines  on  the  Lincoln  Laboratory  library  tape. 

Some  of  the  subscripted  variables  and  their  analytical  counterparts  are  presented  in  Table  VI. 
It  should  be  understood  that  if  normalized  results  are  requested,  affected  expressions  in  the 
third  column  of  Table  VI  should  be  replaced  by  the  corresponding  starred  quantities. 

C.  EXAMPLES 

1.  Axisymmetric  Deformations 

Consider  a  shell  in  the  face-up  position  [4  =  0  °)  which  is  closed  at  the  apex  and  simply 
supported  at  the  upper  edge  r^  with 


*To  be  discussed  in  a  subsequent  chapter. 
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TABLE  VI 

SOME  SUBSCRIPTED 

VARIABLES 

AND 

THEIR 

ANALYTICAL 

COUNTERPARTS 

Subscripted  Variable 

Dimension 

Analytical  Expression 

W 

50  x  10 

w 

WP 

50  X  10 

3w 

3r 

V 

50  X  10 

V 

u 

50  X  10 

u 

CMR 

50  X  10 

M 

r 

CMT 

50  X  10 

Me 

CMRT 

50  X  10 

Mre 

CNR 

50  X  10 

N 

r 

CNT 

50  X  10 

Ne 

CNRT 

50  X  10 

Nre 

QR 

50  X  10 

QT 

50  X  10 

Qe 

WFLEX 

50  X  10 

Wri 

flex 

P 

50 

s 

w 

W1 

50 

,  3w  vs 
{  Or  } 

Q 

50 

s 

u 

EMR 

50 

M  s 
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EMT 

50 

Me 

ENR 

50 

N  S 
r 

ENT 

50 

Ne 

Qi 

50 

QrS 

PWB 

50 

S 

wr, 

flex 

ww 

50 

a 

w 

W3 

50 

/  cftv  »a 

1  3r  ’ 

vv 

50 

a 

V 

uu 

50 

a 

u 

DMR 

50 

M  a 
r 
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TABLE  VI  (Continued) 

Subscripted  Variable 

Dimension 

Analytical  Expression 

DMT 

50 

Moa 

DMRT 

50 

DNR 

50 

N  a 
r 

DNT 

50 

Ne 

DNRT 

50 

Nre 

QQR 

50 

Qra 

QQT 

50 

Qo 

WA5 

50 

a 

wr, 

flex 

WT 

50 

w 

WPT 

50 

0W 

3r 

UT 

50 

u 

V2 

50 

V 

WFLEXT 

50 

wrl 

flex 

R 

50 

r 

THETA 

20 

0 

BRER 

50 

ber 

BBEI 

50 

bei 

A  AKER 

50 

ker 

AAKEI 

50 

kei 

RBRP 

50 

ber' 

BBIP 

50 

bei* 

AAKERP 

50 

ker' 

AAKEIP 

50 

kei' 

SBER 

50 

ber" 

SBEI 

50 

bei" 

SKER 

50 

ker" 

SKEI 

50 

kei" 

TBER 

50 

ber'" 

TREI 

50 

bei"' 

TKER 

50 

ker"' 

TKEI 

50 

kei"’ 
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rho  =  0.1  lb/ in? 

E  107  lb/in.2 
v  =  0.3 
r  j  =  0  in. 

=  300  in. 
f  =  450  in. 
h  =  1  in. 

Since  the  deformation  is  axisymmetric,  we  set  N  0  and  omit  Records  6  and  7.  For  the  purpose 
of  illustration,  we  consider  twelve  y-tabular  points  with  a  finer  grid  in  the  upper  edge  zone 
(MC  =  1).  As  the  shell  is  closed  at  the  apex,  we  avoid  calculating  the  limiting  value  of  the  output 
at  the  apex  by  setting  R3  >  0.  The  results  will  not  be  normalized,  so  we  get  the  actual  physical 
quantities.  Output  will  be  the  zero-superscripted  quantities  for  the  set  of  y-tabular  points,  where 
the  displacement  quantities  are  those  relative  to  the  undeformed  shell. 


INPUT 


Record  1  Control  Parameters  I  (415) 
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Record  2  Control  Parameters  II  (515) 
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Record  3  Material  Parameters  (3E20.8) 
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Record  4  Geometrical  Parameters  I  (4F15.9) 
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Recard  5  Geometrical  Parameters  II  (4F15.9) 
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Remark:  Since  N  =  0,  Records  6  and  7  have  been  omitted  according  to  input  instructions. 
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Record  8  Control  Parameters  III  (FI 5.  9,  15) 
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The  computer  output  is  presented  in  Table  VII. 

2.  Unsymmetric  Deformations 

Consider  a  shell  which  is  clamped  at  the  lower  edge  r ^ ,  free  at  the  upper  edge  r^,  and 
oriented  in  such  a  way  that  ip  ^  0.  We  have 

rho  =  0.1  lb/ in? 

E  -  107  lb/in  2 

v  -  0.3 

r ,  =  100  in. 

1 

r^  =  300  in. 
f  =  500  in. 
h  =  1  in. 

We  now  want  normalized  output  for  seven  0-tabular  points  in  the  interval  — tt/2  <  0  <  7t/2  and 
five  y-tabular  points  in  the  interval  y  <  y  <  y As  far  as  output  is  concerned,  we  will  list  only 
the  actual  displacements  for  the  net  of  tabular  points  relative  to  those  of  the  face-up  position. 
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INPUT 


Record  1  Control  Porometers  I  (415) 
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Record  2  Control  Parameters  II  (515) 
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Record  3  Material  Parameters  (3E20.8) 
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Record  4  Geometrical  Parameters  I  (4F15.9) 
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Record  5  Geometrical  Parameters  It  (4F15.9) 
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Record  6  Variable  Format  Statement  I-FMT1-(12A6) 


(18H  GAMMA/THETA(DEG.),  7F11.2) 
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|  777777  7777777777777  777  7777777777777777777777777777777777777777777777777777777 

8  8  88888888  888  888  8888  8  88)8888888888888888888888888888888888888888888888 
99999999999999999999999999999999939999999999999999999999999999999999999999999999 
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Record  7  Variable  Format  Statement  ll-FMT2-(12A6) 


(FI 8. 4, 7F1 1.6) 

0  00000  000000000000000000000000000000000000000000000000000000000000000000000000 

1  2  3  4  1  *  7  •  I  1*  n  n  13  M  B  n  17  tt  n  »  21  71  n  2«  K  N  n  N  H  *  31  M  13  W  JS  *  J7  M  J«  *0  «'  «2  U  U  «  W  «•  M  M  J1  J?  U  54  ss  5«  r  M  if  *  »i  m  u  m  n  ••  r  n  n  n  7?  w  *  n  n  n  r»  H 
111  111111  11111111111111111111111111111111111111111111111111111111111111111111 

2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2 

33333  3  3333  3333333333333333333333333333333333333333333333333333333333333333333 

4  4444  4444444  44444444444444444444444444444444444444444444444444444444444444444 

55555555555555555555555555555555555555555555555555555555555555555555555555555555 
66  666666  666  666666666666666666666666666666666666666666666666666666666666666666 

8  88  8  8888  8  88888888888888888888888888888818888888888888888818888818188818888 
99999999999999999999999999999999999999999999999999999999999999999999999999999999 

1  ?  3  4  J  (  7  »  IIIM  17  a  M  »  H  17  M  !»  TO  71  77  23  71  »  «  77  7*  71  30  31  J2  33  W  M  X  J7  M  31  40  <1 4?  43  44  4S40 17  41 14  50  J1  S2  1J  W  M  30  *7  M  MM  »1  *2  03  14  03  04  *7  M  M  «  71  77  77  H  TIM  77  71  77  M 

ni"iOiD 


Record  8  Control  Parameters  III  (FI 5. 9,  15) 


16666667 


180000000000000 
i  2  3  n  t  miMinDMB 
111111  11111111 

222222222222222 

33333  333333333 

4444444444444444 

555555555555555 

6666666  66 

7777777777777  7 

18888  888888888 


999999999999999 

1  7  3  4  J  4  7  f  t  to  71  12  13  W  IS 

mvToFn  _ 


M 

G 

O 

T 

O 

3 


0  0  0  0  0 

N  77  II  If  20 

11111 

2  2  2  2  2 

3  3  3  3 

4  4  4  4 

5  5  5  5  5 

6  6  6  6  6 

7  7  7  7  7 

1  8  8  8  8 


9  9  9  9  9 

17  «  If  70 


000000000000000000000000000000000000000000000000000000000000 
lb  77  23  74  n  71  77  71  71  30  31  37  33  34  33  31  37  31  31  W  41  <7  <3  44  45  U  47  41  4f  30  St  12  53  51  IS  SI  37  M  Jf  M  fl  12  U  W  13  M«7  H  If  70  71  7?  77  74  75  71  77  71  71  M 

111111111111111111111111111111111111111111111111111111111111 

222222222222222222222222222222222222222222222222222222222222 

333333333333333333333333333333333333333333333333333333333333 

444444444444444444444444444444444444444444444444444444444444 

555555555555555555555555555555555555555555555555555555555555 

666666666666666666666666666666666666666666666666666666666666 

777777777777777777777777777777777777777777777777777777777777 

888888888888888888888888848888868888888888888888888884888888 

999999999999999999999999999999999999999999999999999999999999 

ITl2I1174?S2l72nn3O  31  33  33  34  3S3l  37  3l3ll0  41  47  43444SI4  47  4l4l30lt5333S43SSIS7  50  5ll00lfl7l3l4ISIM?0l0l70nn77?4n7l  77  7Onie 


The  computer  output  is  presented  in  Table  VIII. 
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LINCOLN  LABORATORY  ANALYSES 
OF  PARABOLOIDAL  SHELLS 
(LLAPS) 


USER'S  MANUAL 

I.  INTRODUCTION 

The  quest  for  more  and  more  precise1  radars  in  the  ultra-high  frequency  regime  has  im¬ 
posed  stringent  requirements  on  the  structural  behavior  of  large  antennas.  The  tolerance 
requirements  for  an  antenna  surface  at  10,000  Mcps,  for  instance,  are  generally  set  at  about 
l/l6  of  the  operating  wavelength.  At  a  frequency  of  10,000  Mcps,  this  is  3/ 1 6  of  a  centimeter 
or  0.074  inch.  Such  a  tolerance  on  the  permissible  distortions  of  a  structure  which  may  be  a 
hundred  feet  or  more  in  over-all  size,  and  which  assumes  different  orientations  with  respect  to 
the  direction  of  gravity,  requires  an  extremely  high  degree  of  sophistication  in  analysis,  design, 
and  construction. 

Structures  such  as  bridges,  buildings,  and  flight  vehicles  are  designed  mainly  by  strength 
considerations,  although  flight  vehicles  rrust  also  have  a  certain  minimum  stiffness  in  order  to 
avoid  aero-elastic  difficulties.  Machine  lools  must  possess  great  stiffness,  but  machine  tools 
are  generally  compact  and  weight  limitations  are  relatively  unimportant.  On  the  other  hand,  the 
primary  design  requirement  of  a  high-performance  antenna  is  that  the  reflecting  surface  remain 
paraboloidal  and,  in  the  case  of  an  antenna  housed  in  a  radome,  strength  considerations  play  a 
minor  role  in  the  design.  The  antenna  must  therefore  have  adequate  structural  stiffness  but, 
since  the  predominant  loads  are  its  own  dead  weight,  the  structural  stiffness  should  be  accom¬ 
panied  by  minimum  weight;  that  is,  the  antenna  design  should  maximize  the  ratio  of  structural 
stiffness  to  weight. 

The  basic  structural  components  of  the  antenna  are  paraboloidal  panels  which,  when  joined 
together,  form  a  surface  of  revolution.  Such  a  structural  configuration,  if  properly  fabricated, 
can  be  considered  as  a  shell.  The  calculation  of  stresses  for  the  strength  design  of  a  shell  can 
be  achieved  with  a  fair  degree  of  confidence,  since  confidence  in  the  integrity  of  a  structure  can 
be  obtained  by  increasing  the  factor  of  safety,  i.  e.  ,  by  putting  more  material  into  the  structure. 
Such  a  course  of  action  may  be  self-defeating  in  an  antenna  which  has  only  to  resist  its  own  dead 
weight.  Moreover,  the  determination  of  the  shape  of  the  antenna  surface  must  be  precise,  and 
cannot  be  approached  with  the  same  philosophy  which  is  attendant  to  a  strength  design,  i.  c.  , 
hidden  under  a  factor  of  safety. 

Shell  structures  derive  many  of  their  attractive  features  from  their  two-dimensional  surface 
nature.  This  two-dimensional  nature  brings  with  it  geometrical  complications  to  the  strain- 
deflection  relations  and  the  equilibrium  equations.  Although  the  deflections  of  trusses,  beams, 
and  space  frameworks  are  well  understood,  well  documented,  and  easily  obtained,  the  exact 
opposite  is  true  for  shells,  as  evidenced  by  the  fact  that  shell  behavior  is  currently  the  major 
topic  of  study  in  the  field  of  structural  mechanics.  The  available  solutions  for  even  simple 
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loadings  of  simple  shells  are  of  such  a  form  that  numerical  results  are  not  easily  obtained.  For 
these  reasons,  Lineoln  Laboratory  has  been  actively  studying  the  deformations  of  paraboloidal 
shells.^  ^ 

This  user's  manual  will  describe  the  capability,  potentiality,  and  idiosyncrasies  of  the 
various  LLAPS  (Lineoln  Laboratory  Analyses  of  Paraboloidal  Shells)  computer  programs  whieh 
are  products  of  the  above  study.  These  programs  are  all  dirceted  at  the  deflection  problem  of 
antennas  although  the  forec  and  moment  resultants  are  also  available.  This  manual  is  open 
ended,  and  additions  will  be  made  as  new  developments  are  completed. 


Accepted  for  the  Air  Force 
Stanley  .1  Wisniewski 
Lt  Colonel,  USAF 
Chief,  Lincoln  Laboratory  Office 
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II.  FORMULATION  OF  PROBLEM 


A.  GEOMETRY  OF  SHELL 

The  middle  surface  of  the  antenna  is  a  paraboloid  of  revolution  whose  geometry  is  described 
by  the  following  figures  and  formulas.  A  point  o  on  the  middle  surface  is  located  by  rectangular 
Cartesian  coordinates  y^,  y^,  y^  or  by  circular  cylindrical  coordinates  r,  0,  y^  (Fig.  1) 


y^  =  r  cos  9 
y^  =  r  sin  O 
2 

y3  =  4f 

where  f  is  the  focal  length  of  the  parabola  (Fig.  2).  Let 


(II-D 

(H-2) 

(H-3) 


(II-4) 


Then  the  slope  of  the  parabola  (and  therefore  the  slope  of  a  meridian  of  the  paraboloidal  surface) 
is 


An  element  of  arc  length  along  a  meridian  (see  Fig.  3)  is 


(II-5) 


ds^  =  2f  \  1  +  y^  dy 

An  element  of  arc  length  along  a  latitude  (see  Fig.  3)  is 
ds^  =  2fyd0 


(II-6) 


(H-7) 


An  element  of  surface  area  on  the  middle  surface  is 


dA  =  4f 1  +  y^  dydG 


(II-8) 


A  point  p  in  the  shell  is  located  by  orthogonal  coordinates  r,  0,  £  (see  Fig.  3).  [£  is  perpendicu¬ 

lar  to  the  middle  surface  (see  Fig.  3)  and  is  positive  inward,  h  is  the  thickness  of  shell.] 

Let  ed  be  of  unit  length,  then  the  directional  cosines  of  £,  that  is,  the  orientation  of  £  (see 
Fig.  4),  are  given  by 


—  y  cos  9 
ec  -  1 


VI 


eb 


+  y 

y  sin  0 


+  y 


(Note  it  is  in  the  negative  y^  direction.) 


(Note  it  is  in  the  negative  y^  direction.) 


(II-9) 


(11-10) 


and 


ad  =■  1 


V7 


+  y 


dl-ii) 
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Fig.l.  Paraboloid  of  revolution. 


Fig.  2.  Focal  length  of  parabola. 
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The  principal  radius  of  curvature  along  a  meridian  (see  Fig.  3)  is 

Rr  =  2f(l  +  y2)3/2  (11-12) 

and  the  principal  radius  of  curvature  along  a  latitude  (see  Fig.  3)  is 

R0  =  2f  -/l  +  y2  (11-13) 

The  boundaries  or  edges  of  shell  are  described  by  (see  Fig.  5) 

r  =  r^  =  the  radius  of  the  inner  edge  (11-14) 

r  =  r^  =  the  radius  of  the  outer  edge  .  (II- 1  5 ) 

The  corresponding  value  of  y  at  the  edges  will  be  denoted  by  y^  and  y^,  respectively.  The 

directional  cosines  of  a  unit  vector  tangent  to  a  meridian  (see  Fig.  6)  are 


y^  component 

cos  0 

(11-16) 

\/l  +  y2 

y^  component 

sin  0 

(11-17) 

\l  1  +  y^ 

y^  component 

y 

(11-18) 

./.  ,  2 

The  directional  cosines  of  a  unit  vector  tangent  to  a  latitude  (see  Fig.  6)  are 


y^  component  —  sin0 

(11-19) 

y^  component  cos  0 

(11-20) 

y^  component  0 

(11-21) 

Finally,  we  define  a  pointing  angle  ip  to  be  the  angle  between  the  axis  of  revolution  of  the  parab¬ 
oloidal  surface  and  the  direction  of  gravity  (see  Fig.  7). 


B.  STRUCTURAL  PARAMETERS  OF  SHELL 


The  parameters  which  describe  the  structural  behavior  of  the  shell  are  the  middle -surface 
displacements,  strains,  rotations,  curvatures,  stresses,  force  resultants,  transverse  shear 
resultants,  and  moment  resultants.  These  will  now  be  defined. 

Let 


u 

v 


w 


displacement  of  point 
(see  Fig.  8) 

displacement  of  point 
(see  Fig.  8) 

displacement  of  point 
Fig.  8). 


o  along  tangent  to  meridian  at  point  o 
o  along  tangent  to  latitude  at  point  o 
o  along  normal  at  point  o  (see 


(11-22) 

(H-23) 

(11-24) 
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1  WI-1919  | 


JV-  FI-1930  | 


Fig.  6.  Base  vector  and  midsurface  normal. 


Fig.  7.  Pointing  angle. 


LATITUDINAL  CROSS  SECTION 


Fig.  8.  Deformation  of  shell  element. 
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Then  the  middle  surface  strains  and  rotations  are  given  by  the  formulas  listed  below: 
strain  along  meridian 

1  8u  w 


€  - 
r 


2f 


dr 


+  y 


2f(l  +  y2)3/2 


strain  along  latitude 
1  8v 

ee  zfy  ae 

shearing  strain 


2fy  ^77 


2f  s/l  +  y2 


re 


2f  nJ  1  +  y2  2fy  1  +  y‘ 

rotation  of  op  in  meridional  direction  (see  Fig.  8) 
u  1 


8v 

ay 


Ji 


1  8u 
2fy  8  0 


CO  =  — 

r 


2f(l  +  y2)3//2 


_  8w 

-z  8? 


2f  V  1  +  y" 

rotation  of  op  in  latitude  direction  (see  Fig.  8) 


Co  =  — 
© 


2f  si  1 


1  8  w 

2fy  8  0 


+  y 


change  in  curvature  along  meridian 

8w 


2f  J~l~  +  y 

change  in  curvature  along  latitude 


ay 


0 


2fy 


V7 


+  y 


twist 


cr0 


2fy  ae 


dco. 


2f 


n/1 


ay 


+  r 


2fy 


V7 


-  dco 
\ _ r 

2fy  80 


+  y 


(H-25) 


(11-26) 


(11-27) 


(11-28) 


(11-29) 


(11-30) 


(H-31) 


(11-32) 


The  stresses  at  point  p  in  the  shell  are  defined  by  means  of  Fig.  9.  It  is  more  convenient 
to  define  a  system  of  force  resultants  and  moment  resultants  acting  at  the  middle  surface  of  the 
shell  which  are  equipollent  to  the  stresses  integrated  over  the  thickness  (Figs.  10  and  11).  Let 


N  =  force  resultant  in  meridional  direction  = 
r 


rh/2 
\  <7 

J-h/2  r 


N._  =  force  resultant  in  latitude  direction  = 
0 


£ 


h/  2 


h/2 


CTe  d£ 


(11-33) 


(H-34) 
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Fig.  10.  Stress  resultants. 

h 


17-71-  1934  | 


]3-7l-  1935] 

TANGENT  TO  LATITUDE  AT  o 

L  » 


M 


r9 


\ 

TANGENT  TO  MERIDIAN  AT  o 


Fig.  11.  Moment  resultants .  (The 
double-headed  arrows  represent 
moments.  Use  right-hand  rule.) 
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=  NQr  =  shear  force  resultant  = 

p  h/2 

3 -h/2  <rr0df 

(H-35) 

Qr  =  transverse  shear  resultant  on 

r  face  =  ^ 

-.h/2 

'-h/2 

Crfd£ 

(11-36) 

=  transverse  shear  resultant  on 

0  face  =  ^ 

ph/2 

-h/2 

a0^df 

(11-37) 

=  bending  moment  about  tangent  to  latitude  at  o 


rh/2 

"  J-h/2 

My  =  bending  moment  about  tangent  to  meridian  at  o 

ph/2 


3. 


h/2 


(11-38) 


(11-39) 


M 


r0  =  ^0r  =  twining  moment  =  j 


h/Z 


h/2 


crrefdf 


(11-40) 


Then  an  element  of  volume  bounded  by  surfaces  r,  r  +  dr,  0,  0  +  d0,  and  ±h/Z  {see  Fig.  3) 
satisfies  the  following  equations  of  force  and  moment  equilibrium. 


9Nr  / - 2  9Nr0 

y  -57  Wl  +Y  W-  +Nr-Ne- 


- —  Q„  +  2fy  J77~yZ  pr  =  0 


1  +  y 


2  r 


3N 


r0 


t  2  r~  z 

+  V  1  +  y  -gy-  +  2Nrg  -  YQq  +  2f>  Vl  +  ?  P  ft  =  0 


3y 


9Q 


3Qf 


Y  *8^  +  */T+  y2  ^  + 


N 


90  '  A  2  r  r  0  ^r 

1  +  y 


+  yNft  +  Q  +  2fy  \  1  +  y2  p  =  0 


9M  i - j  9M 

+  VI  +  Y  -gf2  +  Mr  -  M0  -  2fy  V 1  +  y“  Qr  =  0 


J77  2 


9M 


re 


ay 


+  -Vl  +  y2  +  2M_/->  -  2fy  Ji  +  2 


ae 


+  2Mr0  -  2fy  N/l  +  y  Q0  =  0 


(11-41) 

(11-42) 

(11-43) 

(11-44) 

(11-45) 


where  and  are  the  components  of  the  applied  surface  loads  along  the  tangent  to  meridian 
and  to  a  latitude,  respectively;  p^  is  the  component  along  the  normal  to  the  middle  surface.  For 
a  shell  subjected  to  gravity  load  with  an  arbitrary  pointing  angle  ip, 


p  = —  H  fsin  0  sirup  —  y  cos  ip]  (11-46) 

r  4777 

p 0  =  P  cos  0  sinip  (11-47) 

and 

p  =— - £ - fy  sin  0  sin  ip  —  cos  ip]  (11-48) 

4777 
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where  p  is  the  surface  weight  density  of  the  shell  (lb/in ?“)  and  is  related  to  the  volume  weight 
density  p^  by 

rh/2 

P  =  \  ME  •  (11-49) 

J-h/2  ° 

For  a  homogeneous  shell,  we  have  p  =  p^h. 

For  an  isotropic  shell,  the  stress  resultants  and  the  moment  resultants  are  related  to  the 
shell  middle-surface  strains  and  curvature  changes  as  follows: 


N  =  C(e  +  v  <rQ-  AT) 
r  '  r  m  0 

(11-50) 

N =  C(<r~  +  v  e  -  AT) 

0  '  0  m  r 

(11-51) 

N  ~  =  C(1  —  v  )  e 
rO  m  r0 

(11-52) 

M  =  D(k  +  -  AT) 

r  r  b  0 

(11-53) 

M0  =  D<*e  + Vr"^ 

(11-54) 

Mre  =  D<!  -  V  Kr0  • 

(11-55) 

The  extensional  stiffness  C,  the  bending  stiffness  D,  the  effective  stretching  and  bending 
Poisson's  ratios  and  v ^  are  given  in  terms  of  the  Young's  modulus  E  and  Poisson's  ratio 
v  by 

rh/2  e 

C  =  \  —  7  df  (11-56) 

J-h/2  1  -  v 


i  rh/2  _jn 

C  *  -h/2  1  -  v 


Z 


rh/2 

°=  J_ 


E  ,  £2d£ 

h/2  1  -  i/ 


,  ^h/2  . 

1  (  '  v\ 


"b  =  D  J_ 


h/2  1  —  v 


rd? 


(II-57) 


(11-58) 


(11-59) 


The  membrane  and  bending  thermal  strains  AT  and  AT,  are  given  in  terms  of  E,  v,  the  change 
in  temperature  AT,  and  the  coefficient  of  thermal  expansion  of  the  shell  a  by 


AT 


AT=  ^ 


1 

ph/2 

a  EAT 

C  , 

'-h/2 

1  —  V 

1 

ph/2 

a  EAT 

D  , 

'-h/2 

1  -  V 

dr 


rdr 


(11-60) 


(11-61) 


If  the  shell  is  homogeneous  in  the  thickness  direction,  then 


v  -  v ,  =  v 
m  b 


(11-62) 
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c  = 


Eh 

l-,2 


D  = 


Eh 


12(1  -  O 


AT  = 


g(l  +  v)  f 


h/2 


3. 


^  _  12(1  +  v)  a  C 


h/2 
h/2 


ATd£ 


i 


h/2 


AT£d£ 


(11-63 


(11-64 


(11-65) 


(11-66) 


We  shall  have  occasion  in  the  subsequent  development  to  refer  to  a  parameter  defined  by 

2 

2  r2 
= - - - 


4hz 


(11-67) 


AD 


where  H  is  the  representative  magnitude  of  the  principal  radii  of  curvature.  For  a  paraboloida 
shell  of  revolution,  we  have  H  =  2f.  If,  in  addition,  the  shell  is  isotropic  and  homogeneous,  we 
have 


\  = - —  \l  12(1  -  vZ)  =  O 

\r2fh 


(i) 

'■mJ 


(11-68) 


Equations  (11-25)  through  (11-32),  (11-41)  through  (11-45),  and  (11-50)  through  (11-55)  form  a 
set  of  nineteen  equations  for  the  nineteen  structural  parameters  (N^,  NQ,  Q^,  Q^,  M^, 

Mr0j  c  ,  e  w  ,  c Oq,  *  ,  kq,  *  u,  v,  and  w).  They  are  the  shell  equations  for  a  parab¬ 
oloidal  shell  of  revolution. 

The  components  of  stress  at  a  point  o  are  related  to  the  stress  and  moment  resultants  by 
the  following  formulas. 


N 


12M 


<7  = 

r 

„r- 

r  £ 

3  b 

h 

(11-69) 

°6  = 

12Me  , 

h3 

(11-70) 

II 

© 

u 

b 

12Mr0  , 

3  b 

hi 

(11-71) 

(J  y  ~ 

r£ 

3Q  , 

—  1 

2h  V 

i_ 

h2/ 

(11-72) 

ffef  = 

3%  ( 
2h  V 

4_id\ 

h2) 

(H-73) 

C.  EDGES  OF  SHELL 

Along  a  y  =  constant  edge,  we  may  prescribe  an  appropriate  combination  of: 
(1)  Any  one  of  the  four  quantities 
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(a)  Normal  displacement  w 

(b)  Axial  displacement  uv 

=  w  cos  <p  +  u  sin  <p 

(c)  Effective  transverse  shear  +  (l/2fy)  (aM^/80) 

(d)  Effective  axial  resultant  V 

/  1  9Mr©\ 

V  lQr  +  2f7  ~a0~ )  COS  v  +  Nr  sin  V  * 

(2)  Any  one  of  the  four  quantities 

(a)  Meridional  displacement  u 

(b)  Radial  displacement  u^ 

uh  ~  —  w  s^n  +  u  cos  <P 

(c)  Meridional  stress  resultant 

(d)  Effective  radial  resultant  H 

/  i  aMre\ 

H  =  ~  lQr  +  2f7  ~dO~)  sin  *  +  Nr  cos  * 

(3)  Either  the  circumferential  displacement  v  or  the  effective  shear 
resultant  Nr0  +  (M^/Rq),  and 

(4)  Either  the  moment  M  or  the  rotation  u>  ,  where 

r  r 


(11-74) 


(H-75) 


(11-76) 


(11-77) 


cos  (p  = -  (11-78) 

vl  +  y2 

sin  <p  =  —  y  .  (11-79) 

-s/ 1  +  y2 


A  combination  of  prescribed  conditions  is  appropriate  if  the  structure  is  in  global  static  equilib¬ 
rium  under  these  conditions. 

The  standard  idealized  edges  supports  arc  merely  some  special  combinations  of  the  above 
general  set.  For  a  simply  supported  edge,  wc  have 

w  =  u  =  v  =  Mr  =  0  (11-80) 


For  a  clamped  edge,  we  have 


w  =  cj  =  u  -  v  0 
r 


(11-81) 


For  a  free  edge,  we  have 


Mre  l  aMr0 

=  M  =  Q  +  4-  0 


rG  R 


G 


cr  2fy  80 


(11-82) 


If  the  shell  is  closed  at  the  apex,  we  require  that  the  stresses  and  displacements  be  finite  at 
y  =  0. 
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D.  SOLUTION  TO  SHELL  EQUATIONS  FOR  GRAVITY  LOAD 

For  a  shell  under  gravity  load  and  without  edge  loads  and  temperature  gradients,  the  solution 

4 

to  our  shell  equations  must  take  the  following  form. 

<Nr'  Ne*  Qr’  Mr’  M0’  u'  w’  “r^  =  *Nr’  N0’  ^r’  Mr'  M0'  u&'  wS'  “r'  cos^ 

+  (Na,  Na,  Qra,  M^,  Mg,  ua,  wa,  wa)  sinG  simp  (11-83) 

(Nr0,  Mr0,  Q0,  v)  =  (Na0,  Ma0,  Qa,  va)  cos  G  simp  .  (11-84) 

In  the  subsequent  development,  we  will  not  be  concerned  with  the  other  seven  structural  quanti¬ 
ties  which  can  be  obtained  by  appropriate  combinations  of  those  previously  given.  It  should  be 
said  at  this  point  that,  although  the  form  of  the  solution  to  our  shell  equations  has  been  obtained 
with  relative  ease,  the  exact  solution  to  these  same  equations  does  not  seem  likely.  The  com¬ 
puter  programs  presented  herein  represent  several  approaches  to  an  approximate  solution. 

These  different  approaches  will  be  described  under  the  appropriate  programs  and  more  thoroughly 
discussed  in  Refs.  1-5.  In  the  remaining  portion  of  this  section,  we  shall  discuss  some  general 
features  of  the  output  of  these  programs  in  connection  with  the  shell  behavior  from  the  designer's 
point  of  view. 

The  superscripted  quantities  appearing  in  Eqs.  (11-83)  and  (11-84)  are  independent  of  ip  and 
0.  Their  dependence  on  y  is  generally  complicated.  It  is  clear  from  the  same  equations  that 
these  are  the  key  quantities  to  our  problem.  Once  they  are  determined,  the  physical  quantities 
appearing  on  the  left-hand  sides  of  Eqs.  (11-83)  and  (11-84)  can  be  obtained  for  any  value  of  ip  and 
0  by  straightforward  calculations.  The  computer  programs  discussed  herein  are  designed 
mainly  to  calculate  these  superscripted  quantities  for  a  given  set  of  values  of  y  which  will  be 
referred  to  as  the  y-tabular  points.  However,  since  ip  is  really  a  load  parameter,  its  effect 
on  the  stresses  and  deformations  of  the  shell  is  generally  of  prime  interest  to  the  designers  of 
larger  antennas.  Therefore,  the  programs  for  gravity  load  give  generally  the  following  quantities 
as  its  output. 

,XT  o  — , T  o  o  o.  iXT  s  s  s  s.  , 

(Nr  ,  N0,  .  .  . ,  w  ,  u»r )  =  (Nr  ,  N0,  .  .  . ,  w  ,  o»r )  cos  ip 

+  (Na,  N0,  .  .  .  ,  wa,  u>a)  sin  ip  (11-85) 

(Nr°0,  M°0,  Q°,  v°)  =  (Na0 . va)  sini p  .  (11-86) 

The  various  physical  quantities  themselves  will  also  be  calculated  upon  request.  To  this  end, 
we  must  prescribe,  in  addition  to  the  y-tabular  points,  a  set  of  values  for  0  which  will  be 
referred  to  as  the  0 -tabular  points.  Together  with  the  former,  they  form  a  net  of  points  each 
describing  a  particular  point  on  the  middle  surface  of  the  shell.  The  programs  calculate  the 
stress  and  displacement  of  the  shell  for  this  net  of  tabular  points.  In  all  cases,  the  unit  for 
length  is  inch  and  the  unit  for  force  is  pound.  While  the  shell  is  closed  in  the  circumferential 
direction  so  that  the  range  of  0  is  (0,  27r),  it  is  clearly  sufficient  to  calculate  for  any  particular 
latitude  (i.e.,  for  any  fixed  value  of  y)  the  values  of  the  physical  quantities  for  o  <  0  <  7r.  This 
rather  trivial  observation  may  save  a  great  deal  of  computing  time  if  a  large  number  of  runs 
is  desired. 


13 


It  is  often  desirable  to  present  numerical  results  in  a  form  which  is  independent  of  the 
various  geometrical  and  material  parameters  associated  with  the  shell.  In  general,  this  is  not 
possible  in  shell  analysis.  Thus,  we  can  only  hope  to  find  normalizing  factors  which  yield  di¬ 
mensionless  quantities  of  the  same  order  of  magnitude  for  shells  with  different  geometry  and 
different  materials.  To  this  end,  we  let 


M*,  M* 


q:. 


■*,  w*)  = 

C(1  -  v 

4f2p 

NrO>  - 

Wp  (IV 

Mro\ 

(mk  )2 

°  i 

mko/ 

2  ' 
4fp 

Qo\ 

mko/ 

mk 

Zip  <«r 

A 

1  “ 

C(1  -  v1 

(u,  V,  w) 


V  0’  rG 


r  2fp(mko)  r 


,  4  4 rc 

k  =  — pr- 

o  D 


4  2x 

m  =  (1  —  v  ) 


(11-87) 

(11-88) 

(11-89) 

(11-90) 

(11-91) 

(11-92) 

(H-93) 


Although  the  starred  quantities  are  not  invariant  with  respect  to  the  various  geometrical  and 
material  properties  of  the  shell  (for  instanee,  they  obviously  depend  on  Poisson's  ratio  v),  their 
magnitude  does  not  vary  appreciably  even  with  appreciable  changes  in  the  relevant  geometrical 
and  material  parameters.  Both  the  starred  (normalized)  quantities  and  the  unstarred  (un¬ 
normalized)  quantities  may  be  obtained  from  our  eomputer  programs. 

It  ean  be  seen  from  an  examination  of  the  solution  that  a  portion  of  the  displacements  cor¬ 
responds  to  a  rigid  body  translation  and/or  rigid  body  rotation  about  the  y1  axis.  For  example, 

in  Group  Report  71G-1-II  (Ref.  2),  the  constant  C?  in  Eq.  (6.4.21)  for  w  (membrane  analysis)  is 

L  3 

seen  to  represent  a  translation  of  the  shell  parallel  to  the  y  axis.  Sinee  this  is  a  rigid  body 
motion,  C ^  does  not  appear  in  any  of  the  equations  for  the  force  resultants.  Similarly,  in  the 
asymptotic  analysis,  the  constant  C7  [Eqs.  (8.9.1)  and  (8.9.2),  Group  Report  71G-1,  Part  IV 

L  3 

(Ref.  4)]  denotes  a  rigid  body  translation  parallel  to  the  y  axis  for  symmetric  loads.  For  anti¬ 
symmetric  loads,  the  constants  and  C4  (Eqs.  (8.10.1)  to  (8.10.3),  Group  Report  71G-1,  Part 

3  4  1 

IV  (Ref.  4)]  represent  a  rigid  body  rotation  about  the  y  axis  and  a  rigid  body  translation  parallel 
2 

to  the  y  axis,  respectively.  The  rigid  body  displacements  do  not  affeet  the  shape  of  the  parab¬ 
oloidal  surface  and  henee  can  be  interpreted  as  a  change  in  the  position  of  the  focal  point.  We 

/  .  A  A  A 

will  denote  the  defleetions  due  to  rigid  body  translations  and/or  rotation  by  u,  v,  and  w.  Then, 
the  actual  flexibility,  i.e.,  distortions,  will  be  given  by 


A 

uri  u  —  u 

flex 

(11-94) 

v„  A 

flex  =  v  —  v 

(11-95) 

A 

w„  w  —  w 

flex 

(11-96) 
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From  the  previous  discussion,  it  is  clear  that  these  quantities  should  be  of  considerable  interest 
to  the  designers  of  large  antennas.  Our  computer  program  will  display  these  as  part  of  the 
output. 

Thus  far,  it  has  been  assumed  that  the  shell  has  its  prescribed  shape  before  the  introduction 
of  the  applied  loads.  Once  loaded,  the  distortion  of  the  shell  is  taken  with  respect  to  this  pre¬ 
scribed  (paraboloidal)  shape.  The  antenna  designer  may  also  be  interested  in  the  distortions  of 
the  shell  relative  to  its  shape  when  the  axis  of  the  shell  coincides  with  gravity  vector.  This  is 
usually  the  attitude  of  the  shell  during  erection  and  the  designer  may  elect  to  construct  the  shell 
to  a  prescribed  shape  while  it  is  in  this  attitude.  Thus,  we  are  often  interested  in 


w  =  wS(cos  ip  —  1)  +  wa  sin  0  sin  ^ 

(H-97) 

u  =  uS(cos  ip  —  1)  +  ua  sin  0  sin  ip 

(11-98) 

v  =  va  cos  0  sin  ip 

(11-99) 

where  the  tilde  quantities  are  the  components  of  displacement  relative  to  the  face-up  (i.e.,  ip  =  0) 
position.  The  corresponding  quantities  uflex,  vflex,  wflex,  u°,  v°.  w°.  u°flex,  v°flex,  and  w°flex 
are  defined  in  the  obvious  way. 
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III.  MEMBRANE  ANALYSIS  OF  GRAVITY  LOAD 


A.  MEMBRANE  SOLUTION 

Under  conducive  circumstances,  the  interior  of  the  shell  behaves  like  a  membrane.  Away 
from  the  edges  of  the  structure,  a  good  approximation  of  the  behavior  of  the  shell  can  be  obtained 
by  the  so-called  membrane  (or  momentless)  theory  of  shell.  Such  a  theory  assumes  that  the 
shell  has  no  bending  stiffness,  and  quantities  associated  with  the  bending  action  of  the  shell  vanish 
identically,  i.e., 

Mr=Me=Mre  =  Qr=Qe  =  0  • 


Consequently,  Eqs.  (11-41)  to  (11-43)  become 
8N 


8y 


J 1  + 


— r  3N  Q 

y2  -  ^-0  +  N 
7  80  r 


N0  +  2fy 


77 


+  Y  Pr  =  0 


(IH-2) 


9N 


rO 


r— 2  aN0  + 
n/i  +  y  to  + 


3y 


2Nr0  +  2fy 


A 


2 

+  Y  P, 


=  0 


(III  —  3 ) 


iv  / - j 

+  N„  +  2f  V 1  +  y p  =  0  (III -4) 

.  2  0  *n 

1  +  y 

and  Eqs.  (11-44)  and  (11-45)  are  identically  satisfied. 

Equations  (111-2)  to  (I II — 4 )  contain  only  three  unknowns  and  can  therefore  be  solved  for  the 
stress  resultants  N  ,  N0,  and  N^0.  In  other  words,  the  membrane  problem  is  statically  deter¬ 
minate.  The  solution  for  N  ,  N~,  and  N  „  can  then  be  inserted  in  the  left-hand  side  of  Eqs. 

r  0  r0  1 

(11-25)  to  (11-27)  by  way  of  Eqs.  (11-50)  to  (11-52)  with  the  temperature  terms  omitted.  The  three 
strain-displacement  relations  also  contain  only  three  unknowns  and  can  therefore  be  solved  for 
u,  v,  and  w. 

Note  that  associated  with  the  momentless  assumption  is  a  reduction  of  the  order  of  our  sys¬ 
tem  of  differential  (shell)  equations  from  eight  to  four.  Thus,  only  half  the  conditions  at  each 
edge  (Sec.  Il-C)  can  be  satisfied.  For  a  free  edge,  the  last  two  conditions  [Eq.  (11-82)]  are  satis¬ 
fied  identically;  therefore,  we  are  left  with  the  necessary  two  boundary  conditions.  In  the  case 
of  a  supported  edge,  it  has  been  shown  that  the  appropriate  boundary  conditions  are 

u  =  0  (III-  5) 


and 


v  -  0  .  ( II I -  6) 

Since  the  shell  can  not  be  acted  upon  by  transverse  forces  or  bending  moments,  it  seems  rea¬ 
sonable  that  it  may  not  be  constrained  in  the  normal  direction. 

The  exact  analytical  solution  for  the  membrane  stress  resultants  as  well  as  the  middle- 
surface  displacement  components  in  accordance  with  the  momentless  theory  have  been  obtained 
and  tabulated  in  Ref.  2.  It  can  be  seen  from  these  analytical  solutions  that  the  corresponding 
normalizing  factors  given  in  Sec.  II-D  are  nearly  the  true  scale  factors  for  these  quantities.  For 
a  fixed  value  of  r,  one  set  of  normalized  results  generated  from  our  computer  program  is  valid 
for  the  entire  class  of  geometrically  similar  shells,  i.e.,  shells  with  the  same  y^  and  y^. 

The  analytic  expressions  for  the  various  physical  quantities  determined  by  the  momentless 
theory  of  shell  are  formed  by  complicated  combinations  of  elementary  functions  To  evaluate 
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them  for  a  large  number  of  tabular  points  even  by  a  desk  calculator  is  a  formidable  task.  A 
computer  program  is  written  to  make  this  information  easily  accessible  to  the  designers  of  large 
antennas  and  to  others  who  might  be  interested.  The  program  generates  the  tabular  points,  com¬ 
putes  the  constants  of  integration,  and  calculates  the  stress  resultants  and  displacements  of  the 
middle  surface. 

B.  DIGITAL  COMPUTER  PROGRAM 

The  general  scheme  of  the  computer  program  for  a  membrane  analysis  is  outlined  in  Fig.  12, 
the  master  flow  chart.  At  the  end  of  a  complete  run,  the  program  returns  to  step  (1).  When  it 


Fig.  12.  Master  flow  chart. 

fails  to  locate  additional  input,  the  program  exits  automatically.  This  is  a  special  feature  of 
the  Lincoln  I  aboratory  Express  Runs.  Slight  modifications  of  the  ending  will  be  necessary  if 
the  program  is  to  be  run  differently  or  if  some  indication  of  a  successful  run  is  desired.  The 
same  Express  Runs  also  require  that  the  input  be  prestored  on  machine  tape  A2  and  that  the  out¬ 
put  be  written  on  A3.  In  complying  with  this  restriction,  the  program  reads  all  its  input  from 
A2  and  writes  all  its  output  on  A3.  For  the  purpose  of  the  Express  Runs,  the  object  deck  must 
be  preceded  by  two  cards  each  with  an  asterisk  in  column  one.  The  first  of  these  is  an  identifi¬ 
cation  card  while  the  second  contains  the  work  XEQ  occupying  columns  seven,  eight,  and  nine 
(Sample  Card  1).  Following  the  object  deck  is  another  card  with  an  asterisk  in  column  one  and 
the  word  DATA  occupying  columns  seven,  eight,  nine,  and  ten  (Sample  Card  2).  Then  comes 
the  input  to  the  program  to  be  prestored  on  A2.  The  program  listing  will  be  presented  in 
Table  I. 
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Sample  Card  1 
*  XEQ 


XEQ 


0  0  0  3  0  0  0  0  0  1  0  3  0  0  0  il  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  o  0  0  3  0  0  0  0  0  0  C  0  0  0  0  0  0  o  0  3  0  0  0  0  0  0  0  o  0  0  r,  0  0  o  0  0  0  n  0  0 

’  I  .  :i  *4  tt  'i  1»  tl  It  ??  7174  75  7f  .'f  4*5  31  Jl?7tt‘>4*iX  V  <1  31  40  41  47  4344  -M,  «J  Si  il  .«  V  x  J?  ■  ,<■  --  '  fUJft  F’.  (il  '  .  f  K  <t  •  n  vu 

1111111111111111111111111111111111111111111111111111)111111111111111.11111111111 

?  2  2  2  2  2  2  ?  ?  I  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  ?  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  7  7  7  7  7  2  2  2  2  ?  2  2  2  2  2  ? 

]  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  J  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  1  3  3  3  3  3  3  3  3  3  3  3  3  3  2  2 

4  4  4  4  4  4  4  4  1  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  M  4  4  M  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4 

^55  5  5  5  5  5  5  5  $  5  5  5  5  5  5  5  5  5  5  5  S 5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5 5 5  5  5  5  5  5  5  5  5  5  5  5  5  S 5  5  5  5  5  5 5  5  r  5  5  5  5  5  5  5  5  5  5  5  5 

6  6  P  6  fi  6  6  6  B  B  8  b  6  6  6  6  6  b  6  6  6  6  6  6  6  6  6  C  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  f  6  b  6  6  M  6  6  6  b  6  6  6  6  6  6  C  C  6  f>  I  (5  fi  6  E .6  6  6  P  C  6  C  f  6  6  6 

mm  n  mm  mu  iminiiiiiiiiiiiimiiiiin  in  mnimn  n  inni  in  mm 

8  ?  4  S  i  4  S  S  8  3  8  8  8  !  8  8  8  i  8  8  8  8  8  8  8  MS  1  8  8  8  8  8  8  8  8  8  8  8  8  8  M  8  8  8  8  8  8  8  8  8  8  8  8  8  8  M  8  M  8  8  88  M  g  !  8  ft  8  8  f  6  4 

4  n  9  3  3  9  9  9  9  2  *  3  iS  9  $  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  S  9  M  9  9  9  9  9  5  9  9  9  9  9  9  9  9  *  2  *  0  9  9  9  3  5  j  n  ;:••«?<?  9  2  9  5  5  1  9  9  S  «  S  C 

VC1  i»  i5  .*•  71  7?  *  <•  77  it  71  ’<3  l<  X’  M  Jin  *  J»  Jf  4H«I  c  «  «  «5  i*  «  M* !  *  «  «;• il'l  •  .  MO  f,  t*,.  ‘  ‘  ''.KO"  M  '  ’  * 


Sample  Card  2  DATA 
*  DATA 


onoooooo  OOOOODO&QOOOOOOOOOOOOnooOOOOOOOOOOOOOOOOOOOOODOOOOOOCOOODOOOOCCOOOOaOOO 

/  ’  *  i  t  7  i  i  '0  t  w  t?  *4  t$  s  ■■  i  »  ?<*  ?i  *t  -)  .v"!  ,•  *d  « i;  i: « (*  tt  ir  «  « v  rjj  « y.  *  •  s’  n  r  ’.tvis 

I  M  1  I  1  1  I  11111111111111111111111111111111111111111111111111111111111111111)1111 
?  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2 
33333333  33333333333333333333333333333333333333333333333333333333333333333333333 
4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  1  4  4  4  4  4  4  4  4  4  4  4  t  4  4  4  4  4  4  4  4  A  4  4  4  »  4  4  4  4  4  4 

3S5555555555555555555555555555555555555555555555555555!i555555555555555555555‘555 

$  6  6  fi  6  6  8  6  6  6  6  6  6  b  6  6  6  6  6  6  6  6  6  6  6  1  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  8  6  b  6  6  6  6  f  6  6  6  6  f  6  6  6  6  6  6  F  6  6  8  G  6  6  6  R  F  6  R  C  6 

nnimimmmmmmmnmnmmminmiiiiniiiiiinniiiiiiiiiin 

8  8  8  8  8  8  8  8  M  8  8  8  8  8  8  8  8  8  I  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  b  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  1  8  8  8  8  8  8  8  8  8  8  8  8  !  8 

S  9  9  3  *  9  9  9  9  9  3  9  9  9  9  9  9  9  9  9  9  9  9  9  9  5  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  3  9  9  9  S  5  9  9  9  9  9  9  9  9  9  9  9  5  9  5  9  9  9  S  9  S  1  9  9  9  9  c  c  f  9  9  9 

l  !  >  '  |  I  <7  II  i  7  *4  IS  If  17  H  W  70  .  77  J  7S  JS  77  ?t  73  1C  \  J?  Jt  M  35  *  V  7»  31 10  41  4?  11  M  4'  41  •’  4»  '*  ?f  I  i7 ' 1  V  '*tav  C  ~  U  f  •  ,«  ;  \r 

Hi'  *0«i 
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TABLE  I 

PROGRAM  LISTING 


*  SHEA-WAN  MEMBRANE  SHELL  WITH  GRAVITY  LOAD  9/3/64 

*  XEO 

C  WAN.  FRED  MEMBRANE  SHELL  WITH  GRAVITY  LOAD  9/3/64 

DIMENSION  GAMMA! 150) .THETA! 10) .ENR! 150,10) , ENT! 150,10) , ENRT (150,10 
1),U(150,10),V(150,10),W(150,10),EN1(150),EN2(150),U1(150),W1(150), 
2R(150),A(150)»B(150).C(150),G(150),P(150),Q(150).D(3,3),Y(3),Z(3)» 
3GAMMRI 150) ,S( 150) , FMT 1 ( 12) ,FMT2( 12) 

DIMENSION  UFLEXA! 150) ,VFLEXA( 150) .WFLEXA! 150) .UFLEXS! 150) .WFLEXS 
1(150) ,UFLEX (150,10) ,VFLEX( 150,10) ,WFLEX (150,10) 

COMMON  ENR,ENT»ENRT,V»W,U 

EQUIVALENCE  (GAMMA! 1 ) ,R ( 1 ) ) , ( GAMMR ( 1 ) ,S ( 1 ) ) 

C 

C - (1)  INPUT 

C 

10  READ  INPUT  TAPE  2 , 101 »M , N , KE , NORM ,KD , KP 
READ  INPUT  TAPE  2 , 105 ,R 1 , R2 , R3 , R4 
READ  INPUT  TAPE  2 , 105 , F , PR , PS  I 
X 1 = 1 • +PR 
X6=1.-PR 
X7=X1*X6 

IF ( KP-l ) 320,321,322 

320  READ  INPUT  TAPE  2 , 102  ,RHO , H , E 
RHO=RHO*H 

CC  =  E*H/X  7 
GO  TO  323 

321  READ  INPUT  TAPE  2,102.RHO,CC 
GO  TO  323 

322  READ  INPUT  TAPE  2 , 102 ,RHS ,RHC , T , H , E 
RH0  =  2.*T*RHS+ ( H-T  >*RHC 
CC=2.*T*E/X7 

323  I F ( N ) 7 , 7 , 8 

8  READ  INPUT  TAPE  2 , 105 ,THETA3 , THETA4 
READ  INPUT  TAPE  2 , 205 , FMT 1 , FMT 2 


C 

C 

C 


-  -  (2)  OUTPUT  I 


7  WRITE  OUTPUT  TAPE 
I F ( KP-l )67,68,69 

67  WRITE  OUTPUT  TAPE 
GO  TO  73 

68  WRITE  OUTPUT  TAPE 
GO  TO  73 

69  WRITE  OUTPUT  TAPE 

73  I F ( NORM  174,74,75 

74  WRITE  OUTPUT  TAPE 
GO  TO  302 

75  WRITE  OUTPUT  TAPE 
302  IFIKE-l 193,92,50 

93  WRITE  OUTPUT  TAPE 
GO  TO  91 

92  WRITE  OUTPUT  TAPE 
GO  TO  91 

50  WRITE  OUTPUT  TAPE 

91  IFIKP-l >76,77,78 

76  WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 


3,107 

3.160 

3.161 

3.162 

3.139 

3.140 

3.115 

3.116 
3,118 
3,112 

3,111  ,RHO,E,H,F ,PR 


GO  TO  79 

77  WRITE  OUTPUT  TAPE  3,163 
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TABLE  I  (Continued) 


WRITE 
GO  TO 

OUTPUT 

79 

tape 

3* 

111 ♦ CC  *  RHO  *  F  *  PR 

78 

WRITE 

OUTPUT 

TAPE 

3* 

1  64 

WRITE 

OUTPUT 

TAPE 

3* 

111  ,  E  *  H  *  T  *  RHS 

WRITE 

OUTPUT 

TAPE 

3* 

114 

WRITE 

OUTPUT 

TAPE 

3* 

165 

WR  I  TE 

OUTPUT 

tape 

3* 

1 1 1  *  RHC  *  RHO  *  F *PR 

79 

WR  I  TE 

OUTPUT 

TAPE 

3* 

114 

PI  =3.14159263 
RAD=57. 2957795 
PSIDEG=PSI*PI*RAD 
WRITE  OUTPUT  TAPE  3*113 

WRITE  OUTPUT  TAPE  3  *  1 08  * R 1  * R2 *R 3  * R4 *  PS  I  DEG 
WRITE  OUTPUT  TAPE  3*114 
I F ( N ) 1*1*9 

9  T3DEG=THETA3*P I*RAD 
T4DEG=THETA4*P I*RAD 
WRITE  OUTPUT  TAPE  3*117 
WRITE  OUTPUT  TAPE  3  *  1 1 9  * T 3DEG * T4DEG 
1  WRITE  OUTPUT  TAPE  3*114 
C 

c  -  -  (3)  clearing  storage  for  subscripted  variables 

c 

DO  51  1=1,150 

en i  (  n=o. 

EN2( I)=0. 

U1  (I) =0  . 

W1  (  I  ) =0. 

A( I ) =0. 

B(  I  )  =0. 

C(  I  ) =  0  • 

G( I ) =0. 

P (  I  )  =0. 

Q (  I  )  =0  • 

GAMMR (I  )  =C  • 

51  GAMMA ( I ) =0. 

DO  52  J=l*10 
THETA( J)=0. 

DO  52  1=1*150 
U(  I  * J) =0  . 

W(  I  *  J ) =  0  • 

V  (  I  *  J  )  =  0  . 

ENR (  I  * J  )  =0. 

EN  T (  I *J)=0. 

52  ENR  T (  I  *  J  )  =0 • 

C 

C  -  -  (4)  GENERATE  INDEPENDENT  VARIABLES 
C 

F2  =  2 • *F 

UR=. 00000005 

TWOP I  =  2 • *P  I 

R1=R1/F2 

R2=R2/F2 

R3=R3/F2 

R4=R4/F2 

M 1 =M— 1 

EM  1  =  M 1 

DELTR=R4-R3 
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TABLE  I  (Continued) 


DELTR=DELTR/EM1 
GAMMA { 1  ) =R  3 
DO  30  1=1, Ml 
11=1+1 
E  I  =  I 

30  GAMMA{ I  1 )=R3+DELTR*EI 
I F { NORM ) 630 ,630,632 

630  DO  631  1=1, M 

631  GAMMR { I ) =GAMMA ( I ) 

GO  TO  634 

632  DO  633  1=1, M 

633  GAMMR { I ) =GAMMA ( I ) *F2 

634  I F ( N ) 5,5,40 

40  THETA3=THETA3*PI 
THETA4=THETA4*PI 
N 1 =N-1 
AN  1  =N 1 

THETA5=THETA4-THETA3 
DELTT=THETA5/AN1 
THETA ( 1  )  =THET  A  3 
DO  2  1  =  1,  N1 
11=1+1 

2  THETA { 1 1  )  =THET A (  I ) +DELTT 
IF ( THET  A 4- T WOP  I  )  5*6*6 

6  THETA (N  )  =THETA3 
C 

C  -  -  (5)  CALCULATE  THE  CONSTANTS  OF  INTEGRATION 

C 

5  CS=F2*RHO 

CD=F2**2*RHO/ { X7*CC) 

PHI=PI*PSI 
CPH I =COS  F ( PH  I  ) 

SPHI=$INF(PHI ) 

I F ( PS  I ) 38*31 *32 
38  WRITE  OUTPUT  TAPE  3*143 
GO  TO  10 

32  CSS=CS*CPH I / 3 • 

CDS=CD*CPHI /3. 

CSA=2»*CS*SPHI 
CDA  =  CD*S  PH  I 

31  X2  =  • 5  * ( 3 • -PR ) 

X3  = • 2  5 

X4  =  “2 • * { 1 • -PR  ) 

X5  =  PR 

IF{KE“1 ) 95*96*98 
98  A 1 =- 1 • 

Cl=-1 • / 3 • 

C2=-l • / 1 5 • 

ARGl  =  1 •  +R2*  *2 
ARG1=S0RTF{ ARG1 ) 

X0=- ( 1 • +PR ) *A 1 

A3=“ARG1* { A1*R2+X0/R2 )-XO*R2*LOGF ( 1 ,+ARGl ) +{ XO-X 1 ) *R2*LOGF  <R2 ) -X2* 
1R2**3-X3*R2**5+X 1/R2 
A3=A3/R2 
GO  TO  94 
96  Al=l.+Rl**2 
A 1 =SORT  F { A  1 ) 

A  1  =-A 1**3 
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TABLE  I  (Continued) 


ARG1= 1 • +R2**2 
ARG1=SQRTF( ARG1 ) 

X0=— ( 1 . +  PR ) *A1 

A3=~ARGl*( Al*R2+XO/R2)-XO*R2*LOGF(l • +ARG 1 )+(XO“Xl ) *R2*L0GF ( R2 ) -X2* 
1R2**3-X3*R2**5+X 1/R2 
A3=A3/R2 
ARG1=1.+R1**2 
ARG1=SQRTF( ARG1 ) 

C 1 =-ARG 1 **3 /3 • 

C2  =- (  .5*C1*R1**2  +  ARG1**5/  15. ) 

94  ARG1=SQRTF( l.+R2**2) 

AR  G2  =  ARG 1 - 1 . 

ARG3=ARG2/R2 

H 1 = ARG 1  * ( R  2  *.5-Xl/R2  ) 

H2=X6*.5*R2  *LOGF ( ARG3) +X1*ARG1* ( 1 .-.5*  R2  **2)/R2  **3 

H3=R2**5-( 40.+29.*PR  ) *R 2 * *3  +  X 6* 2 . *R 2 * LOGF (R2)+X1*(4.-20.*R2**2 

1  ) /R2**3 
H3=H3/60. 

G 1 =“X 1  * ( ARG 1  + LOGF (ARG 3)  ) +ARGl**3/6. 

G2=.25*X6*(R2  **2*LOGF ( ARG3 ) + 1 . + ARG 1 ) - . 5*X 1  * ARG 1  * *3 /R 2  **2 

G  3  =  R2  **6-1 .5* ( 40.+29.*PR  )*R2  **4-3.*X6*R2  **2+LOGF(R2  )*( 

16.*X6*R2  **2-Xl*120. ) -X 1  *  12. /R2  **2 

G3=G3/360. 

TERM1=-H2+2.*X1*ARG1/R2**3 
TERM2  =H 1 +  X 1  * ARG 1 /R2 

TERM3  =  -H3  +  2.*X1*(  1.-4  .*R 2**2 ) * ARG 1**4/ (  15.*R2**3  ) 

C3= ( TERM3+C2*TERM1-TERM2*C1 ) /R2 
C4  =  -( »5*C3*R2**2+C2*G2+Cl*Gl+G3 ) 

GO  TO  5  A 

95  ARGll=SORTF(l.+Rl**2) 

ARG 1 2  =5QR  T  F ( 1 .+R2**2 ) 

ARG3= ( 1 .+ARG1  1 ) /R  1 
ARG4  = ( 1 . +ARG 12 ) / R  2 

T1=R1*ARG1 1-Xl* ( Rl*LOGF (ARG3 )+ARGll/Rl ) 

T2=R2*ARG12-X1* (R2*LOGF ( ARGA )+ARGl2/R2 ) 

T3=-X3*R1**5-X2*R1**3-(R1*L0GF(R1 )-l./Rl ) *X1 

TA  =  -X3*R2**5-X2*R2**3-( R2  *LOGF (R2)-l./R2)*Xl 

E 1 =R2*T  3-T A*R 1 

E2=T3*T2-TA*T1 

DEN0M=T1*R2-T2*R1 

A 1 =E 1 /DE  NOM 

A3  =  -E  2/ DENOM 

ARG2 1 = ARG1 1 -1 • 

AR  G22  = ARG 1 2 - 1 • 

ARG31=ARG21 /R1 
ARG32=ARG22/R2 
H11  =  ARG11*(R1*.5-X1/R1  ) 

H12  =  ARG12* ( R2* . 5-X  1/R2  ) 

H2  2  =  X  6* • 5 *R 2* LOGF ( ARG32 ) +X1*ARG12*  (  l.-.5*R2**2>/R2**3 
H2  1=X6*.5*R1*L0GF( ARG31 ) +X 1 *ARG1 1  * (  1  .  -  .  5* R  1  ** 2 ) /R 1 ** 3 
H31=Rl**5-<  AO. +29.*PR  )*R1**3  +  X6*2«*R l*LOGF( R 1)+X1*(4.-20.*R1** 

1 2 ) /R 1 **  3 

H32  =  R2**5-( A0.+29.*PR  )  *R2**3  +  X6*2» *R2*LOGF (R2)+X1*(4.-20.*R2** 

12 ) /R 2 **  3 
H31=H31/60. 

H32=H32 / 60 • 

G 1 1 =  — X 1  * ( ARG 11  +  LOGF ( ARG 31  )  ) + ARG 1 1 ** 3 /6  . 

G12  =  -X1*(ARG 12  +  LOGF ( ARG 3 2  )  ) + ARG 1 2** 3 /6 . 
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TABLE  I  (Continued) 


G21=.25*X6*(R1**2*L0GF(ARG31 )  +  l .+ARG11 ) -. 5*X1*ARG1 1**3 /R 1**2 
G22=.25*X6*(R2**2*L0GF( ARG32 )+l.+ ARG 12 )-.5*X 1*ARG 12**3 /R2**2 
G31=R1**6“1.5*(40.+29.*PR  ) *R 1**4- 3. *X6*R 1 **2+L0GF ( R 1 )*( 6.*X6*R 

11**2-X1 *120. ) -XI* 12. /Rl**2 

G32=R2**6-1 .5*( 40.+29.*PR  ) *R2 **4-3 . *X6*R2**2 +L0GF ( R2 ) * ( 6 . * X6*R 

1 2  **  2-X 1*120. )-Xl*12./R2**2 
G  3 1 =G  3 1 / 360 • 

G  32  =G  32 / 360 . 

D (  M  )  =G  1 1  -G  1 2 
D ( 1*2)=G21-G22 
D  (  1»3)  =  .5*(R1**2-R2**2) 

D  (  2  » 1  ) =  H  1 1+X1*ARG11/R1 
D (  2  *2  ) =H21-2.*X1*ARG11/R1**3 
D ( 2 . 3 ) =R 1 

D  (  3  ♦  1 ) =H12+X1*ARG12/R2 
D  (  3 . 2  ) =H22-2.*X1*ARG12/R2**3 
D ( 3 . 3 ) =  R2 
2(1) =  G3  2 -G3 1 

2(2) =-H31+2.*Xl*ARGl  1**4* ( 1 .-4.*R1**2 ) / ( 15 .*R 1**3) 

2(3) =-H32+2.*Xl*ARG12**4* ( 1 . -4 . *R2** 2 ) / ( 15.*R2**3) 

DET*1 • 

L  =  XS  I MEQF (3*3*1 *  D  *2 ♦ DET  *  Y  ) 

GO  TO  (70.71  .72  )  .L 
72  WRITE  OUTPUT  TAPE  3.103 
GO  TO  10 

71  WRITE  OUTPUT  TAPE  3.104 
GO  TO  10 
70  Cl-DUtl  ) 

C2*D( 2  1 1 ) 

C3  =  D( 3.1  ) 

C4  =  -  ( G3 1 +C 1 *G 1 1+C2*G21+C3*. 5*R 1 **  2 ) 

X0  =  -X 1  * A  1 
C 

C  -  -  (6)  GENERATE  THE  DEPENDENT  VARIABLES 

C 

54  I  F ( PS  I  J38.35.39 
39  DO  4  1=1 .M 

ARG3= 1 • +GAMMA ( I )**2 
ARG4=S0RTF( ARG3) 

ARG 1  =  ARG  3- 1 • 

ARG2=ARG1 /GAMMA ( I ) 

BEBA=( .5*C1*ARG4/R( I  )  )  +  ( C2*ARG4/R ( I  ) ** 3  )  +  ( ARG3** 3/ ( 1 5 . *R ( I  )**3)  ) 
BEBA=-BERA 

BEBC=(-.5*C1/R(  I  )  )  +  ( C2/R (  I  )  **3  )  +(ARG4**3* ( 1 • -4 . *R ( I  ) *  *  2  )  / ( 1 5 • *R (  I 
1  )  **3 )  ) 

BEBB  =  .5*R(  I  )+(  (C2/R(  I  )  **3 )  +  ( .5*C1/R ( I ) ) +ARG4* ARG3**2/ ( 1 5 . *R (  I  ) * *3 ) 
1  )  / ARG4 
A (  I  ) =CS A  *BEBA 
B ( I ) =CS A*BEBB 
C ( I ) =CS  A  *BFRC 
ARG 1 =SQR  TF ( 1 . +R (  I  )  **2  ) 

ARG2  =  ARG 1-1 . 

APG3  =  ARG2/R (  I  ) 

G1=-X1* (ARG1+LOGF (ARG3) ) +ARGl**3/6. 

G2  =  • 2  5*X6* ( R (  I  )**2*L0GF(ARG3)  +  1.+ARG1)-.5*X1*ARG1**3/R(  I  )  * *2 
G3=R( I )**6-1.5*(40.+29.*PR  )*R(  I ) **4-3 . *X6*R ( I )**2+L0GF(R( I ) )*( 

16. *X6*R (  I  )**2-X 1*120. ) -X 1 *  1 2 • /R (  I  ) **2 
G3=G3/360. 


23 


TABLE  I  (Continued) 


VFLEX A (  I  )=2.*CDA*(C1*G1+C2*G2+G3) 

G (  I  )=2.*(C4+.5*C3*R<  I  )**2+C2*G2+Cl*Gl+G3 )*CDA 
HI = ARG1 * ( R (  I  )*.C'-X1/R(  I  )  ) 

H2=X6*.5*R ( I )*LOGF ( ARG3) +X1*ARG1*( 1  .-.5*  R(I)**2)/R(I)**3 

H3=R ( I ) **5- (40.+29.*PR  )*R(I)**3+X6*2.*R(I) *LOGF( R( I ) )+Xl*( 4.-2 

10. *R ( I ) **2 ) /R ( I ) **3 
H3=H3/60. 

G4  =  -X1+(G1-R(  I  )*H1  1/ARG1 
G5=2.*X1*ARG1+R(  I  )**2*(G2-R(  I  ) *H2 ) 

G5=G5/( ARG1*R( I )**2) 

G6  =  ( G3-R ( I ) *H3 ) +2 . *X 1  * ( 1 . -4 . *R ( 1  ) **2 ) * ARG1 **4/ ( 1 5 . *R (  I ) **2 ) 
G6=G6/ARG1 

G7=(G4-ARG1*G1-.5*X1-.5*PR  *R (  I  ) **2 ) /R (  I ) 

G8=( G5-ARG1*G2-PR  -X 1 / R ( I ) **2  )  /R (  I  ) 

G9=-.5*R( I >*<  2.+RI I )**2 1/ARG1 
Gl0  =  -R(  I  )/ARGl 

G1 1=ARG1*( 15.*R  <  I )**4+2 . *ARG 1 **4+2 . *PR  *ARG1**6 )/(30.*R(I)**2) 

G1  1  =  <G6-ARG1*G3-G11 )/R( I) 

P(  I  )  =2.*CDA*(C1*G7  +C2*G8  +C3*G9  +C4*G10+G11) 

0( I ) =2.*CDA*( Cl*G4+C2*G5  +  G6+<  C4-.5*C3*R ( I ) ** 2 ) / ARG1 ) 

WFLEX A (  I  )=2.*CDA*(C1*G7+C2*G8+G11 ) 

UFLEX A (  I  )=2.*CDA*(C1*G4+C2*G5+G6) 

4  CONTINUE 

UB  1  = ABSF (0(1)) 

UR2=ABSF (0(H) ) 

IF (UB1-UR)47  *47 . 53 
47  0(1) =0. 

53  IF (UB2-UB) 58.58.35 

58  0 ( M ) =0 . 

3  5  DO  11  I =  1 .  M 

ARG 1 = 1 • +GAMMA ( I )**2 
ARG1=S0RTF( ARG1 ) 

ARG2=1. /GAMMA!  I  ) 

ARG2=ARG2**2 

EN 1 ( I ) =  CS  *  ARG 1*ARG2* ( ARG 1 **3+Al )/3. 

EN2 ( I ) =CS  *( 3.-ARG2*! ARG1**3+A1 ) /ARG1 ) / 3 . 

ARG1=1./ GAMMA ( I ) 

ARG1=ARG1**2 
ARG2=1.+GAMMA( I ) **2 
ARG2  =  S0RTF( ARG2  ) 

ARG21 =1 .  /ARG2 
ARG3  =  L0GF(GAMMA(  I  )  ) 

ARG4  =  ARG1 /GAMMA (  I  ) 

ARG5=ARG4*ARG1 

U1  (  I  )  =C_D  *(  (  XO*GAMMA  (  I  )*L0GF  (  1.+ARG2  )  +  (  Xl-XO  )  *GAMMA  (  I  )*LOGF  (  GAM 
IMA ( I ) )+A3*GAMMA(  I  )+X2*GAMMA( I  ) ** 3+X 3*GAMMA ( I ) ** 5-X 1 /GAMMA ( I )  )*ARG2 
2 1 + A 1 *GAMMA ( I )+X0 /GAMMA! I )  ) / 3 . 

UFLEXSI  I  )  =U1 ( I  )-CD*A3*GAMMA( I )*ARG21/3. 

Will) =CD  * (  ( XO*LOGF ( 1 . +ARG2 )  +  ( XI -XO ) *LOGF ( GAMMA ( I ) )+X2*GAMMA (  I  ) 
1 **2+X  3*G AMMA ( I  ) **4  +  A3-Xl/GAMMA( I  ) **2 ) *ARG2 1  +  ( X4  +  X5*GAMMA ( I  )**2  +  Xl/ 
2GAMMA ( I ) **2 )*ARG2-X0 ) /3 . 

WFLEXSI 1 ) =W 1 ( I )-CD*A3*ARG21/3. 

11  CONTINUE 

UB 1 =ABSF (U1 (1) ) 

UB2=ABSF (U1 (M)  ) 

IF(UB1-UB)59.59.62 

59  U 1 (  1  )  =0  . 

62  IF (UB2-UBI63.63.80 
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C 

C 

c 


63  U1(M)=0. 


-  -  (7 )  OUTPUT  I  I 


80 

225 

226 


227 

228 

33 

36 


48 

61 

13 

14 

15 


84 

85 

86 


34 

57 


60 


44 


46 


IF (N) 225*225*226 
I  F ( PS  I  ) 38*226*22  7 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
IF  (N) 228*228*227 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
I F ( PS  I  ) 38*33*34 
I F ( NORM ) 36*36*61 
DO  48  I  =  1  * M 
EN1  (  I  ) =  E N 1 ( I  ) /CS 


3*229 

3*231  *  A  1 ♦ A  3 
3*114 

3*230 

3*231 *C1*C2*C3*C4 
3*107 


EN2(  I  )  =  EN2 ( I ) /CS 
UFLEXS (  I  )=UFLEXS(  I  ) /CD 
WFLEXSI  I  )  = WFLE XS (  I  ) /CD 


W 1  (  I  )  = W 1(1) /CD 
U 1  (  I  )  =U 1 (  I ) /CD 
WRITE  OUTPUT  TAPE 
I F ( NORM  )  13*13*14 
WRITE  OUTPUT  TAPE 
GO  TO  15 

WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
I F ( NORM ) 84*84  *85 
WRITE  OUTPUT  TAPE 
GO  TO  86 

WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 


3*109 

3*310 

3*110 

3*811  ♦  (GAMMR  (I)*Ul(I)*Wim*ENl(I)*EN2(I).I  =  l*M) 
3*107 

3*348 

3*148 

3*144* (GAMMR (  I  )  *UFLEXS(  I  )  *WFLEXS (  I  ) *  I = 1  *  M ) 


WRITE  OUTPUT  TAPE  3*107 

GO  TO  10 

IF(N)57, 57*56 

DO  60  I • 1  *  M 

A  (  I  ) = A  (  I  ) +  EN  1 (  I  ) *CPH I 

0 ( I ) =B ( I ) +EN2 ( I ) *CPH I 

WFLEXA (  I  )  =WFLEX A (  I  )+WFLEXS( I  ) *CPHI 

UFLEXA (  I  )=UFLEXA(  I  )+UFLEXS( I  )*CPHI 


P ( I ) =P ( I )+Wl( I ) *CPH I 
Q(  I  ) =0 (  I  )+Ul(  I  )*CPHI 
CNORMs 1 . 

I F ( NORM ) 44  *  44  *  45 
DO  46  I =1 *M 
A (  I  ) = A (  I  )/C$ 

R ( I  ) =R (  I  )  /CS 

C(  I  ) =C (  I  )/CS 

G ( I ) =G ( I ) /CD 

P (  I  )  =P (  I  ) /CD 

0(  I  ) =Q(  I  )  /CD 

UFLEXA (  I  ) =UFLE XA (  I  )/CD 

VFLEX A (  I  )=VFLEXA(  I  ) /CD 

WFLEXA(  I  )=WFLEXA(  I  ) /CD 

CONTINUE 

CNORM=CD 
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TABLE  1  (Continued) 


45 

16 

17 

18 


42 

IP 

20 

21 


87 

88 
89 


43 


4  P 


22 

23 

24 


200 

201 

202 

56 


3 

12 


64 


IF (NORM) 16*16*17 
WRITE  OUTPUT  TAPE 
GO  TO  18 

write  OUTPUT  TAPE 
write  output  tape 

WRITE  OUTPUT  TAPE 
IFOCD-l  )  42  *  4  2  *  4  3 
I F ( NORM ) 19*19*20 
WR I TE  OUTPUT  TAPE 
GO  TO  21 

WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
IF (NORM ) 87  *  87  *88 
WRITE  OUTPUT  TAPE 
GO  TO  89 

WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
IF ( KD-1  )  10*43  *43 
DO  49  I =1  *M 


3*303 

3*103 

3*106*(S(  I  ) *A(  I  ) *  8 (  I  )  *  C  (  I )*I=1*M) 
3*107 


3*304 

3*104 

3*812*(S(I)*0(I)*G(I)*P(I)*I=1*M) 

3*107 

3*349 

3*149 

3* 106  *  ( S(  I)  *UFLEXA (  I  )  * VFLEXA (  I  ) * WFLEX A (  I  ) *  I  =  1 *M ) 
3*107 


P ( I ) =P ( 1 )-Wl ( I )  / CNORM 

0(11=0(1  )-Ul (  I  )  /CNORM 
WFLEX A (  I  ) =WFLEX A (  I  ) -WFLEXS ( I  ) /CNORM 
UFLEXA(  I  )=UFLEXA(  I  )-UFLEXS(  I  1 /CNORM 
CONTINUE 

WRITE  OUTPUT  TAPE  3*142 
IF(NORM)22*22*23 
WRITE  OUTPUT  TAPE 
GO  TO  24 

WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
I F ( NORM ) 200*200*201 
WRITE  OUTPUT  TAPE  3*440 
202 

OUTPUT  TAPE  3*240 

OUTPUT  TAPE  3*106* (S( I  1  *UFLEXA(  I  )  *VFLEXA(  I  ) *WFLEXA(  I  ) 


3*439 

3*239 

3*812*(S( I  )  *Q(  I  ) *  G (  I  ) *P(  I  )* I=1*M) 
3*107 


GO  TO 
WRITE 
WRI  TE 
GO  TO 


1=1 *M) 


10 


DO  3  J=1 * N 


DO  3  1  =  1  *  M 

V (  I  *  J  )  =G (  I )  *COSF ( THETA ( J)  ) 

W (  I  * J ) =  SINF(THETA(J)  )*P(  I  )  +W 1(1)* CP HI 

U (  I  *  J  )  =  SINF(THETA(J)  )*Q(  I  )  +  U1  (  I  )*CPHI 

VFLEX (  I  *  J ) =  VFLEX  A (  I  ) *COSF ( THE T A ( J )  ) 

WFLEX ( I  *J)=WFLEXA(  I  ) * S I NF ( THE T A ( J  )  )+WFLEXS(  I  )*CPHI 

UFLEX ( I*J)=UFLEXA(I)*SINF(THETA(J) )+UFLEXS( I )*CPHI 

ENR (  I  * J ) = A (  I  )  *SINF ( THETA ( J)  )  +  EN1 (  I  )*CPHI 

ENT (  I  * J  )  =B ( I  )  *S INF ( THETA ( J )  )  +EN2 (  I )*CPHI 

ENR T (  I  *  J  )  =C (  I  ) *COSF( THETA ( J)  ) 

CONTINUE 

DO  12  J  =  1  *  N 

THETA(J )=THETA(J)*RAD 


SK  I  P  =  0 


CNORM  =  1  . 

IF(N0RM)64, 64*65 
DO  66  I =1  *M 
DO  66  J  =  1  *  N 
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TABLE  I  (Continued) 

W(  I  * J ) =W ( I ,J)/CD 

U(  I  * J ) =U (  I .J) /CD 

V< I . J ) =V ( I ,J) /CD 
WFLEX ( I »J)=WFLEX(  I 

♦  J)  /CD 

UFLEX ( I  *  J ) =UFLEX (  I 

.J) /CD 

VFLEX ( I  * J ) =VFLEX (  I 

f J) /CD 

ENR (  I , J  )  =  ENR (  I  * J ) /CS 

ENT ( I .J 1 =ENT ( I . J ) /CS 

ENRT (  I  *  J  )  =ENRT ( I.Jl/CS 

66 

CONTINUE 

CNORM=CD 

WRITE  OUTPUT  TAPE 
GO  TO  25 

3*330 

65 

WRITE  OUTPUT  TAPE 

3.130 

25 

write  OUTPUT  TAPE 

3.FMT1.(THETA(J).J=1.N) 

WRITE  OUTPUT  TAPE 

3*114 

WRITE  OUTPUT  TAPE 

3 . FMT  2 . ( GAMMR (I).(ENR(I,J),J  =  1,N), 

1=1 *M) 

WRITE  output  tape 

I F ( NORM ) 26*26*27 

3*107 

26 

WRITE  OUTPUT  TAPE 
GO  TO  28 

3*331 

27 

WRITE  OUTPUT  TAPE 

3*131 

28 

WRITE  OUTPUT  TAPE 

3 . FMT 1 , ( THETA! J) .J=l .N) 

WRITE  OUTPUT  TAPE 

3*114 

WRITE  OUTPUT  TAPE 

3* FMT 2  * (GAMMR (I)*(ENT(I*J)*J=1*N)* 

1=1 *M) 

WRITE  OUTPUT  TAPE 

I  F ( NORM ) 29*29*37 

3*107 

29 

WRITE  OUTPUT  TAPE 
GO  TO  41 

3*332 

37 

WRITE  OUTPUT  TAPE 

3*132 

41 

WRITE  OUTPUT  TAPE 

3  *  FMT 1  * ( THE T  A ( J ) * J=1 *N ) 

WRITE  OUTPUT  TAPE 

3*114 

WRITE  OUTPUT  TAPE 

3*FMT2* (GAMMR (  I  )  * ( ENRJ ( I  ,  J )  * J*1 *N ) 

* I=1*M) 

WRITE  OUTPUT  TAPE 
IFIKD-l 182.82.83 

3*107 

82 

IF (NORM ) 501.501.504 

501 

IFISKIP1502.502.503 

502 

WRITE  OUTPUT  TAPE 
GO  TO  507 

3*333 

503 

WRITE  OUTPUT  TAPE 
GO  TO  507 

3*433 

504 

IF (SKIP ) 505*505*506 

505 

WRITE  OUTPUT  TAPE 
GO  TO  507 

3*133 

506 

WRITE  OUTPUT  TAPE 

3*233 

507 

WRITE  OUTPUT  TAPE 

3  *  FMT 1  * ( THETA (J ) * J=1 *N) 

WRITE  OUTPUT  TAPE 

3*114 

WRITE  OUTPUT  TAPE 

3* FMT 2* (GAMMR ( I)  *  ( U ( I . J )  * J= 1  * N  )  *  I  = 

1  *  M  ) 

WRITE  OUTPUT  TAPE 

3*107 

I F ( NORM  1508.508.511 

508 

IF (SKIP  1 509.509.510 

509 

WRITE  OUTPUT  TAPE 
GO  TO  514 

3*334 

510 

WRITE  OUTPUT  TAPE 
GO  TO  514 

3*434 

511 

IFISKIP 1512.512.513 

512 

WRITE  OUTPUT  TAPE 
GO  TO  514 

3*134 
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TABLE  I  (Continued) 


513  WRITE  OUTPUT  TAPE  3,23A 

51A  WRITE  OUTPUT  TAPE  3,FMT1,(THETA(J),J=1,N) 

WRITE  OUTPUT  TAPE  3,11A 

WRITE  OUTPUT  TAPE  3 , FMT2 . ( GAMMR U ) , ( V U , J ) , J=1 .N ) . I = 1 , M ) 

WRITE  OUTPUT  TAPE  3,107 

IF (NORM  1515,515,518 

515  IF(SKIP)516*516*517 

516  WRITE  OUTPUT  TAPE  3,335 
GO  TO  521 

517  WRITE  OUTPUT  TAPE  3.A35 
GO  TO  521 

518  IF <SK  IP  )  519,519,520 

519  WRITE  OUTPUT  TAPE  3,135 
GO  TO  521 

520  WRITE  OUTPUT  TAPE  3,235 

521  WRITE  OUTPUT  TAPE  3  ,  FMT 1 , ( THE T A ( J )  ,  J  =  1  ,  N  1 
WRITE  OUTPUT  TAPE  3,11A 

WRITE  OUTPUT  TAPE  3 , FMT2 , ( GAMMR (  I  )  , ( W (  I , J )  , J  =  1 ,N ) , I  =  1 , M ) 
WRITE  OUTPUT  TAPE  3,107 
IF ( NORM ) 522,522,525 
52?  IF (SKIP ) 523, 523.52A 
523  WRITE  OUTPUT  TAPE  3,350 
GO  TO  528 

52 A  WRITE  OUTPUT  TAPE  3,A36 
GO  TO  528 

525  IF ( SKIP ) 526,526, 527 

526  WRITE  OUTPUT  TAPE  3,150 
GO  TO  528 

527  WRITE  OUTPUT  TAPE  3,236 

528  WRITE  OUTPUT  TAPE  3 , FMT 1 , ( ThET A ( J ) , J= 1 > N ) 

WRITE  OUTPUT  TAPE  3,11A 

WRITE  OUTPUT  TAPE  3 , FMT2 > ( GAMMR ( I  ) , ( UF LEX (  I  , J  )  ,  J  =  1 ,N ) , I = 1  ,M ) 
WRITE  OUTPUT  TAPE  3,107 
I  F ( NORM  )  529,529,532 

529  IF ( SK IP ) 530,530 ,531 

530  WRITE  OUTPUT  TAPE  3,351 
GO  TO  535 

531  WRITE  OUTPUT  TAPE  3,A37 
GO  TO  535 

532  IF(SKIP)533, 533.53A 

533  WRITE  OUTPUT  TAPE  3,151 
GO  TO  535 

53A  WRITE  OUTPUT  TAPE  3,237 

535  WRITE  OUTPUT  TAPE  3 , FMT 1 , ( THET A ( J  1  » J  =  1  ,  N  ) 

WRITE  OUTPUT  TAPE  3,11A 

WRITE  OUTPUT  TAPE  3 » FMT2 , ( GAMMR (  I  )  ,  ( VFLEX (  I , J )  * J  =  1  ,N ) , I  =  1 ,M ) 
WRITE  OUTPUT  TAPE  3,107 
I  F ( NORM  1536,536,539 

536  IF ( SK IP  1 537,537 ,538 

537  WRITE  OUTPUT  TAPE  3,352 
GO  TO  5  A  ? 

538  WRITE  OUTPUT  TAPE  3.A38 
GO  TO  5A? 

539  IF(SKIP)5A0,5A0,5A1 

SA n  WRITE  OUTPUT  TAPE  3,152 
GO  TO  5  A 2 

5 A 1  WRITE  OUTPUT  TAPE  3,238 

5 A 2  WRITE  OUTPUT  TAPE  3 , FMT 1 , ( THET A ( J 1 , J  =  1 , N 1 
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215 

83 


223 

224 


WRITE  OUTPUT  TAPE  3*114 

WRITE  OUTPUT  TAPE  3  * FMT 2  * ( GAMMR (  I  ) * ( WF LEX (  I  * J )  * J = 1  * N ) *  I = 1  * M ) 

WRITE  OUTPUT  TAPE  3*107 

IF(SKIP)215*215*10 

IFOCD-l ) 10*83*83 

DO  224  J  =  1  *  N 

DO  223  1=1, M 

UFLEX ( I  * J)=UFLEX(  I  ♦  J ) — U 1 (  I )/CNORM 
WFLEX (  I  , J) =WFLEX (  I  * J ) -W1 (  I ) /CNORM 
U(I*J)=U(I*J)-U1(I ) /CNORM 
W(  I*J)=W(I*J)-W1(  I  ) /CNORM 
CONTINUE 

WRITE  OUTPUT  TAPE  3*142 
SK  I  P= 1 
GO  TO  82 


C 

C  -  - 


(8 )  RETURN  TO  ( 1 ) 


101 

FORMAT ( 

8  15) 

102 

FORMAT ( 

5  E  14  *  8  ) 

103 

FORMAT ( 

1  16H 

R  ( 

IN.  ) 

NR ( LB. / I N . ) 

1 

NTHETA(LB./ IN. ) 

NRTHETA! LB. /I N. )  //// ) 

104 

FORMAT 

(  104H 

R  ( 

IN.  ] 

l  U ( IN.) 

1 

V (  IN*) 

W ( IN*)  ////) 

105 

FORMAT 

( 4F15.9) 

106 

formati 

:  F25 • 4  *  3F2  5 • 8 ) 

107 

format 

(  1 H  1  ) 

108 

FORMAT  1 

[ 4F24.5*F18 .2) 

109 

FORMAT 

(75H  THE  LOADING  AS 

WELL 

AS 

THE  DEFORMATION  OF  THE  SHELL  I 

IS  AXISYMMFTRIC 

110 

format ( 

122H 

1 

2)////) 

W  (  IN, 

111 

FORMAT ( 

5F24.5) 

112 

FORMAT ( 

1 2  OH 

12  ) 

2///) 

THICKNI 

113 

FORMAT ( 

1  16H 

1 

R  3 (  IN, 

114 

FORMAT 

( //// ) 

115 

FORMAT 

(  50H 

116 

FORMAT 
IT  R 1 

(  70H 

117 

FORMAT ( 

5  1  H 

118 

FORMAT ( 

75H  T) 

1 ALLY  AT 

R2 

119 

FORMAT ( 

2F24.2) 

130 

FORMAT ( 

5  OH 

131 

FORMAT ( 
1//) 

53H 

132 

FORMAT ( 
1///) 

54H 

133 

FORMAT ( 

5  OH 

134 

FORMAT ( 

5  OH 

135 

FORMAT ( 

5  OH 

139 

FORMAT! 

5  1 H - 1 

140 

FORMAT ( 

5  1 H - 1 

142 

FORMAT ( 

95H  THE 

////  ) 


R (  IN*) 

NR ( LB*/ IN* ) 


WEIGHT  DENSI TY ( LB./ I N.2 )  YOUNGS 

>S( IN. I  FOCAL  LENGTH! IN. ) 

R  1  (  IN.  ) 

)  R4 (  I N  .  ) 


U (  IN.  ) 

NTHETA(LB. /IN. 


MODULUS ( LB • / I N • 
POISSONS  RAT  10/ 

R2( IN. ) 

PSI (DEG. )////) 


THE  SHELL  IS  FIXED  TANGENTIALLY  AT  BOTH  EDGES  ///) 
THE  SHELL  IS  FIXED  TANGENTIALLY  AT  R2  AND  IS  FREE  A 
////  ) 

THE  T  A3 ( DEG • )  THETA4 ( DEG.  )////) 


SHELL  IS 
////) 


CLOSED  AT  THE  APEX  AND  IS  FIXED  TANGENTI 


-STRESS 

-STRESS 


RESULTANT 

resultant 


NR ( LB. / I N . ) - 

NTHETA(LB./IN.L 


///// ) 
- /// 


-STRESS  RESULTANT  NRTHET A ( LB. / I N . ) « 


-// 


—DISPLACEMENT  U(  IN.  ) - 

—DISPLACEMENT  V( IN. ) - 

—DISPLACEMENT  W( IN. ) - 

RESULTS  BY  MEMBRANE  ANALYSIS 
:ED  RESULTS  BY  MEMBRANE  ANALYSIS 


/////) 
/////) 
/////) 
////  ) 
////  ) 


DISTORTION  OF  THE  SHELL  GIVEN  BELOW  IS  MEASURED  REL 
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TABLE  I  (Continued) 


iative  to  that  of  the  face-up  position  ////) 

143  FORMAT  (  55H  the  PROGRAM  DOES  NOT  ACCEPT  A  NEGATIVE  POINTING  ANGLE 
1  //// ) 

144  FORMAT (F24.4*2F24.8) 


[N.  ) 


148  FORMAT ( 1 16H  R ( IN. )  UFLEXUN.) 

1  W  FLE  X (  IN* )  ////) 

149  FORMAT ( 1 1 6H  R ( IN. )  UFLEXUN.) 

1  VFLEX ( 

150  FORMAT ( 50H 

151  FORMAT ( 50H 

152  FORMAT ( 50H 

160  FORMAT ( 56H  A  HOMOGENEOUS  PARABOLOIDAL  SHELL  SUBJECTED  TO 
1  ////  ) 

161  FORMAT ( 56H  A  LAMINAR  PARABOLOIDAL  SHELL  SUBJECTED  TO  GRAVITY 


R(  IN.  ) 

WF  LEX (  IN.  ) 

- DISPLACEMENT  UFLEX (  IN.  ) - 

- DISPLACEMENT  VFLEX (IN.) - 

- DISPLACEMENT  WFLEXUN.) - 

HOMOGENEOUS  PARABOLOIDAL  SHELL  SUBJECTED 


/ ///) 
/////) 
/////) 
/////) 
GRAVITY 


1  //// ) 

162  FORMAT ( 56H  A  SANDWICH  PARABOLOIDAL  SHELL  SUBJECTED  TO  GRAVITY 
1  //// ) 

163  FORMAT ( 1 1 6H  STRETCHING  S T I F F NESS ( LB / I N**2 )  WEIGHT  DENSITY! LB/ IN* 

1*2)  FOCAL  LENGTH ( IN. )  POISSONS  RATIO  ////) 

164  FORMAT (  1  16H  YOUNGS  MODULUS ( LB/ I N** 2 )  CORE  TH I CKNESS (IN.) 

1  SKIN  THICKNESS! IN. )  RHO  OF  SK I N ( LB/ I N** 3 )  ////) 

165  F ORM A T (  1 1 6H  RHO  OF  CORE ( LB/ I N**3 )  WEIGHT  DEN S I T Y ( LB / I N* 

1*2)  FOCAL  LENGTHUN.)  POISSONS  RATIO  ////) 

205  pORMAT ( 1 2A6 ) 

22  9  FORMAT! 19X*2HA1 *  2 3X ♦ 2HA2 ♦ / / / / ) 

230  FORMAT! 1QX*2HC1«23X*2HC2*23X*2HC3*23X*2HC4*//// ) 


231  FORMAT ( 4F25. 6) 

233  FORMAT  (  5  OH  DISPLACEMENT  U  TILDEUN.) -  /////) 

234  FORMAT  (  5  OH  DISPLACEMENT  V  TILDEUN.) -  /////) 

235  FORMAT ( 50H  - DISPLACEMENT  w  TILDE! IN. ) -  /////) 

236  FORMA  T  (  5  OH  DISPLACEMENT  UFLEX  TILDEUN.) - /////) 

237  FORMAT  (  5  OH  DISPLACEMENT  VFLEX  TILDEUN.) - /////) 

238  F  ORM  A  T  (  5  OH  DISPLACEMENT  WFLEX  TILDEUN.) - /////) 

239  FORMAT  (  1  1  6H  RUN.)  U  TILDEUN.) 

1  V  TILDEUN.)  W  TILDEUN.)  ////) 

240  FORMAT  (  1  1 6H  RUN.)  UFLEX  TILDFUN.) 

1  UFLEX  TILDEUN.)  WFLEX  TILDEUN.)  ////) 

303  FORMAT (116H  GAMMA  NR* 

1  NTHETA*  NRTHETA*  ////) 

304  FORMAT  ( 1 04  H  GAMMA  U* 

1  V*  W*  ////) 

310  FORMAT! 121H  GAMMA  U* 


1  W* 

2////) 

330  FORMAT (75H 

1LESS) - 

331  FORMAT ( 7  5H 

1SIONLESS) - 

332  FORMAT ( 7  5H 

INSIONlESS) - 

333  FORMAT ( 60H 


NR* 


NTHETA* 


- NORMALIZED 

///// ) 

- normalized 

/////) 

- NORMAL  I  ZED 


STRESS  resultant 
STRESS  resultant 
STRESS  RESULTANT 


NR* ( D I  MEN  S I  ON 
NTHETA*(DIMEN 
NR  THE  T  A* ( D I  ME 


/////) 

- NORMALIZED  DISPLACEMENT  U*(DIMENSIONLESS) 


1 - /////) 

334  FORMAT (60H 

1 - /////) 

335  FORMAT ( 6 OH 


- NORMALIZED  DISPLACEMENT  V*(DIMENSIONLESS) 

- NORMALIZED  DISPLACEMENT  W*(DIMENSIONLESS) 


1 - /////) 

348  FORMAT! 1 16H 

1  WFLEX* 


GAMMA 


UFLEX* 


////  ) 
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TABLE  I  (Continued) 


349  FORMAT ( 1 16H 

GAMMA 

UFLEX* 

1  VFLFX* 

wflex* 

////  ) 

350  FORMAT (65H 

- NORMALIZED 

displacement 

UFLEX*  ( DIMENS  I ONL 

1 FSS -  /////) 

351  FORMAT ( 65H 

- NORMALIZED 

displacement 

VFLEX* ( D I  MENS  I ONL 

1 FSS -  /////) 

352  FORMAT (65H 

- normalized 

displacement 

WELEX* ( D I  MENS  I ONL 

1 ESS -  /////) 

433  FORMAT ( 7  5H 

- NORMALIZED 

displacement 

U  T  ILDE*(DIMENSIO 

INLESS) - 

/////) 

434  FORMAT ( 7  5H 

- NORMALIZED 

displacement 

V  T  I  LDE  * ( D I  MENS  I 0 

INLESS) - 

/////) 

435  FORMAT ( 75H 

- NORMAL  I  ZED 

displacement 

W  T  I  LDE* ( D I  MENS  I 0 

INLESS ) - 

/////) 

436  FORMAT ( 7  5H 

- NORMAL  I  zed 

displacement 

UFLEX  TILDE*  (DIM 

1ENSIONLESS) - 

/////) 

437  FORMAT ( 7  5H 

- NORMAL  I  ZED 

displacement 

VFLCX  TILDE*  (DIM 

1ENSIONL  ESS ) - 

/////) 

438  F ORMA T ( 7  5H 

- NORMAL  I  ZED 

displacement 

WFLEX  TILDE*  (DIM 

1ENSIONLESS ) - 

/////) 

439  FORMAT ( 1  16H 

GAMMA 

U  TILDE* 

1  V  TILDE* 

W  TILDE* 

////  ) 

440  FORMAT ( 1 16H 

gamma 

UFLEX  TILDE* 

1  VFLEX  TILDE* 

WFLEX 

TILDE* 

////  ) 

811  F0RMAT(r24.4,L24.6,3F24.8) 
81?  FORMAT  (  F25.4*E25*6*2F25*8  ) 
FND 

*  DATA 
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1 .  Input 


The  input  to  our  program  as  indicated  by  step  (1)  in  Fig.  12  consists  of  a  number  of  records 
prestored  on  tape  A2. 

Record  1  Control  Parameters  (615) 

This  line  contains  six  non-negative  fixed-point  variables  (integers). 

I  M  I  N  I  KE  I  NORM  I  KD  I  KP  I 

Number  of  (equally  spaced)  y  —  tabular  points  between  R3  and  R4 
(R3  and  R4  are  to  be  defined  later)  for  any  fixed  value  of  0. 

Number  of  (equally  spaced)  0  —  tabular  points  between  THETA3  and 
THETA4  (THETA3  and  THETA4  are  to  be  defined  later)  for  any 
fixed  value  of  y.  If  only  those  quantities  with  superscript  zero 
are  desired  [see  Eqs.  (11-85)  and  (11-86)],  set  N  equal  to  zero. 

Control  parameter  specifying  the  edge  conditions  at  R1  and  R2 
(R 1  and  R2  will  be  defined  later.  They  are  not  to  be  confused 
with  r^  and  r^  which  are  the  lower  and  upper  edge  of  the  shell). 

KE  <  1  Shell  is  fixed  tangentially  at  both  edges. 

KE  =  1  Shell  is  fixed  tangentially  at  R2  and  is  free 

at  R1 . 

KE  >  1  Shell  is  closed  at  the  apex  and  is  supported 
tangentially  at  R2. 

Control  parameter  specifying  whether  normalized  results  are 
desired. 

NORM  =  0  Normalized  results. 

NORM  >  0  Un-normalized  results. 

Control  parameter  specifying  the  type  of  distortion  to  be  given 
in  the  output. 

KD  <1  u,  v,  w,  and  the  corresponding  u^ex-  x"  Wflex 
are  given  in  the  output. 

KD  >  1  u,  v,  w,  and  the  corresponding  uflex<  vflex,  wflgx 
are  given  in  the  output. 

KD  -  1  u,  v,  w,  u,  v,  w,  as  well  as  the  corresponding  six 
flex  quantities. 

If  N  -  0,  a  superscript  zero  should  be  added  to  these  displacement  quantities. 

KP  Control  parameter  specifying  the  type  of  shell  under  consideration. 

Set 

KP  <1  if  the  shell  is  homogeneous, 

KP  >1  if  the  shell  is  a  sandwiched  construction 

with  a  very  soft  core, 

KP  =1  if  the  shell  is  laminar  in  that  its  material 
properties  vary  across  the  thickness  in  a 
manner  other  than  previously  mentioned. 


M  (<  150) 
N  (,$10) 

KE 


NORM 

KD 
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Uocord  2  Geometrical  Parameters  ( 4 L*'  15.9) 

This  line  contains  four  non-negative  F-type  floating  point  variables. 

HI  '  K2  T  K  3  1  K4  1 

HI  (in.)  Value  of  r  at  the  free  edge  of  the  shell.  If  the  shell  is  closed 
at  the  apex,  then  HI  0.  If  the  shell  is  fixed  tangentially  at 
both  edges,  then  set  HI  equal  to  the  value  of  r  at  the  lower  edge. 

HZ  (in.)  Value  of  r  at  the  tangentially  fixed  edge.  If  both  edges  of  the 

shell  are  fixed  tangentially,  set  HZ  equal  to  the  value  of  r  at  the 
upper  edge. 

H 3  (in.)  Smallest  value  of  r  at  which  the  stresses  and  distortion  of  the 

shell  are  to  be  given  in  the  output.  Generally,  this  is  the  same 
as  r^.  However,  if  the  shell  is  closed  at  the  apex,  it  is  usually 
wise  to  avoid  calculating  the  limiting  value  of  the  desired  output 
at  the  apex  by  setting  R3  >  0. 

K4  (in.)  Largest  value  of  r  at  which  the  stresses  and  distortion  of  the 

shell  are  to  be  given  in  the  output.  Generally,  this  is  the  same 
as  r2. 

Record  3  Geometrical  and  Material  Parameters  (3F15.9) 

This  line  contains  three  non-negative  F-type  floating  point  variables. 

I  F  I  PR  I  PSI 

F  Focal  length  of  the  paraboloidal  surface  (in.). 

PH  Poisson’s  ratio. 

PSI  Pointing  angle  ip  in  fractions  of  tt  (e  g.,  if  4  =  45°,  then  psi  =  0.25). 


Record  4  Geometrical  and  Material  Parameters  (5E14.8) 


KP  <  1  This  line  contains  three  E-type  floating  point  variables. 

r  KHO  T  H  I  E  | 

RHO  Volume  weight  density  of  shell  (lb/in.  ). 

H  Shell  thickness  (in.). 

E  Young's  modulus  (lb/in2). 


KP  >  1 


This  line  contains  five  E-type  floating  point  variables. 

I  RHS  r  RHC  I  T  T  H  1  E  I 

RHS  Volume  weight  density  of  skin  (lb/in.  ). 

RHC  Volume  weight  density  of  core  (lb/in.^). 

T  Thickness  of  skin  (in.). 

H  Thickness  of  core  (in.). 

E  Young's  modulus  of  skin  (lb/in.“). 
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Kl*  t  This  line  contains  two  E-type  floating  point  variables. 

r  KUO  T  CC  1 

R!l()  Surface  weight  density  of  shell  (lb/in.  ). 

(T  Kxtensional  stiffness  of  shell. 

Record  5  Geometrical  Parameters  (2F15.9) 

This  line  contains  two  F-type  floating  point  variables. 

I  TIIETA3  I  THETA4  | 

T11ETA3  Smallest  value  of  G  in  fractions  of  v.  at  whieh  the  stresses 
and  distortion  of  the  shell  are  to  be  given  in  the  output. 

THETA4  Largest  value  of  G  in  fractions  of  ir  at  whieh  the  stresses 
and  distortion  of  the  shell  are  to  be  given  in  the  output. 

Record  6  Variable  Format  Statement  1 

This  line  provides  a  format  statement  for  listing  the  set  of  G-tabular  points  (see  also 
Sec.  lll-B-2),  for  example, 

|  ( 19H  R(IN.)/THETA(DEG  ).  NF11.2)  I  . 

Reeord  7  Variable  Format  Statement  11 

This  line  provides  a  format  statement  for  the  values  of  each  physical  quantity  at  different 
positions  (see  also  Sec.  I1I-B-2),  for  example, 

|~~(F19.4,  NF11.6)  |  . 

Remark—  If  N  0  or  if  PS1  =  0,  then  the  input  must  eonsist  of  only  the  first  four  records 

in  the  preceding  list.  Additional  reeords  will  be  treated  as  a  new  set  of  input  upon  return.  For 

the  purpose  of  an  Express  Run,  eaeh  of  these  reeords  is  punehed  out  on  one  eard.  A  complete 
set  of  input  will  then  eonsist  of  seven  cards  (or  four  eards  if  N  =  0  or  if  tp  -  0).  It  is  to  be  plaeed 
after'  the  DATA  card  which  in  turn  follows  the  object  deck.  Additional  runs  ean  be  made  by  merely 
placing  additional  sets  of  input  after  the  first  set. 

If  some  restriction  on  the  input  is  violated,  the  program  returns  to  step  (1)  to  take  in  a  new 
set  of  data.  It  may  or  may  not  give  a  statement  pointing  out  the  souree  of  trouble. 

2.  Output 

The  first  part  of  the  output,  given  by  step  (2)  in  Fig.  12,  reproduces  for  the  reeord  the  input 
to  the  program.  It  states  the  problem,  the  edge  conditions,  the  various  geometrical  and  material 
parameters  and  whether  the  numerical  results  have  been  normalized. 

The  second  part  of  the  output,  given  by  step  (7)  in  Fig.  12,  gives  the  values  of  the  six  physical 
quantities,  N^,  N(),  N^,  u,  v,  and  w  at  the  tabular  points.  If  4‘  7  0,  the  behavior  of  the  shell 
is  axisymmetric .  Thus,  v  £  N  =  and  Nf,  Nq»  u,  and  w  are  functions  of  y  only.  The  output 
in  this  case  is  given  in  one  block  of  five  columns  of  numbers.  The  first  column  lists  the  set  of 

y-tabular  points.  The  next  four  columns  list  the  values  of  u,  w,  N^,  and  in  that  order,  cor¬ 

responding  to  these  values  of  y. 

If  4>  -/  0,  the  stress  resultants  and  the  displacements  generally  depend  on  both  y  and  O.  The 
output  for  each  of  these  quantities  is  arranged  in  a  rectangular  bloek  so  that  eaeh  eolumn  of  data 
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corresponds  to  the  values  of  the  quantity  for  fixed  value  of  9  and  for  different  values  of  y,  and 
each  row  corresponds  to  the  values  of  the  quantity  for  a  fixed  y  and  for  different  values  of  9. 
Inasmuch  as  N  is  not  constant,  the  number  of  columns  may  vary  for  different  runs.  Therefore, 
if  N  >  0,  two  variable  format  statements  are  needed  for  the  output.  The  first  of  these  is  used 
to  list  all  possible  values  of  9  as  the  first  record  of  the  block.  The  second  is  used  repeatedly 
to  list  records  of  (N  +  1)  numbers.  The  first  number  of  each  record  is  a  value  of  y  while  the 
remaining  numbers  are  the  N  values  of  the  physical  quantity  in  question  corresponding  to  the  N 
positions  given  by  this  value  of  r  and  the  N  9 -tabular  points  listed  in  the  first  record.  Examples 
given  in  a  later  section  will  clarify  this  discussion.  Having  listed  all  the  stress  and  displacement 
quantities,  the  next  block  of  data  gives  the  maximum  values  (in  modulus)  of  u^  ,  vfqex>  anc* 

wflex  °r  Uflex’  ^fiex*  an<^  ^flex  ^or  set  ta^u-^ar  V's- 

If  N  -  0,  so  that  only  the  zero-superscripted  quantities  are  requested,  this  second  part  of 

the  output  will  be  arranged  in  three  blocks  of  data  each  consisting  of  four  columns.  The  first  of 

these  columns  in  each  block  lists  the  set  of  y-tabular  points.  The  remaining  three  columns  of 

the  first  block  list  N  °,  N°  and  N  °_  _  those  of  the  second  block  list  u°,  v°,  and  w°,  and  those 
r  9  r9' 

of  the  third  list  ,  v^c^,  and  wf°ex  f°r  the  given  set  of  tabular  points. 

The  above  discussion  applies  only  when  KD  <  1.  For  KD  >  1,  u  ,  v  ,  and  w  and  the  cor¬ 
responding  ujqex>  v^ex>  and  wQex  will  take  the  place  of  u°,  v°,  and  w°  and  the  corresponding 
uflex’  vHex'  anc*  wflex‘  ^  KD  =  1,  additional  blocks  of  data  will  appear  for  the  obvious  reason. 
Keeping  in  mind  the  discussion  in  the  earlier  sections,  the  output  in  this  case  is  self-explanatory. 

The  numerical  results  in  the  output  are  generally  printed  in  the  F-type  floating  point  format. 
There  are  a  few  exceptions  due  to  the  size  of  the  numbers  involved.  In  these  exceptional  cases, 
the  E-type  floating  point  format  is  used. 

3.  Operational  Information 

The  main  program  is  coded  in  FORTRAN  language  for  a  32K  IBM  7090  (and  IBM  7094)  and 
is  compiled  by  the  FORTRAN  II  compiler.  It  requires  no  routines  other  than  those  appearing 
as  library  routines  on  the  Lincoln  Laboratory  library  tape.  It  should  be  cautioned  that  these 
may  be  different  from  the  routines  under  the  same  name  used  elsewhere.  The  correspondence 
between  logical  tape  units  and  machine  tape  units  is  also  different  at  Lincoln  Laboratory.  Ap¬ 
proximately  11,000  storage  locations  have  not  been  touched  by  the  program.  Neither  sense 
switch  nor  sense  light  is  used.  Some  of  the  subscripted  variables  and  their  analytical  counter¬ 
parts  are  presented  in  Table  II. 

It  should  be  understood  that  if  normalized  results  are  requested,  the  stress  and  displace¬ 
ment  expressions  in  the  third  column  of  Table  II  should  be  replaced  by  the  corresponding  starred 
quantities. 

C.  EXAMPLES 

1.  Axisymmetric  Deformations 

Consider  a  sandwich  shell  supported  at  the  lower  edge  r^  in  a  face-up  position  (4  =  0)  with 

rhe  =  0.01  lb/ in? 
rhs  =  0.1  lb/in.3 
E  =  107  lb/in.2 
v  =  0.3 
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TABLE  11 

SOME  SUBSCRIPTED  VARIABLES 

AND  THEIR  ANALYTICAL  COUNTERPARTS 

Subscripted  Variable 

Dimensionality 

Analytical  Expression 

U 

150  X  10 

u 

V 

150  X  10 

V 

W 

150  X  10 

w 

UFLEX 

150  X  10 

Uflex 

VFLEX 

150  X  10 

Vflex 

WFLEX 

150  X  10 

wr, 

flex 

ENH 

150  X  10 

N 

r 

ENT 

150  x  10 

Ne 

ENRT 

150  X  10 

Nre 

U1 

150 

S 

U 

W1 

150 

s 

w 

G 

150 

a 

V 

P 

150 

a 

w 

Q 

150 

a 

u 

ENl 

150 

NS 

r 

EN2 

150 

NS 

1  0 

A 

150 

Na 

r 

B 

150 

Ne 

C 

150 

Nr0 

UFLEXA 

150 

Uflex 

VFLEXA 

150 

a 

vflex 

WFLEXA 

150 

wn 

flex 

UFLEXS 

150 

uflex 

VFLEXS 

150 

vflex 

WFLEXS 

150 

s 

Wn 

flex 

GAMMR 

150 

r 

GAMMA 

150 

y 

THETA 

150 

e 
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i*t  100  in. 

v  y  300  in. 

f  500  in. 

h  0.5  in. 

t  0. 1  in. 

Since  the  deformation  is  axisy  mmetrie,  we  set  N  =  0  and  omit  records  5,  6,  and  7.  For  the 
purpose  of  illustration,  we  consider  only  five  y-tabular  points  (M  5).  The  output  will  not  he 
normalized  so  that  we  get  the  actual  physical  quantities.  Insofar  as  displacement  quantities  are 
concerned,  we  want  only  u  and  w  (and  of  course  the  corresponding  flex  displacements). 


INPUT 


Record  1  Control  Parameters  (615) 


M 

5 

I 


6  Z  (1  1  0  0  <3  0  KMMQOpmQ 

j  1  (•  :  ■  -M" ,i  ' t 

11111,111111111  1  ii  i 

2  7  2  22  7  7  7^2  7  2  1  1  2  7  7  7  7  7 

s  3  3  3  3  3  3  1  3  3  3  3  3  1  :>  J  3  1  3  3 

i  i  4  4  4  M  4  4  4  *,  ?4  *1  4  4  4  14  14  4  4 


0  9  0  0  3  0  0  0  0  0  3  0  0  0  C  0  G  0  0  0  li  9  0  0  C  2  0  u  3  5  0  0  0  0  9  0  0  0  0  C  C  C  ■  0  i  0 

.  >«?•  i  j»  rs it .t  ji  ::  jh  »%  r  jhhh  *•  >;  ou  *  «  ■  j.m  v  ■ :.« v  «• « :<  mi.  ,i  t>  „ »  a  r  • ,  • 

111111111  1  I  I  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  I  1  M  1  I  1  I 
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Record  2  Geometrical  Parameters  (4F15.9) 
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Record  3  Geometrical  and  Material  Parameters  (3F15.  9) 
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Record  4  Geometrical  and  Material  Parameters  (5E14.  8) 
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Remark:—  Since  PSI  is  zero,  the  last  three  records  must  be  omitted  according  to  input 
instructions. 

The  computer  output  is  presented  in  Table  Ill. 


2.  Unsymmetric  Deformations 

Consider  a  homogeneous  shell  supported  by  both  edges  and  oriented  in  such  a  way  that  i!  ^  0. 
We  have 

E  107  lb/in.^ 
v  =  0.3 
r  -  1 00  in. 

r2  =  300  in. 

f  =  500  in. 
h  =  0.01  in. 

rho  =  0.1  lb/in? 

We  now  want  normalized  outputs  for  seven  0 -tabular  points  in  the  interval  —  (7r / 2 )  <  0  <  (tt/2)  and 
again  five  y-tabular  points  in  the  interval  y  ^  ^  y  <  y^.  Insofar  as  displacements  are  concerned, 
we  will  limit  ourselves  to  u,  v,  w,  and  the  corresponding  flex  quantities.  (In  this  case,  that 
quantity  which  is  listed  in  the  output  as  the  weight  density  is  rho  x  h.) 
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2  Geometrical  Parameters  (4F15.  9) 
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Record  3  Geometrical  and  Material  Parameters  (3F15.  9) 
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Record  5  Geometricol  Porameters  (2F15.  9) 
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Record  6  Variable  Format  Statement  I  (FMT1) 
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Record  7  Variable  Format  Statement  II  (FMT2) 


(F19.4,7F12.6) 


ll  n  I)  I)  0  fl  3  0  0  5  0  0  0  G  0  0  0  0  0  0  0  0  0  C  0  0  J  0  I)  0  0  C  '  r  0  C  3  ,1  0  C  0  P  0  0  c  0  C  0  0  0  0  C  0  M  0  U  0  0  0  (I  0  C  0  C  C  G  0  0  0  0  u  0  Pi)  o  0  0 

•  s  ■  1  '  i  .  i  .  u  i  r.  i  .  j  « ;u ;  .'j ;  ■  .*  •. ,  .?•  w  •• .«  .  *  »  t  m  i:  »j  m  s  ?► '  ,  n  .'i  >/  :■■■.'  *  n  .•«*  w 

ii  n  i  1 1 1  1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ii 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ii 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1  n  1 1 1 1 1 1 1 ;  ii  i 

2  ?  222222  .2  <  2  2  2  2  2  2  2  2  2  2  2  7  2  /  2  2  1  2  2  2  2  2  2  n  2  7  7  2  ?  2  2  7  2  2  2  2  2  2  1  2  2  2  2  2  2  2  7  2  2  2  2  ?  2  2  2  2  2  2  2  2  7  2  1  2  2  2  : 

13  3  3  3  3  3  3  3  3  3  3  3  3  3  :  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3 

4-144  4  4  4  :  4  1  4  444  4  44444444444444444444444444  14444444444  *4444  ^4444444444444444444 

i  5  5  5  S  f»  S  5  5  b  ■)•»  b  *i  %  5  b  5  5  5  5  S  5  5  i»  5  i  i  5  5  5  5  1  S  5  S  5  5  5  5  5  5  5  'a  !»  S  5  5  5  5  5  5  'i  5  5  •>  fj  5  5  5  S  5  5  *>  $  5  5  5  $  5  5  5  5  5  5  !»  M»  5  5 

5  f  ;■  J  i  G  5  $  6  f,  6  5  6  6  5  B  6  6  6  6  5  3  6  6  6  6  &  6  f  6  6  b  6  6  6  6  6  6  S  $  *j  6  6  6  G  6  6  6  6  b  b  s  Ii  C  b  6  6  6  6  6  6  6  6  6  6  b  6  6  6  6  6  5  M.-  b  h 

mini  i  m  m  ;  n  i  mi  in  n  mm  m  in  mm  m  m  mm  i  mn  m  i  n  m  i  i  i  m  i 

$  i  8  4  3  S  3  8  3  3  ii  3  s  d  M  8  U  8  «  8  3  8  fi  8  8  8  8  8  8  8  8  3  8  8  8  8  8  8  8  8  8  8  3  8  8  B  $  £  8  »  8  3  8  t  C  8  8  8  3  8  i  8  8  8  i  «  6  i  8  E  8 

8  3  *  S  »  S  3  M  !i  *  *3  3  '  3  9  0  9  3  9  9  *  9  S  S  S  9  9  9  S  9  S  9  9  9  5  3  rJ  3  9  VJ  5  9  3  9  9  $  3  r  9  9  a  j  3  9  5  ■>  9  9  9  "  •.  v  3  c  9  >j  9 

»  .  f  f  »:  ■)  .  "  <4  i-  «  .  ■  1  -  ’  .i  ?*  ;»  M  *  M  1’  a;  Jj  •-  jVH  i  ll  J1  H  *:  4>  4>  *•  «s«*  *  4‘  «  Sj  v  Vij  >«•  1  ‘M  «  «■  •  t  f.  k>  v,  [  •  .t  >  -  -  K  "  4  )j  4  ;i  «  ••  q 

->  iV4l 


The  computer  output  is  presented  in  Table  IV. 
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